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PREFACE 


The operation of many industrial devices, especially aircraft engines, is based on the 
laws governing the development of turbulent jets. 


The use of high-velocity and high-temperature gas flows for aerodynamic research into 
flow around bodies makes it possible for us to arrive at recommendations for the choice of 
heat-insulating materials to be used in aircraft. Therefore, the solution of jet problems is of 
great scientific and practical importance. 


The present collection contains the results of work in the area of turbulent-jet gas dynam- 
ics. These results have been obtained by the authors of the papers that follow after the appear- 
ance of G. N. Abramovich's monograph [2]. 


The theory of free turbulent jets is extended by V. A. Golubev to the initial region of a 
very-high-temperature jet (with temperatures as high as 20-10° °K). Since physicochemical 
transformations — dissociation and ionization — occur at these temperatures, theoretical investi- 
gations must be based on special forms of the equation of state in which changes in specific heat 
with temperature are taken into account. 


In the solution of this problem, a modification of the Prandtl theory of turbulence is made, 
since Taylor's theory leads to disagreement with the boundary conditions. The author gives 
results of measurements of the outer boundary and the temperature field in a plasma jet which 
are in satisfactory agreement with the theory. 


V. 1. Bakulev considers a jet of real gas at supercritical pressures, describing the state 
of the gas by equations of the vander Waals type, and solves the resulting problem for an arbi- 
trary value of the turbulent Prandt!] number. The equations describing the propagation of the 
jet are obtained for both the initial and main regions. 


The theory of jets is applied in O. V. Yakovlevskii's paper to the calculation of mixing 
processes occurring in closed channels, as well as to the determination of the shape of an 
isobaric chamber. 


A. N. Sekundov investigates the propagation of jets in an oppoSing stream in both a closed 
channel and in the absence of walls and compares experimental results with theory. 


The paper by Chiang Che-haing is devoted to the study of supersonic jets with underexpan- 
sion. A method for the calculation of the gas-dynamic and main regions of the jet is proposed 
and extensive experimental material is presented for a wide range of variation of the degree of 
nonsimilarity and Mach numbers, which confirms the suitability of this method. 


Integral relations are used in the theoretical work described in the last three papers. 
Finally, the papers of the present collection devoted to the study of jets of real gases are 


supplemented by an appendix written by G. N. Abramovich, V. I. Bakulev, I. S. Makarov, and 


V 


B. G. Khudenko in which results of experiments on submerged jets of real gas are given and 
compared with theory. 


G. N. Abramovich 
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HIGH-TEMPERATURE TURBULENT JETS 


V. A. Golubev 


The continuing progress of aviation technology has created the necessity for a wide range 
of investigations into high-velocity and high-temperature gas flows. Of particular interest is 
the study of the propagation of a free plasma jet. 


A large number of papers by Soviet and foreign authors has been devoted to problems as- 
sociated with the study of jets of incompressible liquids, as well as of gas jets at moderate tem- 
peratures [1-8, 10-16, 18-20, 24, 26, 30-33]. The general formulation of the problem of the inves- 
tigation of high-temperature jets was given in [2} and further developed by other authors [13, 14, 
16, 21, 35]. In some of these papers [13, 14, 36], instead of the usual equation of state for iso- 
baric flow, the authors have used an approximate equation connecting the density and enthalpy, 
namely, p=A/i"™, which considerably simplifies the calculations and at the same time makes it 
possible to include physico-chemical transformations (dissociation and ionization) that take 
place in high-temperature gas. In this case, the gas cannot be treated as a one-component fluid 
having a constant specific heat because the reactions occurring in it produce a change in the 
molecular weight and gas constant [28, 29]. Therefore, temperature can no longer be used as a 
characteristic of the gas determining its energy level. Enthalpy which includes both the ther- 
mal energy and the energy of formation of the products of a reaction (dissociation or ionization) 
should now be used as a characteristic of this type. 


The enthalpy, density, and other parameters of the working substance at high temperatures 
were determined by means of a thermodynamic calculation. It was assumed in this that the pro- 
cesses of dissociation, ionization, recombination, etc. taking place in the gas proceed at thermo- 
dynamic equilibrium. The composition of the water-vapor plasma was calculated with the cor- 
rections for the Coulomb interaction between the ionized particles and electrons taken into 
account (i.e., the Debye-Hiickel correction was introduced). 


The thermodynamic calculations were used, in particular, to obtain the variations of 
enthalpy and density with temperature [23], so that it was possible to determine the connection 
between density and enthalpy. This relation for water plasma and air plasmawas approximated 
by the analytical expression a= A/i® (where A and n are coefficients which depend on the type 
of working fluid and which are constant within the temperature interval 300 to 2-10*°K). The 
graphs of the function log p=f(log i) are shown in Fig. 1. 


The maximum deviation of the approximate relation from the data of the thermodynamic 
calculation for water plasma occurs at a temperature of 4000° and is less than 30%, while there 
is an appreciable deviation from the approximating curve in the interval 2000 to 6000°K. Above 
6000°K, the agreement between the approximate relation and the thermodynamic calculation is 
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satisfactory. In the case of air plasma, this 
error is considerably smaller than that for water 
plasma over the whole range of temperature 
under consideration. 


1. Formulation of the Problem 
and Basic Assumptions 


In addition to being of practical interest, 
the study of very-high-temperature jets is also 
of theoretical value since an attempt is made 
here to study the behavior of free turbulence in 
an ionized gas. As arule, the theoretical in- 


tt ff. 8 
/ xX Thermodynamic 


em Ne calculation [23 Seer , 
ee Ae ay eae it vestigation is based on the solution of the usual 
e e e e 
- ee ! equations of the mechanics of a continuous medi- 
ee | wy 
5 ; ; = Tog i um consisting of a one-component gas. How- 


ever, as the result of dissociation and ionization 
processes, the gas becomes a multicomponent 
medium. Therefore, strictly speaking, we should 
also consider the equations of self- and thermal 
diffusion since the existence of concentration 

and temperature gradients leads to the occur- 
rence of mass and heat diffusion fluxes. 


The variations in the composition of the 
working fluid with temperature are taken into 
account by means of a new equation of state. 

On the basis of this, the problem of a jet of 
ionized gas can be solved with the use of the 

| usual equations for a boundary layer neglecting 
© log the equations of diffusion and thermal diffusion, 
Fig. 1. The variation of density with en- provided the following assumptions are valid: 
thalpy: a) air p=104-°°/;9-821. by) water 


vapor p =105+13 749-91, 1. The processes of dissociation and ioni- 


zation are taken into account through the use of 
anew equation of state (9 = A/in); 


2. It is assumed that the jet has the same physical properties as the surrounding medium. 
The justification for this assumption is the fact that the working fluid (water-vapor plasma) and 
the surrounding medium (air) have almost the same values of the exponent n characterizing 
their physical properties (see Fig. 1). 


3. Assuming that the jet is completely turbulent, we neglect in the energy equation terms 
that describe the energy of radiation and in the equations of motion and energy terms that de- 
scribe molecular viscosity and heat conduction because the latter terms are small by compari- 
son With the turbulence terms. 


2. The Derivation of the Basic Equations 


The case under consideration is of a planar high-temperature gas stream which moves 
parallel to the x-axis. Figure 2 gives a schematic diagram of the jet. At the point 0, the gas 
begins to mix with the surrounding medium. Let us suppose that the parameters of the emerg- 
ing gas (subscript 0) and the parameters of the surrounding medium (subscript H) are given. 
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Then, the following equations can be written down for the plane-parallel jet of the nonstationary 
gas stream [2]: 


equation of motion 


equation of continuity 


energy equation 


0 . ii? 4- UP O ; tut ye 0 Lo gk bel 
0 {+ ———— ol) —- {f }--——--—-J+su —fi-+ = 
"Ot ( T 2 }e OX | T 2 -_ dy ( | 2 
Op 7) 7) OFe Jq y 
= — Uz Ut us + Ut —— 
ot + O% (59 yt dy ( vy F yy) 1 ox | ay (3) 


When fully turbulent flow is being considered, molecular transport terms can be neglected 
for an isobaric jet of compressible gas and Eqs. (1)-(3) for the instantaneous quantities can be 
rewritten as follows: 


Ou Ou Ou 
9 — bon -4U = (0; 4 
‘Ot " Ox am Oy ) 
Oe , oO . O 
—_—— -- Gil) - t —(; 
52 (¢4 8 )4 (i 4 Je t += )=0 (6) 
Ot 2 Ox a 2 


In addition, in Eq. (6) we have neglected terms describing the kinetic energy of the transverse- 
velocity component since it is considerably smaller than the longitudinal component. On the 
basis of the equation of continuity, Eqs. (5) and (6) can be written as 


0 7) 0 
—— (pu) -- —— (su-u)+ —— (sv-u)=0; 7 
sr MAG (tet) = (ya) (7) 
(op 4+ Fe) $2 (i-v) + 5 (et +S (0 <0 (8) 
a \ 2 on ae dy * * Ox ( 2 gy \ 2 : 


The solution of the jet problem can be obtained 
with the help of either the Taylor model of a turbulence 
or the Prandtl model. 


Stream 


y 
bos Ip, Qo: Uy *, (iz), ° . ops 
The choice of the Taylor model can be justified 


by the fact that the velocity and temperature profiles 
for an incompressible liquid obtained from it are in 
good agreement with experiment. On the other hand, 


tn Tas Qui when the Prandtl model was used to calculate such a 
uy 4 jet, the theoretical temperature profile did not agree 
Surrounding medium e-(4} with experiment. However, as has been shown in [13], 
the Taylor theory leads to contradictory results when 
Fig. 2. Schematic diagram of the it is applied to compressible fluids. The turbulentfriction 


initial section of the jet. stress onone boundary of the jet is not cqual to zero. 
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On the basis of this result, we can conclude that the use of the Taylor theory of turbulence 
becomes incorrect for the description of jets of compressible gases. On the other hand, the 
Prandtl theory for incompressible fluids can be modified in such a manner that it leads to the 
same theoretical results as the Taylor theory; moreover, as will be shown below, in the case of 
jets of compressible gases the modified Prandtl theory [14] does not have the drawbacks that 
are assaciated with the use of the Taylor model. 


Indeed, if we use the Prandtl turbulence model for an incompressible-fluid jet, we obtain 
the following formulas for the turbulent friction stress ry and the heat transfer wr: 


Ou af ou \2 
Tq SE =nli, (| ’ (9) 
Oy oy 
WW, _ dh, OT —¢ “ alt Ou OT (1 0) 
Oy / Oy oy 


whereas if we use the Taylor model, the expressions for t- and wp are 


nahh (EY (11) 


we co. 2 OF (12) 


where ly, /7, and 1, are the mixing lengths for momentum, heat, and vorticity transfer, respec- 
tively, while ¢ and Ay are the coefficients of turbulent exchange and turbulent heat conduction, 
respectively. 


It has been shown in [20] that whenever the Prandtl model has been used in investigation 
devoted to the study of this problem, the mechanisms of heat and momentum transport were 
identified, i.e., it was assumed that J; =1,, and whenever the Taylor model has been used, the 
mechanisms of vorticity and heat transport were identified, i.e., it was assumed that lp =]. 


Moreover, in order that the values of the experimental constants be the same in both 
models of turbulence, it is usually assumed that 1,=V2ly. Then , formulas (9) and (11) give 
identical values for the friction stress, whereas formulas (10) and (12) give different values for 
the heat flux. Because of this, the velocity profiles obtained on the basis of the two models coin- 
cide, whereas the temperature profiles are different. As was pointed out above, the temperature 
profiles calculated on the basis of the Prandtl theory contradict experimental data, while those 
calculated on the basis of the Taylor theory are in good agreement with experiment. For this 
reason, the Taylor theory has always been favored. In fact, any attempt at a solution of the 
problem should be based on the following considerations: 


1) it follows from formulas (9) and (10) that the turbulent Prandtl number is 


Pr_=—-——_ => (13) 


from which we have 


p® ‘ 
= — =72 Py,; (14) 
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3) it has been found on the basis of Reichardt's experiments that ly/lp=0.5, from 
which We have Pry =0.5. 


Therefore, if we take the above value of the Prandtl number, formulas (13) and (14) con- 
firm the hypothesis that the mechanisms of heat and vorticity transport are the same, i.e., that 
I~ =l,, Whereas the hypothesis that the mechanisms of heat and momentum transport are the 
same are not confirmed since in this case /7T=2l,. 


From what has been said above, it follows that the solution of the problem significantly de- 
pends on the value adopted for the turbulent Prandtl number. If, in solving the problem accord- 
ing to the Prandtl] model we take Pry =1 and, consequently, /y=/T, then the u and T profiles 
coincide,which is correct for the assumptions used in solving the problem. 


The same result is obtained if in using the Taylor model, we take according to (14) 1,/lp =2 
corresponding to Pry =1 and not /,/lT =1 as has been done by all previous authors. Consequent- 
ly, We can obtain a correct solution of the problem of an incompressible-fluid jet on the basis of 
both models of turbulence if we take PrT=0.5. For consistency, when solving the problem ac- 
cording to the Prandtl model we must take the mixing length for heat transport to be twice as 
large as that for momentum transport, whereas when we use the Taylor model, we must take 
the mixing length for heat transport to be equal to that for vorticity transport. 


Proceeding to the case of a compressible gas, we should apply the Prandtl theory because 
calculations carried out for a compressible gas on the basis of the Taylor model have shown 
that the solutions so obtained do not satisfy the boundary conditions: The turbulent friction 
stress must be zero on one of the boundaries of the jet, but this cannot be realized. 


Further, if in using the Taylor model of turbulence we assume that the turbulent friction 
stress is zero on both boundaries, the problem cannot be solved at all. 


A modification of the Prandtl theory will be made below. Using this model, we will solve 
the problem on the basis of the same assumptions that are used in conjuction with the Taylor 
model. In order to do this, we follow Reynolds and replace the instantaneous values of all 
quantities in Eqs. (5), (7), and (8) by the sums of the time-averaged and fluctuating components. 
Here, in addition to introducing the fluctuations of velocity and density, we use van Dreist's as- 
sumption [9] that the current density (pu and pv) and enthalpy (i) fluctuate as a single quantity. 


After performing the averaging, we obtain the following system of equations: 


0 0 


— 1a — ar nT pene Ee Pe 
(ott + 4’) (ott tt (out) wu’) + “— (60-u + (ov)'n’') =0, (15) 
ot Ox Oy 
M4 € (5) + —_— (pv) =0, (16) 
ot Ox 
17 — 1 ial Oo (7 F1 Gani ~3 0) 77 Ty 
° (ni + 9'i’) + ° (pei (pu) i’) + (pu-f + (0) i’) _— d (ot? + ou!” + Quon’ +o wi’? ) 
ot Ox oy 2 at 
| ) aa ane at gel Vn t / ~~. 29 4 7. 7 / 
y 


Differentiating Eq. (15) and making use of the equation of continuity, we can reduce it to the 
following form: 


— Ou — On — Ou 0 1s QO 7, 
6 — +ou ——-bLeov of i!) — —— (sv )'u'. 18 
ary +f Ox +i ( 3 (p ) ( ) 
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Differentiating expression (17) and transforming it with the help of the equations of con- 
tinuity (16) and motion (15), we obtainT 


— i ! —— —-—- (av) — (sv) 19 
dy ys (19) 


In Eqs. (18) and (19) we have neglected the partial derivatives of the fluctuation terms 
Since these are small by comparison with the derivatives of the fluctuations with respect to y 
and we have also neglected derivatives of 0 2/2as being small by comparison with 0°/2 andi. As 
a result of these transformations, Eqs. (18), (16), and (19) for the averaged stationary motion 
become 


— ou ,— gu 0 TTT 
all Ou Vy of (ov)'IL", (20) 
Ox Oy oy 
(710) +-— (0) =0, (21) 
4 
— dO — di 0 cp a! Ou 
tl UC —_ = — 10) t — (ov . 22 
m+ yr OYE — (wou! “> (22) 


Transforming the right-hand sides of Eqs. (20) and (22) and taking into account that vis a 
smal] quantity, we have 


-_ = fa) --——_, 
—~ 2 (oy = 2 =~ Gf) $0) u' = ~—— (son), 
oy Oy Ov 0 
Ope o —_, —~W7 
— —— (pvy t= — (6+e')(u+o')i = — (pv't’). 
Oy oy 
Equations (20) and (22) can now be written as 
a Ou tou Ou _ _¢ (5 vit’) 93 
\ Ox ‘ oy Oy ‘ J ( ) 
_— } —. 
(ou) + — (pv) =0, (24) 
dy 
— — dO > Ou 
ll (i U ‘1 —- 4. Y il 
t Ox TP Oy < — le vt) Oy ‘ (25) 
where _ =~ ~< (s-v’u’) is the derivative of the turbulent friction stress, 
dy y 


(,-v’i’) is the derivative of the turbulent heat flux, 


ee, "1 1 dua . . ° ° oe 
— pu’ - ee is the dissipative term arising from turbulent stresses. 
vy 


In the same way as in the case of incompressible fluids, we assume that in Eqs. (23) and 
(25) we have 


Ou and i~ 1, — , 
Oy 


—v~u'~l, 
where, on the basis of what has been said above, the mixing lengths 1, for momentum transport 
and Jj for heat transport have been taken as ly =cx, 1j=2cx and asconstant over the cross Section 


of the jet. 


T These transformations are given in greater detail in [9]. 
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Equations (23) and (25) can now be written as 
a Ou ” Ou eye 2 [7 Ou \2 
Ul —* UV See“ yu {--— ; 
ar re Ov vy : ay | (26) 
0 —_—, | 0 —_— 27) 
~~~ (pu) + —— (ov) —0, ( 
ON Oy 
— gs — Yb “~ Oy ii sf du 3 
iy ap energy? a | cate vu j. (28) 
Ox Oy Uy Oy oy dy 


In order to have a closed system of equations, we uSe the relation connecting density and 
enthalpy 


(29) 


3. Transformations of the Basic Equations 


In order to solve the system of equations (26)-(29), we will introduce the new variables x 
and y=y/ax and dimensionless quantities of the form 


— — 


F'(2)= ~ , 9@)=+-, and -()=—, (30) 


roy by fy 


Where a iS an experimental constant. 


Using the current function p= | on -dy==p,uaxF (2) and taking into account that 


Ox ye OV ~ ax” (31) 


we obtain the following expressions for the current density components: 


~ y . 
ot 2 ot, F’ (2), 
Uy 
(32) 
— Oy . 
yu == — SE == Ayylt, [ZF (2) — F &)I. 
Since 
0 iy 0 0 da 0 
—— — uf } _-—_ = naan ’ (33) 
Ox Oz ON vy Op OY 
we use the equality 
— -- ~y -- Gp Ou => > Oo Ou 
on =o 0 = 5u + 11, “he pn + 2c? ; (34) 
dy oy dy ay 
and the relation 
Yer = aXx’, (35) 


which also holds as in the case of an incompressible fluid, to obtain the longitudinal velocity 
component 


en ee —<.(+}. (36) 
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Making use of the continuity equation 


Jen, (37) 
we can obtain an approximate form for Eq. (36) valid when a «1T 
~~ PF’, (38) 
Uy 
With relations (31), (82), (83), and (35) taken into account, Eqs. (26) and (28) become 
a 1 0 ds liad ye, ¢g au aru | __ 
yl yP . TD 7 ) -F oe dp? J=o. (39) 
a “ ? ? ? + 
ptt, F tig Ae ae a ae .( du y=0. (40) 
| do do dg ay d3- do dg dag- 2 ao 
For formulas (38) and (30), we have 
anus" [F’ 4. nF’ —| | 
de = 
Ou Ue" Jet nb” ane (1° —n) fF’ (—)} + nh’ — ; (41) 
doe 5 t 7 
i __ Q/ 4 Tt di taf d-i —j 7 
de =f) 7 = 70) A+ ) de — "Oo" » do" O°: 
Substituting these expressions into Eqs. (89) and (40), we obtain after some transformation 
prs Bap 24 (Ban) Fr (LY enh EP =0 (42) 
Ww no noo / 7’ \2 , 0 > us l d / 3 
Fd nF" 24 FY 2 nF (= 1F 2-4 F 42 — | oF ") —0. (43) 
: </ as < 2iy wh dy 


For subsonic jets with high initial preheating, i.e., jets for which Uy” << 2i), the last term in 
Eq. (43) can be neglected. Subtracting Eq. (48) without its last term from Eq. (42) and trans- 
forming the expression obtained, we have 


t 


= Hn (FY pn’ <~)-<- (Pi baF —}=0 
< ad * 
or 
2 5 a T 


(+ Pan -=) [Pepa —)=0. (44) 


An examination of Eq. (44) shows that if try to satisfy this equation by taking the second 
factor equal to zero, we will obtain a trivial solution which is physically meaningless. Indeed, 
the second factor of Eq. (44) is identical to the expression for the velocity gradient du/dg which 
cannot be equal to zero over the whole of the cross section of the jet. 


T The exact solution is given in Section 6. 
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Therefore, let us equate to zero the first factor of Eq. (44), namely, 


~ - 


Nn —_—-— —(), (45) 


| 
2 


On the basis of the energy equation (45) and (42) or (43), we obtain the equation of motion 


FUL nk! St (np) F! [J + F=0. (46) 


Setting p=const. and taking the equation 


into account, we will reduce Eqs. (45) and (46) to those for an incompressible fluid 


vf 


7" =-Q, (47) 
PY 4 F—0. (48) 
From the equation 

do Ou 


(VU = 40 -- pV! = pu + 2¢7x°—— 
dy ay 


we will obtain the expression for the transverse velocity component 


ds dau 4 

_ _ ax? —— —— —— 

vu pv dp ao arx- 

a 

Hy Olly oll, 

or, after some transformations, 
= eh’ —F) +n < (FY par’ —)|. (49) 
ally . “ 


The turbulent friction stress across the jet cross-section can be obtained from Eq. (26), 


Oz; 9.) 6) ~~ ( Ou y 
—*_ = cx? —_| p ; 
dy oy oy 


which on the basis of Eqs. (29), (30), (31), (33), (37), and (41) becomes after an integration 


ian je" (9) nF’ (| , (50) 


Dll 
2 


Q 


Here, the integration constant is zero if the integration is started on the boundary where Tt7=0 
and au/dy = 0. 


The turbulent heat flux in the transverse direction with the frictional heat neglected is 
given by Eq. (28), namely, 


10 V. A. GOLUBEV 


We can obtain the equation of heat transfer by performing analogous transformations and an 
integration, 
W 


AI oLolo 


=| F'@)+nF (a) 


~?f 
v 
- 
. 


|. (51) 


In the case of constant density (2 =1), the expressions for the longitudinal and transverse 
velocity components and for the turbulent friction stress and heat transfer reduce to the well- 
known [1] expressions for an incompressible fluid, 


—=F' (2), (52) 
Uy 
—— =¢F’ ()— F(s), (53) 
ally 
= 9 =([F')P, (54) 
Pile) 
a we 
2? 
Wr _2/F'(¢). (55) 
AGllolo 


By solving Eqs. (45) and (46), we can obtain the velocity and enthalpy profiles for a plane- 
parallel boundary layer of a turbulent gas jet. 


4. Solutions of the Energy Equation 
Let us take 


Z= ; (5 6) 


Consequently, taking 


we find after an integration that 


Qn 


—2* =c+e. (58) 


2—An 


9 


The boundary conditions for the inner boundary of the boundary layer are i=ij, i.e., 


t=1 when g=q), (59) 


and for the outer boundary they are i=iy, i.e., 


T=Tr When g~=@o. (60) 
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On the basis of conditions (59) and (60), we have 


2 [— 
Cp-= ) (61) 
2—n f1— 2 
2 l—<,° 
2— A Fl — F2 
Substituting these expressions into Eq. (58), we obtain the expression for the enthalpy pro- 
file 
[poo 6s 
where 
en 
¥1 — $2 
With n=2 we have 
—=¢1, (64) 


and after an integration We obtain 
Int=cyoten. 


With the help of the boundary conditions (59) and (60) we will find the integration constants 


C= — tw ; (65) 
1 — Fe 
c= Hg, (66) 
71 2 
and the function 
ct”. (67) 


With n=0, Eq. (45) reduces to the equation for an incompressible fluid and its solution 
then has the form 
t= + (1L—n) tr. (68) 


5. Solution of the Equation of Motion 


Let us study several particular solutions of the equation of motion as we were unable to 
find a general solution (for arbitrary n). 


1. With n=2, Eq. (46) becomes 


F" (@)-3F" (2) —+2F" (2) (-) + F @)=0. (69) 


om 
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Introducing the new variable 


4 — + F2 
: Th — 
we obtain in place of (69) 
F" (y) 4 3t-F" (1) 420 F'(n) +B F (q)=0, (70) 
where 
Ly 
, In = In — } 
—__ Jee = and b=% — Ga. 
+ 71 — $2 TL — Fs # 


The characteristic equation of the equation of motion 
R341 3th? -+9Pk+b,=0 


can be reduced by the substitution z=k+t to the form 


z*+dpze+2q=0, (71) 
where 


The roots of Eq. (71), are 


where 


It should be noted that the solution of Eq. (70) is meaningful only when D > 0, where D=¢ + p. 


Setting 
a,=-2,—f, y= —— 21, B= v3 (§ — 4), 
we find that the roots of the characteristic equation are 
kya, Ro=a,+13, &, =a, —i3. (72) 


The solution of the differential equation (70) corresponding to the roots (72) is of the form 


F (4) = )e" + Cy" COS By + €3€7" Sin By. (73) 


The integration constants c,, c., and cs, as well as the inner and outer ordinates, g, and 
@, of the boundary layer of the jet are found from the same five boundary conditions as in the 
case of an incompressible fluid. 
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On the inner boundary of the jet (where g= 9, or 7 =1) we have 


U=Uy, 9=6,Le., F’(2)=1 or F’(1)=),; 
f 


7=0, Le. FY (a)= —2—— or F"(l)= — 2b,¢ 
} 


v=0, Le., F(z,)=2, or F(1)=*,. 
On the outer boundary of the jet (where y=q or 7 =0) we have (74) 
u=O0, p=, l.e., F’(9,)==-0 or F’(0)=0; 


Ou 
OY 


=0,i.e., F’(2)=0 or F’(0)=0. 


On the basis of these boundary conditions and Eq. (73), we can write 


FY A) = eye 94 4 C90" YE Oye” f= O,, 
PF" (1) 0," af C9 (4271 — Pa) + Cg (AT. + 811) = — 20, t, 
FP (A)=cye +c ,e% cos p+ ce” sin 3=~, (75) 
FY? (O) = 012, + Cot. C33 =0, 
F" (0)= e024 cy (08 — F) +.0,2228 =0, 
where 


“,s=2,c0s 8 —B sin §, 


Yo=a,sin 8-+-8 cos 8. 


It should be noted that the solution of the system of transcendental equations (75) is as- 
sociated with a number of difficulties since the coefficients a,, a2, 8, etc. are themselves func- 
tions of the initial heating ty and the width b, of the jet. Therefore, we will propose the follow- 
ing method of solution: The fourth equation of system (75) yields 

As rc 


—c,—-e. 
ay ay 


Substituting this expression into the fifth equation of the system, we obtain 


Bt 


a; —Q2a | 
a C38 d 2 and Co = C3-—— 
ay 


ae — @}a, — 8? as —a,%, — $2 


Next, after the substitution of the expressions for c, and c, into the second equation of system 
(75) and some transformations, We find that 


26,¢ (a5 — ay — f°) 
qa ee, (76) 
(25 + 6°) M 


after which we can determine c, and cy 


_ 26,3 (a, — 2a) 
(as + 5°) M 
26,3 


C= — —* 78) 
a, M ( 


(77) 
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where 


Masafie {| e™ (ay, 472 fly). 


Substituting the values of e;, cy, and cy into the first equation of system (75), we obtain the 
formula 


on ec 6 ‘2 od 
ep Pe remy Mr Ya) Bt] Ga yay BY gL Pad, (79) 
* | 
from which we find the solution (i.c., the value of b,) of this system. In general, the determina- 
tion of €;, ey, C3, gy, and g, ean be carrted out in the following manner. 


For a given initial heating of the jet my or t In(ip/iy]), we take several values of b, = gy, — go 
and find the values of £, ¢, ay, Go, 2, V4. yo, ete. from the formulas given above for each of the 
values of b, . Substituting these values into formula (79), we find the value of b, which satis- 
fies Eq. (79). Next, using the value of by found in this way we reealculate the values of &, ¢, 
ete. and then use formulas (76), (77), and (78) to find the values of cy, ey, ¢3. ‘The inner boun- 
dary of the jet is determined from the third equation of system (75), 

2hyt [3 


C\ — | 


y : 13 (a, — Zaz) cos B -f (a2 ~ aya — fl) stn 8 |. (80) 
i 7 1 | 42 


| 


Substituting the values of the constants c,, cy, ele. into tq. (70), we obtain an expression from 
Whieh we ean find the funetion (7) and its derivatives. 


We ean determine the profile of the longitudinal velocity component in the boundary layer 
of a high-temperature jet by means of formula (38), 


J Pi (e) <5, (81) 


the transverse velocity component by means of Kq. (49), 


U ) ( ov . re f nd ye? f ‘ 
ee ° — | (« ) to) -}~ i} ro) — 
= leh) — FOAL [P+ 27 7 || (82) 
and the velocity gradient by means of formula (27), 
{ lig ") </f . f 
— ee | OPM (wy © s: 
mM, dg \! (9) | 20" (2) b. (83) 


In order to determine the apparent turbulent friction stress, we make use of Eq. (50) 
Which can be written as 


t af jon opr, LF | 
as G (~) -f- 24°" (9) | . (841) 
folly by. 

Qu 


Cd 


The calculations performed with the help of these formulas are given in Tables! and2. These 
tables show that the turbulent frietion stress on both boundaries of the jets is zero. 


According to the Taylor model of turbulenee, the apparent turbulent frietion stress in the 
boundary layer of a jet is given by the expression 


Os oie ON Ol 
Ts. D2 yal (85) 


Oy oy? 
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TABLE J 
f) 2: hy 1.98 « cal 4.3. aT ] wie 10 on Cy 3.420. Cy $355. Cs AQ LAS 


| 7 ) oda U m1 W/ 

f # * ry) ry) ‘¢) mp Hy iy dy (thy i lly atattols 
2 
an 2.60 1,757.1) 3 35.98 | 0.0 0.0 O6 97 0.0 QO,0 I.tb3-10 4 0.0 0,0) 
i” 2.052 | 3,31-10 4 36,09 | 3,58 BOOK 8 0,079 | 3.93910 5 TK 10 4 GO RU-10 4° 286-10 4 0.0566 
. 1. WA ),0063 $2.02 | 8.23 6.33 O.18%8 | 328-10 1 9.5-10 4 O81 10 5 002246 O15) 
al. (), 856 Q.0118 - 20.7 O.9 2.99 6.1% 153-10 % 60,0036 0,0] 0.0922 0322 
0.4, ),25% 0.0224 IQ.88 11.42 0.45260 V4N7 | 5.73-10 4 0,019 0.0335 983 0.5628 
0.5) 0.34 0,442 13.29 10.48 | — 2.485 2.48 O,O18% | 0,048 0,103 1}, 685 , 9,875 
0,6 Y.94 (079-4 - 7,05 | 6,06 3.448 Q.912 0,045 0,094 0, 290% } 392 | ) 25 
0.7, 1.53 0, 1495 3.0 | 6,48 3.7 Q,1216) oO. 1445 0,224 0,76 2,25 1 O84 
Q,4} 2.13 (),282 0,295 | 4.34 | - 3.3 779 (0345 ().47] 1.7 2.76 1,76 
gg! 2,73 0,55 2 ans) 2.47 ~ 2.43 1.13 0691 0.074 2.67 1.613 1,348 
10) 4.33 1.0 3.53 1 1.0 | 2.42 11g) AD) 0.0 0 0.0 0,0 
| 


TABLE 2 
N- 2; b, - 4,986. ¥\ 2.78; Ty 2 J 049-1072, Cy =z—] 2. Cy~—11,85, cya —I0,77 


Tt 
ray ery | ery |e |e fe 8D gag | 
thy) iy dy aity i AQ lyl y 
O.0 |- 2.206 | 0.01049 | - 13.57 | 0.0 a) 13,4 0,0 0.0) 1,475.10 4 0.0) 0.0 
0,1 —1.707 | 0.0166 | —13.18 1.057 5,08 2,14 2.91-10 4] 00,0015 5.81 8] 7.77-10 5 (),. 0803 
0,2 | —1.209 | 0.0261 -J2,2] 2.76 3.354 | —1,63 |) 1.885-10 47) 00058 (), Q165 (),0478 (),2 
0,3 | —0,7) ,04)] --J0.43 4.185 2.21 -— 2.63 7.07 10 #] 0.0166 (),.043 164 0), 3685 
4 —0.21) | 0,0649 8.44 4.96 0,885 | —2.59 (),0209 (), 0418 0, 1062 0,416 (),588 
8 0.249 | 0.1022 - 5.604, 4.34 (0,299 | - 2.9] 0.05575 | 0,099 0) 251 0,932 (), 879 
0.6 (787 | O.1G15 —5 VG 4.7 —1.184 1.474 (), 1228 0,493 0,53 1.422 1 O87 
),7 1.284 | 0,255 — (955) 3.92 — 1,705 () 624 0,254 0,551 1.028 1.884 1252 
0,8 1.784 | 0,403 0,767) 3.0 - 2.122 (),249 (),482 0), 90 1.75 1,914 ] 258 
0,9 2.284 | 0,644 2.012) 1,99 2.4025 0,213 (), 80) () 047 2.2 1.034 (),928 


1.0 2.78 1.0 2,77 1.0 1,43 1), 548 1.0 0.0 4) 0.0 0.0 
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which after some transformations can be written as follows for the case n=1: 


—t=-2\re[Perre 
Gully ° 
a —- . 

2 


\de (36) 
b, 

Taking the lower limit of integration to be the outer boundary of the jet g, where 8u/dy = 0, 
T= 0, and introducing the variable 7, we can rewrite Eq. (96) as 


- yp 
°T —_ 


= 2 (IFO OD EP (4 an, (87) 
Pylo 


2 


Q 


The expression for my" can be rewritten as 


ot yin (= = ,67" In n=", 
Substituting this expression, as well as the equations for the function F(7) and its deriva- 


tives into Eq. (87) [the function F(n) with n=1 for the Taylor turbulence model is given in [18], 
we obtain after the necessary transformations 


“T 


cd ° 2 1 ‘ r, 
= im, [eh dri} [ep (m, +1) +e] | e*" cos Bud + 
0 


2 
2 


ui 


+ Cy [C3 (77, + mz) — Comms] \ e"" sin Budy “f Co (CoMte + Cys 


0 


) ( e*" cos? Brady + [C3 (Cont, + C33) + 
0 


7 ui 
+ C2 (C3ilx — Com )] \ e**" sin By cos Brady + Cy (CyM— Cos) | e*" sin? Bad, (88) 
0 0 


where 
M, =2, (%,-+ 4), ki = 22, +1, 


Mz= a3 —B? + ant, ko=4, + %-+#, 
ms—=f (2a.+ 2), ky=2a,+1. 


The expressions for c,, Co, cz for n=1 are given in [13]. If in Eq. (88) we perform the 
transformations 


cos? p43 = — (1 + cos 287), 
sin By, cos By = — sin 28%, 


sin? =— (1 — cos 23%), 


then Eq. (88) can be integrated to give 


o~- 


“T 


2: Ray | okey 
alle b, ky Ry Ro -1- 32 ke 4. 


’ i 
a 
2 


- 


eben 3 Ca (Cutts +- Cgity) fF e®% | 
+ Cy [Cx (7, 4 my) — C2m3] | ————— (Asin 34, — 3 cos By)-+ $a 
kAyt+ © ; + 3? 2 k; ky 
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TABLE 3 
ee oe n=l] . a 
y= 1.75710 8 | tO fy T.230 y= 1049-10 2 
" 1 du 1 du tt 1 du tT __ ] du TT 
wile | | ea | ae | ae |g a 
2 ) 9 
0.0 0.0 0.0 0.0 0.0 0.0 O.Q 0.0 0.0 
0.1 0.01136 O.0488 0.01415 0.0478 (), 0226 (). 0-467 0.0249 0.0424 
0.2 | 0.0329 0.2367 039 0.2237 0.059 0).236 0.063 0.1955 
0.3 | 0.0668 0,626 0.0794 0.6017 0.1115 0.58 0.1183 0.4975 
0.4 | 0.1223 1.25 0.1406 (), 2067 0), 1896 1,123 0.1935 0.94 
0.5 | 0.201 2.08 ().2265 1.9697 (0.2865 1,81 ().286 1,48 
0.6 | 0.304 2.972 0.3306 2.7547 0,394 2.49 0.366 2.001 
0.7 | 0.419 3.75 0,438 3.3697 ), 488 2,69 0.49 2,36 
0.8 | 0.48 4.12 0),.493 3.6057 0.505 3.17 0.503 2.55 
0.9 | 0.405 3.892 04035 3.4357 0.424 3,01 0.495 2.415 
1.0 0.0 3.695 0.0 3,251 0,0 2.83 0.0 | 2.328 
] | 
+ 2k, cos 234-+ 23 sin 237) — eee 284, — 
k3 + 442 kt Ae? 2 : k3 + 4 
hy kyr, 
es Le ee ee eee eer | 


Table 3 contains the results of the calculations of the velocity gradient and turbulent fric- 
tion stress carried out according to the Taylor turbulence model. It can be seen from the table 
that despite the fact that the velocity gradient is zero on both boundaries of the jet, the turbulent 
friction stress on one of the boundaries is nonzero. Similar results have been obtained with 
other values of n and, in particular, with n=2/3 and n=1/3. This result contradicts the boun- 
dary conditions of the problem. 


In order to obtain the solutions of the equation of motion (46) for other values of n, let us 


introduce a new independent variable 7 and perform the necessary transformations. On the 
basis of Eq. (57), we have 
az ' nQ 
=—C\~ , 
ae 
so that 
F' (< —7F PF de ia nips (2), 
ay dz as 
v/,, aF- _- 2 and t (- nn Fal! 
Pi (a= a (AA) [bP +k], (90) 
3 3 no 3 2 3, 
F" (9)=c [set F’ ()-bS ne? UP (ye Fe «| 
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Substituting expression (90) into Eq. (46), we have 


Son y . “n= | ., 3 3 Sno? 
0 PY (e)+ Bane FP” (2) as nN (— — | 7 F" (*)-+-cF (<)=0, (91) 
where 
oan 
C= ) = state 
C1 2—n b, 


Wewill seek a solution of Eq. (91) in the form of a generalized series 


PF G=>» a,v** (92) 


k-0 


An analysis of Eq. (91) shows that the solution of this equation with the help of series (92) 
can only be obtained when n=2/3 or 1. 


Let us investigate some particular cases of the solution of Eq. (91). 


2. With n=2/3, Eq. (91) becomes 
sh OF" + cF=0. (93) 
On the basis of series (92), we have 


F’ (= Say (rare, 
FY’ (s)= NV a, (rth) (r +R—-1) 748-2, 
roar n ( + )( + ) (94) 


FY @= a,(r tk) (r $R—1) (r+ —2)vte-3, 


k=0 


Substituting this expression, as well as series (92) into Eq. (93), we have 


bM: 


a,(r +k) (r+k—1) reetey Qa,wtk =, (95) 
k=0 


a 
II 
cS 


Equality (95) must be satisfied identically for any 7, so that all of the coefficients of the series 
on the left-hand side of this equation must be equal to zero. In order to determine the constants 
&), 44, a, etc., we will equate to zcro coefficients of like powers of 7, 


-* with k=0, apr? (r—1)=0, 

7"! with k=1, a,(r +1)r=0, 

* with k=2, a,(r+2)?(r+1)+ca=0, (96) 
w+! with k=3, a,(r-+3)(r4-2)+ca,—0, 

wt? with k=4, a,(r+4)(r+3)+ca,=0 


and so on. 


From the first equation of system (96) on the assumption that a, = 0, we find the roots of 
Eq. (99), 
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r, =O, ro=1, r,—0, 


It can be Seen from system (96) that in order for the system to have a unique solution, all 
of the odd coefficients (a,, a3, etc.) must be set equal to zero. The even coefficients of the 
series corresponding to these roots are determined by the recurrence relation 


Ao, = Ok (97) 
(r -- 2k)? (r + 2k — 1) 


Which is obtained from system (96). Here, we have k=1, 2, 3,.... 


Having substituted the values of the coefficients when r, =0 into series (92), we obtain the 
first particular solution of Eq. (91) as 


F, (2)=1—0.23c2? +0,00521c?=* —0.000029c*=" + 0.0000000646c?x? ~.... (98) 


Carrying out the analogous procedure for the second root r=1, we obtain the second particular 
solution of Eq. (91), 


F(t) =t —0.0556cr?+- 0,000556c2t5 — 0.00000189c%x? + 0.00000000292c479— . . (99) 


The third particular solution of Eq. (91) cannot be obtained by means of series (92) be- 
cause the first and third solutions obtained with the help of series (92) will be linearly depen- 
dent. We can obtain the third linearly independent particular solution by lowering the order of 
differential equation (93). 


If two particular solutions of a third-order differential equation are known, then we can 
reduce the order of the differential equation to the first and solve the resulting equation by 
simple integration. Todo this, we take 


FP ()=Fi(s)Y, (100) 
Where Y is the new unknown function, so that 
F3=F\Y+F,Y’, 
Fa Fi iY +2F iY’ +P’, (101) 
F3 =F Y+3F1Y’4+3F iY" +F YY”. 
Substituting Eq. (100) and its derivatives (101) into expression (93), we obtain 
oF LY" + (3tF i 42F,) ¥"4 (BtF i +4F 1) VY’ +: + 2Fi+cF,) Y=0. (102) 


Because of Eq. (93), the term in the last set of brackets is equal to zero. Moreover, if in 
Eq. (102) we set 


U= y'=(2). (103) 
F, 


then.the order of Eq. (102) will be reduced by unity. As a result of this, we have 


Fy | a Fy l 
n — 19 — , te —_—_—— —_ =(), 4 
ur+(3 —+ a (3 ftp Lu (104) 
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Let us now Set 


U=U,W, 
U’=U, WLUW’, 
U" =U W+2U, WU, Ww”, 


where W is a new unknown function of U and 
F.\' 
U, = (=| . 
F, 


When we substitute these expressions into Eq. (104), we obtain 


~- 
» 


yyir y' ! r Fy C | "Tt ° ' Fy 7! | / F Fy | r 
UW" el 201-30) a 20) W + {Oi +301 — + 2U1 — +30, = +4U, — —] W=0. 
1 . l . 4 1 


Because of (104), the term in the last set of brackets of Eq. (107) is equal to zero. 


If we now set 


then Eq. (108) yields. 


or 


Integrating this equation, we obtain 


InV=—21InU,—3In F,—2Int, 
which after an exponential transformation reduces to 


| 
2 3 
US Fy3 


V= 


where 


F,\' FoF im FF 
y,-(2) 
1 Fy 


Substituting this into (110), we obtain the equation 


Fy 


yo ss 
(FoF, — FoF) 2 


(105) 


(106) 


(107) 


(108) 


(109) 


(110) 


(111) 
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Further, from the definitions of V and U we have 


U=U,\V 
and 


Fy=F,\U=F,\U,\V. (112) 


Omitting cumbersome transformations which are necessary for the derivation of expres- 
sions for U and V and the subsequent integrations, we obtain the following expression for the 
third particular solution of Eq. (91): 


F's (s)= —In<F, —0.5¢x? +0,0147¢2<4 + 0,0000962c3<* +-0.0000027 1cts*#—.... (113) 


The result of all this is that the general solution of Eq. (93) will be the sum of the three particu- 
lar solutions, 


F(t) = Fy (t) 4 CaP (t) + 3 F 3 (1). (114) 


The constants c;, Cy, C3, aS Well as the jet boundaries g, and g are given by the same boundary 
conditions as thos for the case n=2. 


On the inner boundary of the jet (with @=, or tT = 1) we have 


, | 
w= Uo, =, i.e. P(r )=1 or F'(I)= 
Cy 
Ou . ' 3. Un 
=—(Q, 1.€., F”(o,)=n10 or FF’ (l)=-— } 
oy ’ (21) 1 ( ) 2 C; 


v=0, i.e., F(o,)=2, or F (1)=¢%1- 
le. F (2))=% (1) 1 (115) 


On the outer boundary of the jet (with g@=q@ or T = TH) we have 
uw=0, O=O.,; 1.€., F* (9.)=0 or F’ (z,)=0; 
——=0, 1.€., F’” (94) ==0 or FY” (=,)=0. 
Making use of these five boundary conditions, we obtain five equations for the determina- 


tion of Cy, Cy, C3, Qi, and qo, 
F’(l)=c,A, +B, +¢,),= 


' ) 


Cy 
y 3n 
PY (1) Se, Az + ¢,B,-+- C30. = — —: 
2¢ (116) 
F()=¢,A3+ 0B; +¢;D3=% 
PY Gy)=C, Ag+ OB,+¢,D, =0, 


fF" (*,)=¢,A;+¢,8,+¢,D,=0. 


where A;, By, ..., Bs, D3, correspond to the series Fj, F),..., F,, F3; with r =1 and A,, By, ..., 
B;, D, correspond to the series Fj, F),..., Fy, Fy with 7 = 7 yy. 


It should be noted that the solution of system (116) is associated with a number of compu- 
tational difficulties since the coefficients A,, B,,..., Bs, Ds are functions of the initial heating 
and the jet width b,. Therefore, the following procedure is proposed for the solution of system 
(116). 
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For a given initial heating of the jet 7,,, we take a number of values of the jet width b, and 
calculate the coefficients A;, B,, etc. for all values of b, specified. Then, having obtained from 
the fourth and fifth equations of system (116) two expressions for c,; and having equated them, 
we obtain 


D; D. 

Cy =o, Ee . (117) 
i) 4 
As 4 


C,=C; (—£ = —<t =e. (118) 


As a result of the substitution of the expressions for c, and c, into the first equation of system 
(116), we can obtain an expression for the determination of c3(,), 


l l 


Cx) = 7 7oLRPePLn” (119) 
a) C1 Aik + BE + D, 


If the same expressions for c, and c; are substituted into the second equation, we obtain 


C.. =e i (120) 
3(2) 2 ¢f Ak + BoE + D, 
We can now determine c, and co: 
Cy = R53 Q4)3 | 

Cy io = RC,,>,' 

1(2) 3(2)? 
__ ; (121) 

Co1) = ECs); 

C5 (2) = EC). 


The points of intersection of identically labelled curves for some definite value of b, 
which is the same for all curves will then correspond to a unique solution of system (116). Sub- 
situting the values of the constants obtained (cj, co, etc.) into Eq. (114), we obtain an expres- 
sion which can be used to determine the function F(r) and its derivatives. The values of F'(9) 
and F"(g) can be obtained with the help of formulas (90). The longitudinal and transverse veloc- 
ity components, as well as the turbulent friction stress and heat flux are determined from 
Eqs. (38), (49), (50), and (51), respectively (see Table 4). 


3. When n=1, the equation of motion (46) can be written as 


3 1 l 
oP FM (4307 F"()- +2 7 FI) 0F ()=0 (122) 
and solved with the help of series (92) in the same way as the equation with n=2/3. 


Let us substitute series (92) and its derivatives (94) into Eq. (122) to obtain 


@ 


> ag(r-FA)| (FER I) (Fk 2)43 (04 Ky 2)e +c a, ** =0. (123) 


k k=0 


oO 
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TABLE 4 
n=; 6428.62; 121. 16 5 ty=1.757+ 1073; ¢)==—0.72 ; eg=—1.85-10-4 ; ¢3=1 005-1075 
, tl 1 du U “ We ou2 
7 9 t Fs) |F' (9) F°(e) | F's) ly ode | aot | q et | aaouoéy | ou? 
2 
0.0 |} —2.46 1.757-10—3 | —0,.72 [0.0 0.0 0.72 0.0 0.0 0.0105 | 0.00 0.0 0.0 
0.1 | —2.098 0.0302 —0.715 0.0276) 0.133 0.315 | 9.0031 0.0225 | 0.1285 | 0.0045 | 0.0273 | 8.6-1079 
0.2 | —1.736 0.0973 —0.695 {0.0956} 0.243 0.2994, 0.0202 | 0.0768 | 0.2085 | 0.028 0.0668 j1.83-10—3 
0.3 | —1.,376 0.173 —0,.641 |0.2045} 0.348 0.2603] 0.0636 | 0.1644 0.257 0.0858 0.1206 | 0.013 
0.4 | —1.012 0.2618 —0.542 (0,347 0.431 0.1882] 0.142 0.2710 | 0.258 0.1795 | 0.173 0.0493 
0.5 | --0.65 0.3618 —0.388 |0.513 0.481 0.0752} 0.2605 | 0.384 0.233 0,29 0.2203 | 0.134 
0.6 | —0,29 0.4715 —0.1715,0.686 0.48 —0.083 | 0.4155 | 0.478 0.1985 | 0.376 0,251 0.285 
0.7 0.07 0.5909 0.107 |0.855 0.416 | —0.284 | 0.601 0.525 0.17 0.393 0.256 0.513 
0.8 0.435 0.7194 0.44 |0.977 0.271 | —0.522 | 0.784 0.483 0.153 0.291 0.220 0.766 
0.9 | 0.798 | 0.8559 0.806 ;1.032 | 0.032 | —9.7759, 0.931 | 0.315 0.1113 0.1105 | 0.136 | 0.96 
1.0 | 1.16 ' 1,0 | 1,160 1.0 ! —0.288 | 1.022 | 1.0 ! 0.0 0.0 0.0 0.0 1.0 
| | 
Equating to zero the coefficients of like powers of 7, we obtain a system of equations 
. 2 for k=0 aor [(r— I(r —2) $3 — = | =, 
6 rte nal +3)llta\(—3+ 
$3 (r+ )——-] +ea.=0. 
ry (124) 


- ? for k=3 as(r-+3)|(7+2)0-+1) + 

$3 (r-+3)—F] + cay.=0, 
Jt* for k=k ag(r- +A) [FAI r+ 2)+ 
+3(¢r+4)——-|-+ea, ,=0. 


2 


From the first equation of system (124) we find the roots of Eq. (123) to be 


l l 
r=0; re=—; r3=—-— 


Next, from system (124) we find the values of the constants a3 /), a3, etc., for each of the 


roots rj, r2, and, r3. Substituting the values of a3/2, a3, etc. found for r;=0, rp = 17,, and r,=—¥, 


we obtain the three particular solutions of Eq. (122), 
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TABLE 5 
n=13; T,=1.757-10-3; 6,=4.068; 9, =1.327- PEN 17107 1095 040685 $1 = 1 327; Cy =—1. 156 5 Co=—0.09388 5 Cg=— 3.401079 156 ; co=—0.03388 ; cz3=—3.45-1075 
Tr 
Fe) Ir @| reg) | er@ | = | —S} gous || 
Se ee 9 ug 49 Agu 5 agouoig | QotG 
0.0) —2.741 | 1.757-10—3 | —1,158 {0.0 0.0 1.158 0.0 0.0 2.04-10—3 0.0 0.0 0.0 
0.1; —2.334 0.01896 | —1.15 {0.0409} 0.167 0.295 |7.76-10—4]5,82-10—3 0.0399 {1.78-1073} 0.0199 |3.17-10—5 
0,2) —1,.927 0.05452 | —1.117 |0.1325} 0.281 0.268 |7.22-10—3) 0.0299 0.1072 0.0164 0.0603 19.56-10—4 
0.3 —1.521 0.10847 | —1.038 0.268 0.38] 0.222 | 0.0291 | 0.0829 0,187 0.0633 0.1186 | 0.0078 
0.4 —1.111 0.18075 | —-0. cosa. 44 0,456 0.1399} 0.0794 ; 0.1702 0,252 0.1605 0.189 | 0.0349 
0.5 —0.706 0.2714 —(, sist 632 0.483 0.0207; 0.168 | 0.283 0.311 0.295 0.2565 | 0.104 
0.6, —0.301 0.3804 —(). sre 829 0.463 | —0.1528] 0.315 | 0.416 0.367 0.4565 0.318 | 0.261 
0.7 0.109 0.5978 —(). anne. 9975} 0.35 | —0.338 | 0.507 | 0.513 0.4 0.518 0.34 0.504 
0.8, 0.519 0.6535 0.4234/1.112 0.175 | —0.5764| 0.726 | 0.538 0.414 0.4435 0.3135 | 0.811 
0.9 0.919: 0.8176 ().8778)1.122 | —0.1072} —0.7939] 0.918 | 0.391 0.329 0.187 (). 203 1,022 
1.0 1.327 1.0 1.327 0 —().472 | --1.001 1.0 0.0 (),0 0.0 0.0 1.0 
TABLE 6 
n=1; 1,=0,01049 ; 6,=3.87; ¢)=1.314; cy =—1.065; cg=—0.2315 5 cz3=—0.001185 
, ; n u 1 du v ‘ 5 ; ou2 
uF ? t F (9) F’(gy | Fe(g) | F'(¢) uy (| a dp au, ; Ooty aOotigls oot 
0.0) —2,556 | 0.01049 | —1.0972) 0.0 0.0 1.0972) 0.0 0,0 0.0115 | 0.0 0.0 0.0 
0.1] —2.169 | 0.03706 | —1.0896} 0.0483 0.2085} 0.3396] 0.0018 | 0.012 0.0651 | 0.0039 | 0.0288 | 8.65-10—5 
0.2; —1.782 | 0.07964 | —1.0524] 0,1522 0.32 0.26741 0.0121 0.0453 | 0.136 0.0258 | 0.0741 | 1.84-10-3 
0.3} —1,396 | 0.13838 | —0.9653} 0,293 (), 409 0.2033) 0,0405 | 0.107 0,2] 0.0827 | 0.1333 0.0119 
0.4; —1.009 | 0.21317 | —0.8183} 0.466 0.47 0.1113} 0.0992 | 0.2 0.274 0,187 0.2045 0.0462 
0.5} —0,623 | 0.304¢4 | —0.6029] 0.652 0,489 | —0.0131] 0.198 | 0.315 | 0.324 | 0.3265 | 0.2645 0.129 
0.6] —0.236 | 0.41101 | —0.3109] 0.837 0.452 | —0.1791} 0.344 | 0.435 | 0.361 | 0.482 | 0.314 0.280 
0.7), 0.15 | 0.5341 0.0483] 0,992 0.374 | —0.3783) 0.529 0.52 0,384 0.5075 | 0.33 0.5245 
0.8) 0.537 | 0.6733 0.4531) 1.093 0.16 | —0.5883} 0.736 0.523 0.385 0.407 0.295 0.814 
0.9, 0.924 | 0.8286 0.8933) 1.105 | —0,1096) —0.8099) 0.915 | 0.375 | 0,296 | 0.17 0.1915 1.01 
1, 1.314 / 1.1 1.314] 1.0 —0,463 | —0.989 | 1.0 0.0 0.0 0.0 0.0 1,00 


to 
on 
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Fy (7) =1—0.334es8? + 0.012773 — 0.00014 1c729? 4+ 0.000000658c!<" — ..., (125) 
Fy (*)= 7? —0.133c2" + 0,00318c¢7s7.2 — 0.0000257c3=* +0 000000084 Ici — ..., (126) 
F(s) =t-” — 1.333¢7 +-0.089c7<5/? — 0.0014 1c3<4 -- 0.00C00856c!t? — (127) 
The general solution of Eq. (122) is thus 
F()=OF 1 ()boF s(t) + 0,F 3 (2). (128) 


The arbitrary constants c,, cy, C3, aS well as the inner and outer boundaries of the jet, 
gy, and g, are found from the boundary conditions (115). Subsequent calculations involve the 
same considerations and equations as in the case with n=2/3. 


The convergence of the series obtained was checked with the help of D'Alembert's princi- 
ple. The results of the calculations are shown in Tables 5 and 6. The variation of the width b, 
and the inner boundary g, of the boundary layer of a gas jet as a function of the value of n for the 
working medium for two values of the heating ty is shown in Figs. 3 and 4. 


The profiles of the enthalpy 7 and the velocity u/uy in the boundary layer of the jet are 
shown in Fig. 5 for various values of 7, and n. 


4, When n=0, the equation of motion (46) reduces to the equation of motion of an incom- 
pressible fluid 


F(a) + F (2)=0, 
whose solution is well known [1] and can be rewritten as 


ae 
F (¢)=c,e-? + ce * cos 


V3 
2 
An analysis of Eq. (129) and its derivatives shows that this system can have an infinite 
number of solutions. Therefore, we naturally have the problem of determining which of the 
solutions of the systems are physically meaningful. The answer to this question can only be ob- 
tained on the basis of an investigation of these equations with boundary conditions taken into ac- 
count. 


as V ; 
otce° sin = o. (129) 


On the basis of these equations and 
boundary conditions for an incompressible 
fluid, we have obtained numerically two solu- 
tions (two roots of the equation corresponding 
to two values of by = @1 — @_) Which satisfied 
the system of equations and the boundary con- 
ditions of the problem. The first solution 
corresponded to by, =3.02, the second to by, = 
6.655. 


The velocity profiles shown in Fig. 6 
were calculated for these two solutions. As 
can be seen from Fig. 6, the velocities ob- 
tained from the second solution are found to 
be negative in one Section of the jet and this 


: a5 10 situation, of course, should not occur. 
Fig. 3. Variation of the width by, of the boun- Because of this, the second and subse- 
dary layer of a gas jet as a function of n quent solutions must be rejected as being 


(points O represent theoretical results). physically meaningless. 
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It is obvious that the same conclusions 
can be drawn for other values of n, in parti- 
cular, for n=2. 


It should be noted that with n=0, the 
equations for the longitudinal and transverse 
velocity components, as well as those for 
the turbulent friction stress and heat trans- 
fer reduce to the corresponding equations for 
an incompressible fluid. 


6. Exact Solution of the Equa- 


tions of Moti 
Fig. 4. Variation of the position of the inner tions of Motion 


jet boundary g, as a function of n (points O In the solutions of the system of equa- 
represent theoretical results). tions (26)-(29) presented above, we have 

made uSe of an approximate expression for 
the transverse velocity component ti/u)~ F'/3. At the same time, Eq. (36) shows that the exact 
expression for u/uy should be written as 


tof _- qo“ (4), (130) 


so that the solution of the system of equations (26)-(29) can be made more accurate [6]. More- 
over, we assume that the mixing length for heat transport 1; is ~atimes as large as that for the 
momentum transport l,, where a =3V/2¢, 


Let us write the equations of motion and energy as 
F4 > 9'z'+02'"=0, (131) 


F492’ +4 0274 4 627 54 OF (23 220, (132) 
2 2 2 </ 2io <’ 


ff 


Af if 


Samer aa 
Zs Hf nn2 


| | Lays 
EELS Ze 74a 
BS” 22 2anee 


10n 


caeee> ae 
| A 


Fig. 5. Profiles of the enthalpy 7 and the velocity u/uy in the 
boundary layer of the jet for several values of Ty and n: 
Ty = 1.757 107°; — - — Ty = 1.049 107-2 
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The evaluation of the last term of Eq. (132) carried out with 
the help of the approximate expression for u/uy has shown 
that for any value of the heating and any n, the value of 

(2 /t')(z')° is less than unity. In the case of a subsonic high- 
temperature jet, the ratio Up?/2iy is much less than unity. 
Because of this, the last term in Eq. (132) can be neglected. 


re 
——) 43.02 Vy 
= arts 7 


pt 

Pit eT YE Ta 
Pitt TA Pe 
Pitta yy 
Pita tt tg 
_. Pt tt | 

? SSE veee naa Substituting into this equation the expression for the trans- 


ft fsks7] | yy verse velocity component (130), we obtain the equation of 


Fig. 6. Velocity profiles in a motion in the form 
jet of incompressible fluid. 


Differentiating the equation of motion (131) and dividing 
it by °, we obtain 


1 9” 30 oon 
|. g/t. 3” a" —(), (133) 
t 2 4% . i 


, 3. , 1 9%” 2 | 
MSR 4(getyag)} tao 034) 


If we subtract (132) from Eq. (131), we obtain the following expression: 
] a 4 2 : 
— — —]$/2’ — — | dz” — — 2’ — =0. 1 
(— >| +(1 > | = 2! = (135) 


We will not give the solution of these equations in their general form, i.e., for arbitrary values 
of n. 


When @w=2, Eqs. (134) and (135) coincide with the equations derived on the basis of the 
modified Prandtl] model and can be written as 


3.8’ 1 9” y 
Mm Ss __ W\ _ ___ ‘_[+—0 
| re t= tase ; (136) 
1 9) os" 
te pt \yHo. 137 
(— 3 + ;)# ( 


An analysis of Eq. (137) shows that it is meaningless when z'=0 becauSe in this case the velocity 
gradient in the jet is everywhere zero. Therefore, we must equate the bracketed terms to zero. 
The equation of motion (136) and the energy equation (137) can now be integrated. 


Transforming Eqs. (136) and (137) with the help of expression (37), we obtain 
mW 3 <i n n 2’ \2 “71, 
zu — i 2 += ( 1+4)\(=) — an =|: +2=0, (138) 


where 
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On the basis of this, the equation of motion (138) can be written as 


tts I} 4 t ? 
2" (2) 3-2" @) 425 (Sa )2' (2) +2 @)=0. (140) 


Transforming to the variable 7, we obtain 
2” (1) — 3tz” (4) + 2t (tf —ab,) 2’ (x) + Bz (4) =0. (141) 


Its characteristic equation is 


hk? — 3th? +21 (t—ab,)k+62=0 (142) 
and after the substitution of z=k—t, it can be put into the form 


23+3 pz+2q=0, 


where 
2g= —2abt+bs, 3p—=—L—2ab,t. 


The roots of the characteristic equation (142) can be found analogously to those of (72), 
while the solution of Eq. (141) is 


Ul a, a3 a» 
2 (1)=— (1) = ce "t c,e"" cos By +c3,e"" sin By, (143) 
0 


In order to determine c,, Cy, and cz, as Well as the inner and outer jet boundaries, g; and 
wy), we make uSe of the same five boundary conditions as in the case of the incompressible fluid. 


On the basis of the boundary conditions, we can write down the following system of equa- 
tions: 
z (1) =c,e +c,e% cos B+ c,e” sin B=—1, 
2’ (1)=c,e%a, + coe777, + C3e7272=0, 
2" (1)= cyesat + Cx€% (A241 — By2) + C3€™ (A272-+ BI) = — 91D. (144) 
2(0)=c,+¢,=0, 
2! (0) = e401 + C405 -+ c48=0. 


From the first, fourth, and fifth equations of this system we find that 


l a; — Ao ] 


a ) C ’ 
H H 3 


where 


H=e — e% (cos p+ A= ** Sin 8), 
\ 


while g, and @ can be found from the second and third equations. 


The calculations carried out with the help of the improved theory indicate a greater width 
of the jet than the value obtained on the basis of the approximate theory. Thus, with a = 0.08, 
n=2, and heating 7}7=1.757 the value for the width of the jet increases by 9%. For n< 2, the 
percentage increase in the width value obtained from the improved theory is less. This can be 


n=2 
t/’ 
n—z2' 

T 
0.0 0.0 
0.1 0.1136 
0.2 0.3183 
0.3 0.645 
O.-+4 1.128 
0.5 1.76 
0.6 2.505 
0.7 3.183 
0.8 3.52 
(1).9 2.7 
1.0 0.0 


demonstrated in the following manner: 


follows: 


Table 7 contains the values of n(7'/r)z' along the width of the jet and it can be seen that 
this quantity decreases with decreasing n, whereas the magnitude of F'/8= F's, or u/ug, in- 
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It can be seen from formula (130) that the improve- 
ment of the solution depends on the magnitude of the second term which can be written as 


sf 


- 


TABLE 7 
Ty==1.757-10-3 
2 
n= n= 
3 
TO Tt’ 
n—~2 i— 2 
14 T 
0.0 0.0 
0.01985 0.01825 
0.06045 0.04465 
0.1186 0.0803 
0.1885 0.1158 
0.256 0.147 
0.318 0.167 
0.339 0.1704 
0.314 0.1467 
0.2032 0.0903 
0.0 0.0 


anh — 2". 


T= 1.049. 10-2 
n=2 n=! 
t’ t’ 

fl——_ 2 n——2 
20 T 
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creases with decreasing n. As a result of this, the importance of the second term in Eq. (130) 


decreases with decreasing n. An increase in a leads to a somewhat greater increase in the 
width of the jet. However, the velocity profile in the boundary layer of the jet is practically 


independent of the theory used. An analysis of Eq. (136) shows that the latter quantity becomes 


of the form 


i.e., the velocity and enthalpy profiles will be the same (solution of the problem according to 


Prandtl's theory). 


Inthecase of an incompressible fluid, Eq. (134) coincides with Tolmin's equation and can 


be written as 


2’’+2=0. 


The energy equation (135) for an incompressible fluid with ~=1 coincides with Prandtl's 


equation obtained on the basis of his "old" theory, 7"/rT'=z"/z' while with a@=2, it coincides 


with Tolmin's equation 7"/7T'=0. 


In the second case, the enthalpy (temperature) profile in the boundary layer of the jet is a 


straight line. 
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Fig. 7. Plasma generator: 1), 6) current and water 
feed pipes, 3) rear electrode, 4) intermediate housing, 
5) cooling coil, B) places of busbar attachment. 


7. Analysis of Experimental Data and Comparison 
with Theory 


Experimental investigations of high-temperature plasma jets were carried out by means 
of a plasma generator (plasmotron) with graphite electrodes 2 and 3 (See Fig. 7). Arc stabiliza- 
tion was achieved by means of cooling-water which was supplied along tube 6, fed through the 
helical channel 5, and extracted by means of pipe 1. Brass pipes 1 and 6 also served as busbars 
for the supply of electric current to the electrodes 2 and 3. The device was put into operation 
by means of graphite rods 2 mm in diameter. The power supply consisted of a d.c. generator 
with a power output of 300 kW. The stability of the plasma-jet parameters was achieved through 
the stability of the electric and geometric parameters of the source. The temperature along 
transverse cross sections and along the length of the plasma jet was determined spectroscopical- 
ly with an ISP-28 spectrograph by the method based onthe measurement ofthe relative intensities of 
the Balmer-series line Hg of atomic hydrogen [17, 22, 25, 27, 34, 37, 38, 39]. A modification of 
this method of temperature determination was also used in which the whole jet was photographed 
in the light of only the Hg line. Light filters were used to isolate this line from the total spec- 
trum. 


A shadow instrument IAB-451 was used for the determination of the outer boundary of the 
zone in Which the plasma jet intermixed with the surrounding medium. Figure 8 shows shadow 
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Fig. 8. Shadow photographs: a) plasma jet, b) "cold" jet. 


photographs of a plasma jet with nozzle-exit temperature of 17.5-10°°K and a "cold" jet with an 
exit temperature of 340°K. A comparison of these photographs shows that the aperture angle 
and, consequently, the width of the plasma jet is greater than those of a jet of incompressible 
fluid. The increase of the jet width with increasing initial heating is also confirmed by direct 
measurements of temperature and velocity head. Thus, a comparison of the dimensionless 
enthalpy profiles obtained in the initial section of the plasma jet with the temperature profiles 
of the "cold" jet shows that the width of the plasma jet is greater than that of the "cold" jet 
(Fig. 9). Analogous results are obtained from a comparison of the dimensionless velocity—head 
profiles in the initial section of a jet heated to a temperature of 4000°K with those of an isother- 
mal jet (Fig. 10). A proof of the correctness of such an assertion is also provided by experi- 
ments carried out in the main Section of the jet. This can be seen from from a comparison of 
the velocity-head profiles of a "cold" and « hot (T =4000°K) jet (Fig. 11). As can be seen from 
Fig. 12, an increase in the initial heating leads to an increase in the width of the jet. 


The increaSe in the jet width with increasing initial heating is in agreement with the re- 
sults of the theoretical calculations reported in this paper. 


An increase in the intermixing Zone with increasing initial heating leads to a decrease in 
the length of the initial section (Figs. 13, 14). It can also be seen from Fig. 14 that even on the 
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Fig. 9. Comparison of the enthalpy profiles in the boundary 
layer of a jet: —— theory (plasma jet py;/p)=245, a=0.143), 
— — — theory (incompressible fluid, a=0.09). 
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Fig. 10. Comparison of velocity-head profiles in the boundary 
layer of a jet: ——-theory (isothermal jet p;,/p)=15, a=0.127), 
— — —theory (incompressible fluid, @=0.1). 
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Fig. 11. Velocity~head profiles in the main section of a jet: 
Data for the nonisothermal jet are from V. Ya. Bezmenov 
and V. S. Borisov's experiments, those for the incompres- 
sible-fluid jet are from the author's experiments. 


axis of the initial section of the jet, the temperature and, consequently, the enthalpy do not re- 
main constant, but fall to some extent. This temperature decrease can apparently be explained 
by noting that a part of the energy contained in the jet is radiated away. The small temperature 
decrease at such high values of temperature can only be explained if it is assumed that the 
blackness of such a plasma jet is very low and is only 0.001-0.002. This value of the blackness 
was determined as the ratio of the radiant energy (measured experimentally with the help of a 
"Tera-50" radiometer) tothe amount of energy radiated by a black body whose temperature is 
equal to that on the axis of the jet. The shadow photographs show that the outer boundary of the 
plasma jet remains rectilinear, as is the case for a jet of incompressible fluid. In the case of 
strongly heated jets and, in particular, plasma jets, the velocity-head, velocity, and enthalpy 
(temperature) profiles remain affine as can be seen from Figs. 9, 10, 11,12. The linearity of 
the jet boundaries and the affineness of the parameters verify the assumption, adopted for the 
calculation, that the mixing length increases linearly along the jet. 
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Fig. 12. Velocity profiles in the main section of a jet: Data 
for the plasma jet are from G. G. Smolin's experiments, those 
for the nonisothermal jet are from V. Ya. Bezmenov and 

V. 5. Borisov's experiments, and those for the incompressible- 
fluid jet are from the author's experiments. 
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Fig. 14. Variation of the enthalpy along the jet axis; 
O plasma jet, pyj/pp = 245, author's experiments; 
—~-— jet of incompressible fluid. 


Figure 15 shows the dimensionless velocity and temperature profiles obtained in experi- 
ments with a jet heated to 4000°K, and it can be seen from the figure that the temperature pro- 
file lies above the velocity profile. This result is in agreement with the assumption used in the 
calculations that the mixing length for heat transport is greater than that for momentum trans- 
port. Figure 9 shows comparisons of the experimental enthalpy profiles in the initial section 
of the plasma jet and the temperature profiles in a "cold" jet with theoretical data. A compari- 
son of these results show that the experimental data are in good agreement with theoretical 
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Fig. 15. Velocity and temperature profiles in the initial sec- 
tion of a nonisothermal jet, 0;;/p)=15, according to V. Ya. Bez- 
menov and V. S. Borisov's experiments. 


calculations and that the enthalpy profiles are not universal but depend on the initial heating. 
Figure 10 gives the experimental data on the dimensionless velocity-head profiles obtained in 
the initial section of both an isothermal jet and a jet heated to 4000°K, together with the theo- 
retical results; it can be seen from the figure that there is satisfactory agreement between 
theory and experiment. 


The results obtained in this paper show that it is possible to calculate the parameters of 
a very-high-temperature jet with an accuracy that is sufficient for practical purposes. 
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THE CALCULATION OF TURBULENT JETS OF REAL GASES 


V.I. Bakulev 


Much attention is being devoted at the present time to problems of propagation of turbu- 
lent jets of incompressible and compressible gases, supersonic jets, as well as to problems of 
the propagation of very-high-temperature jets. 


One of the interesting problems in this area concerns the propagation of turbulent jets of 
real gases. 


By a real gas (See Fig. 1) we mean a substance situated at a pressure above the critical 
pressure, as well as to the right of the upper boundary curve 2, i.e., a one-phase working sub- 
stance. 


For processes taking place at constant pressure, we give an approximate equation of state 
for a real gas, or, more accurately, the equation connecting the enthalpy with specific volume. 


In the following, we consider three problems of the propagation of jets of real gases: 
1) The boundary layer of an infinite plane-parallel jet (initial region of the jet); 
2) A plane-parallel jet emerging from a very narrow slit (main region of a flat jet); 


3) An axially symmetric jet emerging from a very small aperture (main region of an 
axially symmetric jet). 


The calculations are based on the Prandtl theory of turbulence and the assumption that 
the mixing length for vector quantities is 1/n times as large as those of scalar quantities: 
nly =lj, where n=1/Pr7 (Pr7 is the turbulent Prandtl number). In solving these problems, we 
assume that the thicknesses of the dynamic and thermal 
boundary layers are equal and that J,,=cx for both sub- 
merged jets and jets in a parallel stream. 


1. Thermodynamic Processes in Free 
Turbulent Jets of Real Gases 


The discharge of the working fluid through the 
jet nozzle into a chamber and the subsequent intermix- 
ing of it with the surrounding medium occurs at appro- 
ximately constant static pressure p=const. if the jet 
does not interact with the chamber wall or other jets. 


Fig. 1. pV diagram of a real gas. If the working pressure in the chamber is above 
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6.67°1074i, J/ke the critical pressure for the given sub- 
stance, then the thermodynamic process 
of nonisothermal] jet mixing will be repre- 

—— sented on a pV diagram by a straight line 
which will lie above the liquid-vapor phase 
of the substance (See Fig. 1). 


6000 


In this process, the working sub- 
stances will have no heat of vaporization, 
i.e., the parameters in the boundary layer 
of the jet will vary continuously without 
the formation of droplets and their gradual 
evaporation. In other words, at super- 
critical pressures the substance will have 


4000 


2000 only a single phase [1-4]. 
An analytical expression describing 
the state of a real gas is necessary for 
jet calculations. Qualitatively, a real gas 
can be described by van der Waals' equa- 
0 20 40 60 tion [5, 6, 7], but there is no single quanti- 
9.81 -10°v, m°/N tative equation of state valid for all sub- 
stances. The literature [8, 9] contains 
Fig. 2. Variation of the enthalpy of carbon di- equations of state for real gases, but these 
oxide with specific volume at p=981 -104 N/m’: are different for different substances. In 
tabulated data,-—--: - approximating addition, these equations are very com- 
curve. plex. 


We will attempt to select an approximating function to describe the variation of the enthal- 
py of a real gas with specific volume i=f(V). 


The function giving the variation of enthalpy at pressures above the critical pressure [10- 
15] for processes that take place at p=const. is of the form of a cubic parabola for all substances 
(see Figs. 2, 3). Let us take the approximating function to be also a cubic parabola, namely, 


| C ,D 
j=AV+B424—, (1.1) 


where the coefficients A, B, C, and D for any given substance are functions only of the pressure. 


This equation describes with sufficient accuracy the behavior of the enthalpy of various 
real gases (see Figs. 2, 3) at temperatures up to 2500-3500°K. 


If the accuracy of Eq. (1.1) is insufficient, we can replace this equation by a sum of func- 
tions involving any noncomplex powers of the specific volume, 


jxuan ; 


jank 


where k and n are any real positive numbers forming a Sequence j. 


With the help of Eq. (1.2), we can approximate with high accuracy the real behavior of the 
enthalpy of a real gas in any temperature range. However, the values of the coefficients A; in 
this equation cannot be easily determined since for this we require the solution of a system of 
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Fig. 3. Variation of the enthalpy of nitrogen with specific 
volume at p=393 -10* N/m?: ——tabulated data, — - — ap- 
proximating curve. 


(n + k) linear algebraic equations. For further calculations of jets of real gases we can use 

any of the equations given above. The simple Eq. (1.1) is fully adequate for the study of jet 

flows at pressures above the critical pressure and at moderatetemperatures, so that in the follow- 
ing we will only make use of this equation. 


In order to calculate the temperature field in the boundary layer of a jet of real gas, it is 
also necessary to know the variation of the specific volume with temperature V=f(T) for the 
given working substance at p=const. (Figs. 4, 5). It is not necessary for us to have an analyti- 
cal expression for this dependence, since the temperature field in the boundary layer of a jet 
can be obtained graphically. 


Equation (1.1) should be fitted as closely as possible to the graphical dependence i = f(V). 
We will therefore use the least-squares method [16] for the determination of the coefficients 
A, B, C, and D. 


TH The sum of squares of the deviations 
ou EERE EE A between the approximating curve and the 
TAT actual points (see Fig. 2) should be reduced 
pan to a minimum: 
AEE 


n 


ls (AV y+ B+ tae)) min, 


j=1 


Pt ttt tty 
Se 
500 ptt ttt tT tT tT tt ttt . , 
PPT Prd TPP re where n is the number of points selected 
0 5 10 15 20 from the actual curve. 
9.81°10°v, m°/N 
Fig. 4. Variation of the specific volume of Differentiating this equation with 
carbon dioxide with temperature at p=981 - respect to A, B, C, and D, we obtain the 


10* N/m. system of equations 
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Substituting these coefficients into system (1.3), we obtain four linear algebraic equations 


for the determination of A, B, C, and D. 


a,,=aA+aiB--C 4t-a" 1D, 


awi=aAtB+La iC+a oD, 


(1.4) 


a1; Ata_jB+a.C-+a™ sD, 


* * * 
ES aiA 4- a_oB 


+ a_3C + a_4D. 


It should be emphasized that the values of the coefficients A, B, C, and D obtained inthis 
way will be valid only for a given working substance and a process which takes place at constant 


pressure. 
ficients. 


pp tt EE ET 
pi | ttt Per 
at | | 
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50, 
9.81-10°v, m°/N 
Fig. 5. Variation of the specific volume of ni- 


trogen with temperature at p=393 10! N/m’. 


A change in the pressure or another working substance give other values of the coef- 


2. Initial Region of a Plane- 


Parallel Turbulent Jet of a 
Real Gas 


Formulation of the Problem. 
Let us consider the case of mixing of two 
plane-parallel infinite streams of real 
gases (See Fig. 6). The main stream flows 
with a velocity ug (density py), temperature 
To), While the other stream is flowing in 
the same direction with a velocity uy, 
(density, py, temperature Ty). We will 
place the coordinate origin at the point at 
which the intermixing starts. 
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uy ly ll, Starting at the point of joining of the two streams, a 
free turbulent boundary layer is formed which gradually widens 
! in the direction of motion. The inner boundary of the boun- 
- << ; “ary layer will be denoted by line 1, the outer boundary by 
ine 2. 


In the boundary layer, the velocity, density, and temper- 
ature change smoothly from their values on the inner boundary 
(Ug, Po, To) to the values on the outer boundary (Upy, Py, TH)- 


In order to determine the profiles of the velocity, den- 
sity, and other parameters in the boundary layer, we will 
Fig. 6. Boundary layer of a make use of the equations of motion, continuity, energy, and 

jet. the equation of state [17-26]. 


The averaging of the equations of motion, continuity, and energy is carried out according 
to Reynolds's method [24, 25], while the Prandtl theory [17, 24, 26] will be used for the descrip- 
tion of turbulent friction. 


Equations of Motion and Continuity. The equation for a plane-parallel nonstationary mo- 
tion of a viscous gas in the direction of the axis of abscissas is of the form [17] 


Ou | Ou Ou dp, Wee OF ey 
Ol GU = —_——-— es _ 
ar} ox aa oy Ox | + (2.1) 


Ox Oy 
where u and v are the instantaneous velocity components, p and p are the instantaneous values 
of the density and pressure, Tx x and Txy are the components of the normal and tangential stress. 


Let us also write down the equation of continuity for the instantaneous values of the veloc- 
ity and density of a viscous gas in the nonstationary plane-parallel flow, 


V5 ; ) 
0: 4 0 (ait) -+- ( 


Ot OX Oy 


(sv) =0. (2.2) 


We will replace the instantaneous values of all quantities by the sums of their averaged 
values and their fluctuations. In doing this, we will make use of van Dreist's hypothesis that in 
addition to velocity, density, and pressure fluctuations, there are fluctuations of the current 
densities pu and pv as single entities [17, 27, 28]. 


In this case, the equations of continuity for the averaged quasi-stationary motion (8p/at = 0) 
will be 


—_— 2] — 
—°_ (pi) + (50) =0. (2.3) 
OX Oy 


Making use of the equation of continuity (2.2), we can write Eq. (2.1) as follows: 


0 
at 


° (ou-u) + u (vu) — SP Ox 
Ox ay ax x 


(9-u) 


Oy 


In this equation we now replace all instantaneous values of various quantities by their 
averaged values and obtain the equation of motion for a plane-parallel quasi-stationary motion, 


Op Ory Oxy 


(ou’) + 


“ (-u)+ 


O 
Ot Ot 


2) — / ' ) — ) , 
° (san) + & [uy] -- (0-0) + © [oyu] = 
x OX Oy Oy 
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Assuming that in turbulent jet flow the normal and tangential viscous stresses are small 
by comparison with the turbulent friction stress, that the pressure is constant, that the deriva- 
tive of the correlation of fluctuation quantities with respect to x is small by comparison with the 
derivative with respect to y, and that the motion expressed in terms of the average quantities is 
7 0 
Ot (4) =: ot 
aged two-dimensional turbulent stream 


O 


a quasi-stationary one (Wu) =0] , we obtain the equation of motion of an aver- 


QO =>. 8 O —— 
ait+tt) -- —— (ov-u)== — — [(nv)''). 
ott) 2 (6 tt) = = (GRY 


Using the continuity equation (2.3), we can rewrite this equation as 


=~ [aye], (2.4) 


where (pv)'u'=77 is the turbulent friction stress. 


The expression for the turbulent friction stress can be put into the following form [17]: 
p= (oto )(otv')u' eyo, 


since here pvv'=0 and the terms p'v'u' and vp'u' are negligibly small. 


In the case of an incompressible-fluid jet, we will assume as in the Prandtl theory that 


According to dimension theory, we find that the so-called mixing length 1, has the dimensions 
of length. In view of the absence of walls [ly(y) =const.], we take 1, =cx for both a submerged 
jet and a jet in a stream moving in the same direction (17, 18, 28]. 


The formula for the turbulent friction stress is 


Ou 
Oy 


s Ou 
t= + c?x*9 — 
Oy 


Here, the plus sign is taken when 0u/dy > 0. 


Substituting into Eq. (2.4) the expression for turbulent friction, we obtain the equation of 
motion for a plane-parallel quasi-stationary jet flow, 


=e 


— Qu — ou oo 0 [—-/ du \2 
f Maa —— — HO XX —|!. . 
ary ve ax dy  ( )] (2.5) 


The energy equation for nonstationary plane-parallel flow is of the form [17, 26] 


7 (pp eee) _ op 8 3 09x , Og 
(i +3 j= +2 (Wet ay) + (rye Oy) + (E+), (2.6) 


where Txy =Tyx are the components of tangential stress, T,,, Tyy are the components of normal 
stress, i=c)T is the enthalpy of the gas,.and qy, and qy are terms characterizing the amount of 
heat being transferred per unit time through a unit area. 
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Let us transform the cnergy cquation with thc help of the continuity cquation (2.2), 


2 7] rn ar 
StS) (£5 2)+ 
ox \ 2 2 


() og 
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65a Buty Lees eyp (te) 

vy OX oy 

Let us replace the instantaneous values of quantities appcaring in thc cncrgy equation by the 
sums of their average values and fluctuations. We will also assume that thcre exist fluctuations 


of enthalpy as a single quantity [28]. 


Assuming that the cffects due to molecular viscosity and molccular heat conduction are 
small by comparison with turbulent-mixing effccts, that the molecular friction strcsscs arc small 
by comparison with turbulent friction stress, that prcssure is constant, that thctransverscvclocities 
are considerably smaller than the longitudinal vclocities and that the motion is quasi-stationary 
[17], we obtain the following approximate energy equation for an avcraged two-dimensional tur- 
bulent jet: 


— — (sui) +— (i) + — [(cvy i] + (sv) a’ i (). 
( 
In this equation we have ncglected derivatives of fluetuations with rcspcct to x as being 
small by comparison with the derivatives with respcct to y, we have neglccted the terms 
av: v-av/ay and (v"/2) -a(ev)/ay (since vis small by comparison with &), and wc have also sub- 
tracted the equation of motion and the continuity equation (2.3) from it. 


Taking into account the smallness of the quantity v, we can write 
(ov) t" A 


Then, making uSe of this relation and the transformations made above for the equation of motion, 
we will obtain the energy equation 


— Jo 0 Ou 
Ot 90 Ties — [pol pou! , (2.7) 
dy oy 


— Oa 
VU 
| OX 


The heat flux W7 in a two-dimensional turbulcnt stream is governed hy the correlation 
between the transverse component of the velocity fluctuations and the enthalpy fluctuations 


W, = oUt’. 


Making use of the Prandtl turbulence model, we can obtain the following cxpression for 
the heat flux: 


Ou 


“oy. 


Wo + oh as (2.8) 


— iru 


where J; is the so-called mixing length for heat transfer. 
The plus sign corresponds to 0u/dy > 0. 


Lect us express the turbulent Prandtl] numbcr in terms of the mixing lengths: Pry =l,/lj 
(28, 29]. 


Let us express lj in terms of ly and Pry and substitute it into expression (2.8) aftcr sct- 
ting n=1/Prz7. 


44 V. 1. BAKULEV 


We obtain an expression for the heat flux in turbulent flow 


Oi 


—— ee 


Uy 


Ou 


——— 
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W - 3 ali 
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Substituting into Eq. (2.7) the expression for the heat flux, turbulent friction, and ly, we 
obtain the energy equation for turbulent plane-parallel quasistationary jet flow, 


— J ~,— Jd - >» O |> du dt » oof Ou \3 
Ou — j-l OV —— (== +NC-X- — a —_—- CoXW( (=) . 2.9 
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Transformation of the Equations for the Turbulent Boundary Layer 
of a Jet of Real Gas. In order to solve the equations of motion and energy we will intro- 
duce the new variables 


XxX; c=yfar, 
3 y) o e 
where a= v2c2 tS an experimental constant. 


In this coordinate system, the velocity, current density, gas density, specific volume, and 
enthalpy will be functions only of 


— —_— — — 


t= Z (2), far), £=9); —=8@): —=+(2). 


uy) — tty, ayy 40) Vo by 


Here, Vp) and i, are the values of the specific volume and enthalpy on the inner boundary 
of the jet and Vy, and iyy are the values of these quantities on the outer boundary; F'(g) =dF(g)/dq@ 
is the derivative of a dimensionless function of one independent variable g. 


In this case, the current function is 


v= | sudy=ttyaxF (2). 


a 


From this, it is easy to obtain the transverse component of the current density 


— Oy 


C= o_O 
OX 


oy _dF (eg) oO 
= — figllod | (¢)-+. Te nal 
or, 1n final form, 


OV = Alglly (¢F' — F). 


Substituting into the equation of motion the relations given above, and taking into account 
the relations 


omer eee eee OC 
—— —S Ol 


=m+(l—m)Z (9); ma, 


9] 


we obtain 


~F=+(1 —m) (0 Z' +292"). (2.10) 
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Let us write down the relationship between the main velocity component u and the current 
density pu 


ation -Lo'u’, 


Following the Prandtl model of turbulence, we assume that the correlation p'u' can be 
expressed in terms of the derivatives of these quantities with respect to y, while the coefficient 
of proportionality 7, for the density can be taken equal to 1;. We can then write 


— 77 Oo Ou 
ou ou + Lf, — . 
7) dy 


Substituting into this equation the relations obtained earlier, we finally find 


F’=9 [m-(1—m) Z| +a (lm) 9'2", (2.11) 


Differentiating the equation of motion (2.10) with respect to ~ and replacing F' by its value 
from expressions (2.11), we obtain the equation of motion in the following form: 


, rr ae noof Et Uw ye re 
VA ef (an po ft —_..-— = 
75 "a tla y ' 9 > [2+ 247" 0. (2.12) 
Let us now replace the density function @ in Eq. (2.12) by the specific volume function 
9(3 =1/9) and transform to the new independent variable n =(~ — @)/(@, — @), Where @; — q) =bx 
is the dimensionless width of the boundary layer. We will then obtain the equation of motion 


for a boundary layer of a plane-parallel turbulent jet of a real gas 


_ zu 3 0" Wy) — lq oo 1  /v',\ TP! (4 37 (4 mn 3 
L(Z)=2"(9) — > 2") — a> 85 +|= (<-)\2 (1) + O22 (1) + B= 0. (2.18) 


In solving the equation of motion (2.13), we make use of the following five boundary condi- 
tions. 


On the inner boundary of the boundary layer we have [n=1, (y= q)]: 


1 u=Ny, Z(1)=1; 


2. 7 —0, Z’(1)=0; 


3. pu=0, 2" (1)= —eb2— 


| — m 


(2.13a) 


On the outer boundary of the boundary layer we have (n=0, (y= q)): 


4, u=u,, Z(0)=0; 


~ Ou 


o.— =0; 2’ (0)=0. 
dy 


Let us transform the energy equation (2.9) using the same relations as were used in the 
equation of motion. As the result of such a transformation we will obtain the energy equation 
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in terms of the independent variable g, 


— Fl = (lam) [210 $92" 2 92Z73"| (1 — my a d(Z')", (2.14) 
< by 
In the case of appreciably subsonic flow, the last term in this equation can be neglected 
since u’, << 2i) and, consequently, this term will be small by comparison with the others. 


Multiplying Eq. (2.10) by r' and subtracting it from Eq. (2.14), we obtain 
RY L AR SIZ T 4-0Z's" =0, (2.15) 


where Ky =1-1/n, ky =1~2/n. 


If we now divide Eq. (2.15) by $Z't', we will obtain the energy equation for a boundary 
layer of a plane-parallel turbulent jet in the form 


9’ ZL raid 
Ip StS 0 (2.16) 


We will make use of the following boundary conditions for the solution of Eq. (2.16). 


On the inner boundary of the boundary layer we 
have o—9, ,i=iy, b=, W—=Uy and, consequently, 
t(a)=1, ¥()=1, 2)=1; 


On the outer boundary of the boundary layer we (2.1Ga) 
have o=o,, i=i,, >=, “=U,, and, consequently, 
(2) = Ht ty, Z (2)-=0. 


= (22) = —Tt H? 
by 


HO? 
a0 


Solution of the Equations of Motion and Energy. Let us rewrite the equa- 
tion of motion (2.13) as follows: 


— IM 


3 
5 — 7" 
L (Z) 7 


where 6'/@ and ho" /6 — (a'/6)° =o(n) are unknown functions of 7, to be determined during the 
solution of the energy equation. 


We will seek the general solution of this equation with the help of Galerkin's method [380, 


31], 
jak 
Z,4=)> Li yy (j=1, 2,3 see. k) 
4-1 


For this purpose, we will consider the system of functions 
o == 243 — 377; oy =(I —4)? 47; Og3=(1 —7)"7P; 20... 2,=(1 — 7)°7/. (2.18) 


For any Qj» this system will automatically satisfy three of the boundary conditions (2.13a): 
Z'(1)=0, Z(0)=0, Z'(0) =0. 


If we take Oly =—], then system (2.18) will also satisfy the fourth boundary condition Z(1) =1. 
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We will seek the solution of the equation of motion in the second approximation, namely, 


Z,= — (24) — 37") +a, (1 "yn". 


The unknowns in the second approximation are a», by, and 9. 


Let us write down three equations for the determination of these unknowns 


1 1 
| L(Z2)ed4=0, | L (22) p.d4=0, 
Q 0 


Z3(1)= — 9,67. 


e 
1 — om 


Integrating the first two equations, we obtain 


n 


6 1 A, 4 a 6,43 + B,— 0.37269 — 0.5 F 


2 
ay, = 


m 


fl 
2.44 Ay-+ a 6,Ay + By + 0.016768 


n 


Q, = 


Ay -+ a 


In these equations we have 


—O.4-+As+ a > 6 A. + By +0.01676° + 0.0333 ; 


ft 4 3 
_ * 


i 


— b,Ay + By — 0.015903 


2 


A, = 36/, — 72/, + 27/,, 

A, ==36/, — 90, + 60J,-- Vp, 
A,=12/, — 30/, + 18/3, 
A,=8/,—24/, + 22/, — 6/;, 
A,=18/, —+5/,+ 36/;, —9/, 

A ,= 184, -— 54/,-+ 574, — 247, + 3/0, 
A,=6J,— 18J,-+ 18, — BJs, 
Ay=4J, — 14/, + 184; — 107, + 2/5. 
By =12), — 30/,+- 18), 

B,=8/), — 24), 22), — 6fs, 
B,=6/,;—18/;+ 18/,—6/3, 
B,=4j,—14/, + 18/; —10/,+2/,, 


where 


1 1 
6’ ; > 

I=) Std flo Cdatdn (b= 2, 3, 4...) 
0 0 


(2.19) 


(2.20) 


(2.21) 


The integrals J), and jy, can be evaluated analytically or they may be obtained graphically 


after the energy equation has been solved. 


Solving system (2.20) graphically, we find the coefficient Ol and the dimensionless width 


of the boundary layer b, - 


In order to find the third unknown qg,, we substitute the value of Z"» at n=1 into the fifth 
of the boundary conditions (2.13a). As a result of this, we obtain 


On I 
e1= (6 — 24,) — (Im). 
b 


(2.22) 
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In order to determine the convergence of the method, we uSe analogous procedures to find 
the solution of Kq. (2.17) in the third approximation, 


ZL — (298 — 307) + a, (1 = 4)?7?- ay 1 = 4)? 
Let us write down four equations for the determination of the unknowns Op; Q3; by , and gy. 


1 1 

1. (L(Zs)ady=0, 2. \L(Zs)9.d4=03 
0 0 
, ~ re 

3, \ L(Z3)¢,dy,=:0; 4. Z3(1)= — 9,0" 


0 


l—m_ 
Substituting into these equations the values of Zs, Ze, Z3', 23", O13 Do 03; we obtain 
after some transformations the following system of equations: 


— (6 4 Aya 6, Ag+ By ~ 0.3726) 0.56} )+ 


* 
nM 


413, (2.4- A.ta ~. b A, + By-} 0.01678 )+ 
7 (22+ Ay-Fa 2 b, Ay +B; +0,009926) )=0; 


2. —(—0.4+ As-+a— b,A,+ By-|-0.01670, + 


-+.0.0333 —* 43 )+4, (Ac +a b Ay} B , — 0.001596, }+ 


l1— m 


. 2.23 
1a, (0.02857 + Ay +a b,Ay+ Bo 0.00086, |= 0 12.28) 
3. - (9.2 — Ayy— a b, Aig— B; ~ 0.009920} — 
— 0.01667 "9 )+4, (0.02858 — Ayy—a 26, Ay ~~ Bob 
-— m 
+.0.00088° +4, [— Aisa 2b Ayy— By +.0.0004330)) -=0; 
4. 9, = (6—2a,—24;)— (1 - m). 
b° 
In these equations, we have 
A,=60/,— 162/, -+ 126/, — 27Js: 
Aip= 10, —31/g + 304, — 9/4; 
Ay = 30/, — 96/, -- ] ] lJ, — 54) ,-+ 9/3; 
Ayn = 5, — 18J74- 24), — 14/43) 
Ais= 18/5 — 45; 4-36, — 9: 
Ay, = 18J; -- 54) ,-+ 57, — 24d, + 343 
Ay, = 30/,—- 96J, 4- 111, — 54/5 4-9; 
A,,=6/, — 18/,+ 18/, — 6/,; (2.24) 


Aj =4/g— 14/, + 18/, — 10/, + 2/,; 
Aig=OJ/ >" — 18/,4-24/,— 14), 4-3/,; 
B,=10/,—31/,+30/;— 9/,; 
Be=Sjxy— 18), +24 j¢— 14 ji, 4+3)/ 4; 
B,=6),—18/,+ 18/;—6/,; 
By=4tjy—147,4 18),— 1Oj,+ 2/5; 
By =9fy—18/,+ 24),—147,+3);. 
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The coefficients A; + A3, B; + B, not listed here, as well as the expressions for J}, and jj, 
have been given above when we considered the second approximation (2.21). 


We eliminate Ql and Qs from the first three equations of system (2.23). As a result of 
this, we obtain one equation for the variable b,. The solution of this equation for by is found 
graphically. Having determined b, , we can use system (2.23) for the calculation of @, 023» and 
71° 

The analysis of the approximate solutions given above can be obtained by comparing them 
with the exact Prandtl—Tolmin solution. For this purpose, we will first calculate the velocity 
profile in the boundary layer of an incompressible-gas jet in the second approximation after 
determining the values of a, by, and g, (@,=—1.532, b, =3.037, g,=0.983), 


Z =: — (243 — 377) — 1.582 (1 — 4)? 4. 
Let us also calculate the velocity profile inthe third approximation (a, =—1.479, Ol =~ 0.051, 
bx =3.03, gm, =0.988), 
Zy== — (243 — 3747) — 1.479 C — 4) 7 — 0.051 (1 — 4)? 7°. 
Finally, we will give the results of the Prandtl—Tolmin solution [17, 18] (bx =3.021, 
y, = 9.981). 


The results of these calculations are given in Table 1 which shows that the agreement of 
the second and third approximations with the Prandtl—Tolmin solution is satisfactory. More- 
over, it should be noted that the third approximation does not lead to any significant improve- 
ment by comparison with the second. 


To check the convergence of the solutions obtained, let us also calculate the initial region 
of a jet of carbon dioxide (CO,) at p=981 -10! N/m, T) =248°K, up =0. 


The velocity profile in the second approximation (@=1.725, by =1.86, g,=0.737) is 
Zo= —(2n8—3n?) + 1.725 (1—n) 27? 


TABLE 1 
1 0 0.1 0.2 0.3 04 | 0.8 0.6 
Zo=— 0 0.0156 | 0.0648 | 0.1484] 0.2637 | 0.4042 | 0.56 
ug 
Z,= — 0 0.016 | 0.0661 | 0.1501 | 0.2656} 0.406 | 0.5612 
ig \ 
ii | 
_— 0 0.0162! 0,065 | 0.147 | 0.261 | 0.405 0.56 
Th | | : | 
| | . 
Continuation 
y 0.7 0.8 0.9 l 
Zo=—- | 0.7165] 0.856 | 0.961 Second approximation 
ig 
Zy=— | 0.7173 | 0.95°6 | 0.9611 1 Third approximation 
UQ 


Fi=-— 0.716 0.855 0.961 | ] PrandtlL— Tolmin solution 
| 
| 
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TABLE 2 
q 0 Q.1 0,2 0.3 O.+4 0.5 0.6 
Zo=— 0 0.012 | 1.1481 | 0.292 | 0.4512 | 0.6078 | 0.7471 
ig 
Z3=— 0 0.0417 | 0.1473 | 0,291 | 0.45 | 0.6064] 0.7466 
iy 
Continuation 
r 0.7 0.8 0.9 1.0 
Z=— 0.8599 | 0,9409 | 0.9804 l Second approximation 
lg 
Z4=-— 0.8595 | 0.9407 | 0.9803 l Third approximation 


ilo 


and the velocity profile in the third approximation (a> =1.685, 05 = 0.0452, bx =1.85, g, =0.742) is 
Z3= — (2n3—3n’) + 1.685 (1 —n)) 2n? + 0.0452 (1 —n?) 3. 


The calculations are given in Table 2. 


As shown by the calculations, the convergence of the solutions is satisfactory. The sec- 
ond approximation is very little different from the third, so that it is recommended that the 


second approximation be used in the future. 


Let us now proceed to the solution of the energy equation (2.16), 
9’ Z" a“ 
“| = 0 
ky 5 +h, 7 4 7 


with boundary conditions (2.1Ga). 


As the result of the integration of this equation, we obtain 


Dic! a= ¢, (Z’)- (2.25) 


Replacing ¥ by 69 (8=1/6), » by 7 [7 =(~ — gd) /(G1—@)] in Eq. (2.25), and introducing the 
equation of state i=/(V) after putting it into the dimensionless form 7+ = f(9), and differentiating 


the latter: 


<= f" (8) = 


ar, 


Wwe obtain the energy equation (2.25) in the form 


i) 0 =) (Z’)- #4. (2.26) 
(j*a 


Integrating Eq. (2.26), we obtain 


\aE d9==C) \ (2) + C2. (2.27) 


The left-hand side of this equation can be easily integrated if the equation of state is taken as a 
polynomial in V {Eq. (1.2)]. Recommendations governing the choice of the approximating func- 
tion are given in Section 1. 


V. I. BAKULEV ol 


For jets of real gases at moderate temperatures we will use Eq. (1.1) which after dif- 
ferentiation becomes 


9 
ie’ =(AV, 2-8 _\v (2.28) 
V yf? 503 


Substituting expression (2.28) into Eq. (2.27) and integrating the left-hand side of this 
equation, we obtain 


AV 


0 ————  — LE sO 
l—ky Vo(l +4) 68 | Ve(2+h,) & oe / 
=; \ (Z’)-*edy, 4-9. (2.29) 


gPa . 


There is a remaining integral on the right-hand side of (2.29) and this can be evaluated during 
the simultaneous solution of the equations of motion and energy as follows. Let us take a pro- 
file for an incompressible fluid or a profile obtained from the solution with n=2 (k,=0). Using 
Eq. (2.29), we determine the 9 profile, while the improved profile of the velocity Z can be found 
from the equation of motion. Substituting the new Z profile into Eq. (2.29), we repeat the calcu- 
lations. 


In this way, we use the method of successive approximations to determine the Z and @ 
profiles. The constants c,;* and c)* are determined from the boundary conditions (2.164). 


Making use of the solutions of the energy equation (2.29), we can find the expressions for 
the integrands of the integrals J), and jy. 


The function 9'/6, appearing in the J), and jj integrals is of the form 


ee ee (2.30) 
Cc 2D 1 
AV — —— —- = — 
VY, yz: 
while the function 9"/9 appearing in the integral jj, can be written as 
we ct" (k, +1) yeh (Z')72* Ret; Qs (2’)— Fat) Vil ct g(2Aitl)( 7/y—2ke Y’ (0) (2 31) 
oS y2 (0) Y (8) ¥3 (6) | 
where 
Y())=AV p-C 4122 1 
Vp 4 Ve 02 
Cc 4D 
‘(0)= AV, + —— — + —. 
VO= Mot ete a 


In order for us to calculate the integrals J, and jj, we have to construct graphs of the 
functions 9'/@ and o(n) = 46"/e —(9'/6)’. In general these integrals can only be evaluated 
graphically [32]. 


If we assume that Pry =0.5 (n=2, k,=0), then Eq. (2.29) and the functions (2.30), (2.31) 
can be simplified. 
The solution of the energy equation (2.16) will be 


l C 
l— ak, Vo(l + &) 


qs 


2D l 


| 
G Ve(24+ &) © 


AVy 
=c17,+ 79. (2.32) 
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When Eq. (2.32) is used, the @ profile is found directly from this equation. 
Expressions (2.30) and (2.31) with Pry =0.5 become 


* Ry 
l 9 1? “ 
AV’, 4 — — — — ; - 
" ( Vo ()- 


—— 


(" cr (ky + 1) 0% VP ED yr 4) 
ad Oe (2.34) 
v2 (W) ¥3 (0) 


It can be seen from Eq. (2.32) and expressions (2.33) and (2.34) that the calculations be- 
come appreciably simpler when Pr..=0.5. The equations given above allow us to obtain the 
variations of velocity, density, enthalpy, and temperature across the initial section of the jet. 


Sample Calculations 


1. Let us consider the initial region of a plane-parallel turbulent submerged jet of nitro- 
gen at a pressure of 392 10‘ N/m? and a temperature T)=75°K. 


The propagation of the jet occurs in a medium consisting of nitrogen. 


The turbulent Prandtl number is equal to 0.5 (k, =0.5, k, =0) and the coefficient ais equal 
to 0.07. 


The coefficients appearing in the equation of state are: A=94.74, B=1017.28, C=106.99, 
=—1361.23. 


The approximate equation of state is shown in Fig. 3. 


The parameters of nitrogen on exit from the nozzle are: T)=75°K, V)=0.124-107? m3/N 
ip =45.8-10° J/kg. 


The parameters of the Surrounding medium are: 


1" 3.64, V,=2.065-10-3 mYn, i,=436-10" J/kg; 


9) 


1. 
—" = 6,31, V,=3,57-10-3 mYn, i, = 642-10" J/kg; 


_ "8.97, V,=4.8-10-3 m¥n, i, = 816-108 J/kg. 


9) 


Selecting a number of values of V in the range Vy to Vy, we determine 9= V/Vp. 


Substituting the values of 9 into the energy equation (2.32), we find the corresponding 
values of 7. The profile of the relative specific volume is shown graphically in Fig. 7. 


Making use of the graphs of Figs. 2 and 4, we find the values of i and T corresponding to 
the selected values of V. 


We now determine the profiles of + =i/ig (Fig. 8) and y =T/Ty (Fig. 9). Using Eqs. (2.33) 
and (2.34). we calculate the functions 0'/o (Fig. 10) and 9"/06 and determine the function o(n) = 
ho"/a—(60)* (Fig. 11). 


The results of all calculations are given in Table 3. In the graphs of 0'/9=f(n) and o(n) 
(Figs. 10 and 11), we subdivide the interval of 7 (0 — 1) into r equal subdivisions of width 1/r. 
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ig. 7. The profiles of the specific Fig. 8. Enthalpy profiles in 
volume in the boundary layer of the the boundary layer of the 
jet. 


ject. 


Fig. 10. Profiles of the func- 
tion 0'/o in the boundary lay- 
er of the jet. 


Fig. 9. Temperature profiles 
in the initial section of the jet. 


Fig.11. Profiles of the func- 


tion a(n) in the boundary layer 


of the jet. 
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The formula for the evaluation of the integral J, by the 


I, 


| 
—_— | >.” k —-j_ 7 R 
r [3% oh 


_| 


trapezoidal-rule method is of the form 


‘iy k fy" k 
—}[— } qt (—) qe te. + 
fj | f) 9 


while the corresponding formula for the integral j,, is 


fy’ h ' ~ fy’ k 
a tp 174 0.5 { —\ 4, 
ly r jl t) r . 


(2.35) 


op h ~ 
Ip 71 9.53,4% |. 


(2.36) 


The integral J); has been calculated with r—10, the 


integral j, with r=20. 


The values of the integrals are given in Table 4. 


Having determined the coefficients A,, Ao, 
... (from Eqs. (2.21)] for the second approximation, we 
solve the system of equations (2.20), find the valucs of hz 
and ‘a, and determine q, and qg from Eq. (2.22). 


veny Ps 


The results of the calculations are given in Tahle 5. 


The velocity profile (Fig. 12) is caleulated from 
Eq. (2.19). 


From Eys. (2.11) and (2.19) we obtain a formula for 
the determination of the total head in terms of the independent variahle n, 


i 


——— 


2 


( 


t 
fe 1 
2 


| 


t 


3 (1 —m), 


(2.37) 


We now determine the profile of the total thrust with the condition that uyy=0 and m=0 
(Fig. 13). The results of the caleulations of the velocity and total-head profiles are given in 


Tahle 6. 


TABLE 4 


2 
—),5573 
0.0149 
—().7819 
—().3257 
—().927] 
—(). 6863 


7 
—0, 1257 
(),0774 
—(). 1676 
(),1168 
—(), 1988 
0.144] 


QO 


— 03480 
0,125] 
—().4938 
() 0228 
—() 5938 
—(),1251 


—(), JO86 
0.0656 
— 0.1428 
(0.1018 
— 0.168) 


4 
— 03,2443 
(), 1287 
—(),3435 
(), 1227 
—().4)47 
0.0767 


5 6 
—(), ]860 —(),]498 
Q.Jd11 (),0926 
—(),2570 —(), 2032 
0.1403 0.1319 
—0,30903 | —0.2430 
(,1409 (), 1522 
Continuation 
Sg) 
— 0.0960 Tw 3 64 
(0566 f To 
—(), 124 T 
), 1248 l r 6.3] 
() 0887 J Ty 
—(), 14548 , Tug 97 
0.1119 if To 


9.1302 
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TABLE 5 
mp t 
Ty ao * F1 72 
3.61 ne 2.26 0.728 | 1.539 
6.31 ' 1.93 1.945 0.566 | —1.379 
8.27 | 2.35 1.695 0.452 | —1.243 


As a result of all this, we obtain the variations of the main parameters across the boun- 
dary layer of a submerged jet of nitrogen at a pressure of p=392 104 N/m’ at several tempera- 
tures of the surrounding medium. 


2. Let us consider the initial region of a turbulent plane-parallel jet of incompressible 
fluid in a medium moving in the same direction. In this case 9=1. 


We select several values of the quantity m=uyj/u) in the range 0 to 1. 


Solving the system of equations (2.20), we determine the values of DB , by, and gm, asa 
function of m (Table 7). 


The results of the calculation of velocity profiles in the boundary layer of a jet for values 
of the relative velocity of the parallel current of m=0, m=0.3, and m=0.9 are given in Table 8 
and are plotted graphically in the form of curves of Z=f(n), Z=f(@), and by [Dx wy =o =F (m™) in 
Figs. 14, 15, and 16. 


For comparison purposes, let us plot on the graph of bx /bx ,, —)=f(m) the analogous de- 
pendences obtained by G. N. Abramovich [17]: 


bo. _ il—m 
Dem 0 —Y Lm 
z 
] 
05 
= 
| 
ae 
0 05 
Fig. 12. Velocity profiles in Fig. 13. Total-head profiles 


the boundary layer of the jet. in the boundary layer of the jet. 
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TABLE 6 


Z | 0 0.0372 | 0.1332 | 0.2663.) 0.4177 | 0.5714 | 0.7137 | 0.8343 0.9252 | 0.9812 H | 
— - | 
CC, 0 0.0002 | 0.0032 | 0.0216 | 0.070 | 0.1782) 0.3277] 0.5142 0.7112 | 0.8793 l 3.64 
00° Ho 
% —1.5317 '—1.3057 |—1.0797 '|—0.8537 |—0.6277 |—0.4017 |~0.1757 | 0.0503 0.2763! 0.5023! 0.7283 
| | | | 
| 
Z | 0 0.0436 | 0.1534 | 0.3012 | 0.463 0.6205 | 0.7590 0.2602 | 0.9454 | 0.9876 ] 
aes 0 0.0002 | 0,0019 ] 0.0131 | 0.0479 | 0.1479 | 0.3215 |; 0.5260 | 0.7467 | 0.8923 ] 6.31 
Q0° 49 
% —1.3793 ees —0,9903 '~-0.7953 '—0.601 '—0.4073 '~-0.2123 pos me 0.3701 a 
Z 0 | 0.0470 | 0.1642 | 0.3196 | 0.4874] 0.6469 0.7837 | 0.8876 | 0.9562 | 0.9910 1 
Queu () 0.0002 | 0.0014 | 0.0077 | 0.0884 | 0.1207 0.2865 | 0.5017 | 0.7679 | 0.8987 ] 8.27 
Qo: Lo 
© -- 1.2425 !}—1.0730 |~-0.9035 '~.0,7340 '—0,5645 |—0,3050 ~-0,2255 (— 0.0560 | 0.1135 | 0.2830 | 0.4525 


TABLE 7 
m 0 0.15 0.3 0.5 0.7 0.9 
Qo —1,532 | —1.0 | —0.705 | —0.43 | —0.22 0 
b, 3.037 | 2.67 2.37 2.01 1.62 1.08 
91 0.983 | 0.9551 0.92 0,856 | 0.736] 0.525 


= 
= 
a 
a 
a 
= 
a 
Y 
a 
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Fig. 14. Profiles of the excess Fig. 15. Variation of the excess ve- 
velocity in the boundary layer of locity in the boundary layer of a jet 
a jet ofincompressible fluid. along the coordinate g. 
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TABLE 8 
7 0 0.1 0.2 0.3 0.4 0.5 
Z 0 0,0156 (). 0648 0.1484 (). 2637 0, 4042 
© —2.054 | —1.75 | —1.446 —1.1441 —0.838 —0.534 
Z 0 0.0223} 0.086 0.185 0.3114 0.456 
| 
9 —1.454 1.217 —0.979 —0,742 —0,505 + —0.269 
Z , O 0.028 0.104 0.216 © 0,352 ' 0.5 
2 0.555 OT 0.339 0.281 0.123. | 0.015 
Continuation 
r 0.6 07 ' os | 0.9 1.0 
] | 
Z 0.56 , 0.7165 0.856 0.961 1.0 ; 
| . m= 
p 0,23 0.07 0.376 0.686 0.983 J 
Z 0.6074 0.753 ().878 0.966 1.0 ) 0.3 
MEV oc 
9 —0.032: 0.206 0.444 0.686 | 0.92 J 
Z 0.648 0.784 0.896 0.97 1.0 | ) 
m=0.9 
¢ ().093 0.201 0.309 0.418 0.525 + J 
| | 


(Fig. 16) that the relative width of the boundary layer as 
boundary layer of a jet of incompressible fluid with 
the velocity of the surrounding medium. 


and O. V. Yakovlevskii [19]: 


b l 


* 


bameg 1 +2-36m 


It can be seen from the graph (See Fig. 16) that the relative width of the boundary layer as 
obtained by us is greater than that yielded by the other theories. As the velocity of the parallel 
flow increases, the velocity profiles becomes fuller (Figs. 14, 15). 

3. The Main Region of a Plane-~Parallel Turbulent Jet 


of a Real Gas 


Formulation of the Problem. Let us consider the case of the discharge of a jet 
of real gas from a long and infinitely narrow slit into a parallel stream (Fig. 17). We will place 
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the coordinate origin at the point of emergence of the jet from 
the slit, in other words, at the source singularity [17]. The 
velocity on the jet axis will be denoted by u,,, the density by 

Pm: and the temperature by Ty. These parameters will vary 
along the length of the jet. The velocity of the surrounding 
medium uy, density py, and temperature Ty are constant param- 
eters along the length of the jet. 


The transition section of the jet is taken to be the trans- 
verse section where the central velocity u,, is equal to the 
velocity uy, at which the jet issues from the source. We will 
match the solutions for the initial main regions of the jet across 
this transverse Section. 


Fig. 17. Turbulent source. 


The emerging jet in mixing with the surrounding medium forms a boundary layer whose 
width will gradually increase in the direction of motion of the jet. The variation of the param- 
eters across the mixing zone will be smooth and they will vary from their values on the jet 
boundary 1 to their values on the jet axis. The width of the boundary layer is taken to be the 
distance from the jet axis to the boundary. 


Let us subdivide the jet into an infinitely large number of intervals Ax. We will assume 
that over the interval Ax the density and the value of the ratio m=uj,/Uy, are constants. We 
will solve the problem for an arbitrary interval Ax. 


We will make use of Eqs. (2.3), (2.5), (2.9), and (1.1) for this purpose, as well as the 
equation of conservation of momentum 


Vv 


boo 
( att (Ut —,) dy! =, (9 — Uy) b, (3.1) 
P 


and the equation of conservation of excess enthalpy 


‘b 


\ pu (§—4,) dy = bgtly (ig — iy) Oo, (3.2) 
0 


where 2b, is the width of the initial section of the actual jet. 


Transformation of the Equations. Let us introduce new variables x, 9 =y/ax, 


where a = V2c2 is an experimental constant, for the solution of the system of equations (2.3), (2.9), 
(2.9), (1.1), (8.1), and (3.2). 


We will introduce dimensionless velocities, current density, gas density, specific vol- 
ume, enthalpy, and temperature defined by 


at 7), HO =F (2), + =9(9), 
Um — Uy emiUm Om 


VY _6(s), =3(¢), ——- == 7 (%). (3 .3) 


ry e 
Vin lm 


In the interval Ax (where we assume p,, =const. and m=const.), these functions will be 
independent of the variable Ax. 
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Let us transform the equation of the conservation of excess momentum (3.1), making use 
of the functions (3.3) and introducing the variable m. We obtain the equation 


Tb 
AX (Um Uy) Fe) Z (2) de> sth (ty Uy) dy. (3.4) 
0 


From Ea. (3.4) we find that 


—_—— ee 


M M — 


m —_—_-—______ — 


}d—m J] x } x 


; (3.5) 
where 


poll, (1; — 3) 5, ] 


M= _ ,; - 
/ *p ., am 
\ F’ (9) Z (2) dz 


ww 


0 


Let us express the current in terms of the principal component of the current density 


= \ ody = AX gl m \ F' (2) de. 


. 


Substituting into this equation the expression for um from (3.5), we obtain 


Af 


e 
¢ ——_ 
J=AXG,, 


F (9). (3.6) 


} (1 — m) } x 
Using the continuity equation, We find the transverse current-density component to be 


(0 = — ap, ——— | oF) SF 0) (3.7) 


This expression is correct only when pm and m are constants. 


Let us write down the equation of motion (2.3) as follows: 


— O(u — ty) a 0 (u — Uy) 2.9 7) (p [2 if 


ott —————_——__ + 6U $i 0 OX . 
J { dy J 


Ox Oy Oy (3.8) 
Substituting into this equation expressions (3.3), (3.5), and (3.7) and taking into account 


that 
Od (i — uy) _ l d (u — ty) 0 (u — uy) FP d (tt — ty) 
oy ax dz Ox i x de 
we obtain the equation of motion in the form 
FZ=(1—m)v(2'/. (3.9) 


We can eliminate F between Eqs. (3.9) and (2.11) to obtain the equation of motion in the 
following form: 


M74 Zt g A 8 gg <—(Z'°Z—222'2" 1.(Z')*=0. 


]— om 
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Replacing in this equation the density ° by the specific volume (< =1/9) and transforming to a 
new independent variable 7 (n = ~/bsx), we obtain the equation of motion for the boundary layer 
of a plane-parallel turbulent jet valid in an infinitesimal interval Ax (p,,=const. and m=const.), 


- »73 013 n {y’ 
L(Z) =~ --— £70, 2*b,- a - b, — 


|—m 2 fj 


BL "ZZ ULL AZ, (3.10) 
where bx is the dimensionless width of the boundary layer. In solving these equations, we use 
the following boundary conditions: 

1. On the jet axis with 7 =0 (g=0) 


a) the velocity is equal to u,, and, consequently, 


Z(0)=1; (3.11) 


b) the transverse component of the current density must be zero on the jet axis (pV=0), 
so that 


Z’(0) =0. (3.11a) 


2. On the boundary of the jet where 7 =1 (y= gp) the velocity is equal to uy and, conse- 
quently, 


Z(1) =0. (3.11b) 


Let us investigate the behavior of Eq. (3.10) on the boundaries. In order for the solution 
of this equation to satisfy the boundary conditions on the axis, the following relation must hold 
when 7 = 0: 

3 


b 
——*_ = Z' (0) Z” (0). (3.12) 


l1—w» 


Since Z'(0) = 0, whereas the left-hand side of expression (3.12) is not zero, then Z"(0) ~©. 
On the jet boundary when 7 =1, the solution of Eq. (3.10) will satisfy boundary condition 
(3.11) provided that 
Z°(1) =0. (3.12a) 


The function 9'/6@ appears in Eq. (3.10) and the form of this function will be determined 
when the equations of motion and energy are solved simultaneously. 


Let us transform the equation of excess enthalpy (3.2) making use of functions (3.3), the 
variable g, and the dimensionless enthalpy 


< (2) == isin ME 
im — bu Ain 


The equation of conservation of excess enthalpy will then become 


b 
AXt gl mA | eF'dE= oly (iy — by) Bo. (3.13) 
0 
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Equation (3.13) is now used to find Aiy;: 


_. rll) (4), ~~ in) by, ] 


2b ax 
omlUn | eF'd> 
b 


Substituting into this equation the expression for um (3.5), we obtain 


Ai,=—- ] d—m) —— , (3.14) 
where 


pz Gilly, (fy - bi) by 


*b 
sma \ eF'ds 
0 
Let us write the energy equation (2.9) in the following form: 


— O(i--i,) , — dG-7 
oy Oy OE i) 


— = — NCX* —~— 


OX Oy ay 


0 E O(u—uty,) Oi --in) |. (3.15) 


Oy Ou 
where 
n==1/Pr,. 


In the case of appreciably subsonic flow, the last term in Fq. (2.9) can be neglected since it is 
small by comparison with the other terms of Eq. (3.15). Substituting relations (3.3), (3.5), (3.7), 
and (3.14) into Fq. (3.15), we obtain the energy equation in the following form: 


Fe=n(1—m) 02’e’. (3.15a) 
Dividing the energy equation (3.15) by the equation of motion (3.9), we obtain an energy 
equation that is valid for both variable wand variable n, 


r 
£ 


Zz 
Pr, = 
Z 


r (3.16) 


ps 
cS 


with the boundary conditions 
a) on the jet axis p=0 (n=0); Z(0)=1. The enthalpy on the jet axis is equal to i,, and, 
consequently, 
e (0) =]; - (3.17) 


b) on the jet boundary y= gp(7 =1); Z(@p) =0. The enthalpy on the boundary is equal to 
izjj and, consequently, 


E (Pp) =0. (3.17a) 


Making use of the boundary conditions, we will obtain the solution of the energy equation 
(3.16) in the form 


eZ", (3.18) 
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Solution of the Equation of Motion. We will solve the equation of motion (3.10) 
with the main boundary conditions (3.11), (3.1la), and (3.11b) and supplementary conditions (3.12) 
and (3.12a) for the general case, by Galerkin's method [30] 


i oR 
ZyN ajo; (-=0, 1, 2,3,...8). (3.19) 


} 0 
The system of functions 9} must satisfy the main and supplementary conditions on the 
boundaries. 


Let us consider the following system of functions: 


~~ 


yh g=(4? — 24"); oo=(1 42)? yPlae 
3(&—1) 


O,== (L—7f)?x3, .. Op=(L—y)2y 7 (3.19a) 


If we take @)=1 and @,=1, then the sum of functions (3.19a) satisfies the boundary condi- 
tions for any values of the coefficients @,,@3, ...,@,. 

We will seek the solution of the equation of motion (3.10) with the help of the first three 
functions of system (3.19a). The solution will then be of the form 


Zy= 1+ (M9 —2n") +441 — YP). (3.20) 


In this case, @, and b, = g, Will be the unknowns to be found. 
Using Galerkin's method, we will write down two equations for the determination of these 


unknowns 


1 1 
\ L(Z.)¢an=0; \L(Z,)%d4=0. (3.21) 
0 


0 


Transforming Eqs. (3.21), we obtain two equations for the determination of the unknowns 
ao and Dy 


Ai + A3.-6°B}—6,B3=0; Ay+ Aj+02B)—6,B\=0. (3.22) 


In the first equation of system (3.22), the functions A* and B* are given by 


Aj = 0.58478 + 0.38136a, + 0.24747a3 + 0,01891 a3; 
Ag== A,-+ 0,43 4243; + a2A;; 


Bi = (0.0394 + 0,065 ——- 4 (0.02167-+ 0.0307 | @o-+ 
l—am !—m (3.23) 


M1 


4. (0.00574 +-0.0051 +) a3 + 0,00064a3: 


B= Q > (B, -aB; +a3B, +03B.), 
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where 
AL, = 9S — 54S, -+ 1267, — 144 J/u,+ 814, - 1845,; 


Ay= 36Je., — 225J,, + 558Je,— 702J; + 468Ju,- 153, + 1843 


A. =—— J 3 _ 306S:5:, -|- ji Sy) — “ Ju, + ee J. — 3833/0, a J, ‘ Js .: 
A= SL Jan, 135 Jy Wr Jen. en Jig ce Jn, —198J, $< Inj, — Sys 
By = 3d, - 2104-604, — 90S, + 75d, -- 332, + 6; (3.23a) 
By Ju, -- 75d, Oh te — 345d, + Oe J, = 153, | 2, — 38d: 
2 : 2 2 2 
B,=12),, —87Ja + 261J,, — 420J» ,-+ 390, — 207, +57), — be, 
Bo= ~ Jor, — 334,44 " Js, —165J,, + “= Jn, 87S, +— In. — 3d, 
In the second equation of system (3.22), the functions A* and B* are given by 
Aj= — (0.19612 +. 0,06797a, -+ 0.0514a3 — 0.00023.a3); 
Ai=(A,+a,A,+ aA, + a3A,); 
* [ m m ~ 
>= |(0.01458+ 0.0218 —— ) + f0,00683 +..0.0084 
B} (0.0 + 0.021 4 (0 00683 +0 i) et 3.24 
+ (0.00155 +-0.00117 —" ) + 0.00015@) ; 
— mm 
By=a ~~ [By + a,B, + 0B, + Bul, 
where 
Ay =a ,— BAS ;, + 135Je,, —- 180J, 4- 135Ju), — 54, + Ws, 
Ay = 36J;4 — 225Ju), +594 J, — 855Jir,,-+ 720J; — 351Ju, + 90), — 9Js,,5 
Ay =a Js, — 306J 15 + 846 J), — 1296, + — J, — 66S; +216 Ju, — 36 , + = Jes 
81 270 - 9). 
Ag Sig a Sj, $387 jy — 621 Sin, + 6075S po — 3869S, 135Sq — 27 ys I 
= BJs, — 21S, + 63/0, — 1054 + 105Ju0., — 63/, “+213, — 3/3 (3.24a) 
B,= — Jy, — 75d, 2314, ~399J 0, + 420J, — 273 Je, +105; — 21s, + = Jy; 
By = 12S, — 87 pat 273), — 483Jq, +525 Je, — 357 9+ 147 Ju;, — 33) ,-+ 3/2, 
By=— Syy — 33Soy, 4-105) 4 — 185 ),-+ 210S,, — 14 7S, + 68g — 15, as J. 
The quantity J, appearing in these functions is defined by 
1 0" 
Iy=\ afd, (R=1, Ys Fy 3d) (3.25) 


In order to solve the first equation of system (3.22) for b,, we will write it in the form 
[33] 
b°4+-3P,6, +-2Q,=0 (3.26) 


V. I. BAKULEV 65 


where 


om 


Ay +- A, 
3 Bt 2By 


Since according to its physical meaning by is a real quantity then, rejecting the irrational 
roots, we obtain the solution of Eq. (3.26) in the following form: 


b= Nay Va, (3 .26a) 


where 


LQ + VO+ PIs Ve =V —AaA-) Qi+ PI. 


Wy == J 


The second equation of system (3.22) is solved in an analogous manner. Let us write 
down the second equation of system (3.22) in a manner analogous to (3.26) 


3 
b.--3P.6, -|- 20,=0, (3 27) 
where 
1 2B At 4 A" 
P,= — — : , QO= oe ; 
3 BS 285 


The real root of this equation will be 


6, =Uygg + Vxr (3.27a) 


where 
Vf We Von Wao 
go} —Q.4+)/ Q3-+ P3 ) Viz | —~Q,—-V @+P3 . 


The system of equations (3.26) and (3.27) is solved graphically. The point of intersection of two 
curves gives the solution of the system. 


The Calculation of a Jet of Real Gas. Let us specify the value of the density 
Pm (or the specific volume Vy) on the jet axis and let us assume that the axial density does not 
change over the interval Ax. In the first approximation, we assume that the relative specific 
volume can be described by 


i= Va = 6 (3.28) 
Vin 
In this case 

——=In §,,= const. (3 .2 8a) 


In the first equation of system (3.22), the functions A¥ and Bf of the first approximation 
will be 


At =—In 0,(0.21273+0.04485a, -+ 0.02243 + 0,00397a3), 


~ ~9 ~ (3.29) 
B3=a > 1n 0, (0.08335 -+ 0.027792,-+ 00055623 + 0.000523), 
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While the functions A* and Bf of the second approximation will be 


A= In 9, (0.0621 +-0.003a3 + 0.00068a3), 
Bi= a In 4, (0.02778 +-0.005570, + 0.0095 1a). (3.30) 


We will solve the system of two equations (3.26) and (3.27), making use of the expressions 
(3.23) and (3.24) for Aj, Bf» Ag > Bs , and expressions (3.29) and (3.30) for As, Bs, Aq, B;. 


We first find @ and b, . 


From Eq. (3.20) we determine the velocity profile Z=f(7) in the boundary layer of the jet 
corresponding to the specific-volume profile adopted in the first approximation. Then, we de- 
termine the enthalpy profile from the solution of the energy equation (3.18). The dimensionless 
enthalpy profile can be used to calculate the enthalpy i at any point of the transverse section of 
the boundary layer 


ise (im—in) Hin. (3.31) 


In this equation, the enthalpy of the parallel flow i; is known, while i,, is found from the 
graph of i=/f(V) since py(V)m is known. 


Knowing the enthalpy at any point on a transverse Section of the boundary layer, we make 
use of the graph of i=f(V) to determine the specific volume V and to find 9 =V/Vin> i.e., the 
new specific-volume profile across the transverse Section of the boundary layer. 


Next, we write down the equation of state (1.1) in the following form: 


(in — by) tiy= AV, + Bf-S- 42 | (3.32) 
V m4 ve ne 
Differentiating Eq. (3.32), we determine the function 9'/9 
ee e’ (Lyn — iy) (3 3) 
= Oo 
; AV, — C 2D 
V m9 v7 02 
In this equation, the derivative with respect to cis 
3! — Pr Zen zZ, (3.34) 


Making use of the expressions (3.33) and (3.34), we can find the function 9'/9=f(n) and cal- 
culate the integrals (3.25) in the following manner. On the graph of 9/%@=f/(n), we subdivide the 
range of 7(0—1) into an equal number r of intervals each of length 1/r. 


The formula for the calculation of the integral J), by the trapezoidal method is 


__ | ’ re, { W k (8 _k ~({W\ k 
1, — (—-) of nie (--) —L eee (+), Mr-1 +0.5 (+), Vir | . (3.35) 


We now determine the functions A,, Ag, ... and B,, By,.. .[Eqs. (3.23a) and (3.24a)] and, 
solving the system of two equations (3.26) and (3.27), we determine the new values of the quanti- 
ties a, and b, . 


In this manner we determine the variation of the parameters across the jet. It now re- 
mains for us to calculate the variation of the parameters along the jet. 
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Let us write down the equation of conservation of excess momentum (3.4) in the following 


manner 
1 — m, | J ay 
im — Vee Ve (3.36) 
—— 


Uy, 


where my) =uyq/Uy and by are parameters evaluated at the jet nozzle. 
In this equation, F' is 


Fl =——[m+(1—m)Z b, a (lm) 2). (3.37) 


* 


The equation of conservation of excess enthalpy can be written as follows: 


(3.38) 


We determine the magnitude of the parameter ax/b, from Eq. (3.38) and, substituting its 
value into Eq. (3.26), we find the relative velocity u,;,/u) on the jet axis. 


If the jet is propagating in a parallel stream, then the calculation has to be carried out 
by a trial-and-error method with m the adjustable parameter. 


In order to match the initial and main section of the jet, we calculate the transitional 
transverse Section of the jet (pm =p), Aip, =Aiy). The relative velocity up)/uy on the axis in the 
transitional section (py, =/,) is found to be greater than unity, i.e., the transitional section of 
scalar quantities lies closer to the nozzle than that for vector quantities. The intersection of 
the curve Up)/uUy) =f(@x/by) with the line um/uy)=1 gives the coordinate of the dynamic transitional 
section. 


The calculations are performed for several transverse Sections. 


In this way, we can obtain the variation of parameters along the length of the jet. 


Sample Calculations 


1. Incompressible Jet. Let us compare the approximate solution with the exact Prandtl— 
Tolmin solution for the main region of a plane-parallel jet [17]. 


Using Eqs. (3.26) and (3.27), we find the values a,=—0.5, by =2.477. 
The velocity profile is calculated from Fq. (3.20). 
The results of the calculation are given in Table 9. 


For comparison, Table 9 also contains the Prandtl— Tolmin data on the velocity profile 
F'. The dimensionless width of the boundary layer according to the Prandtl—Tolmin theory is 
Dx =2.4. 
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TABLE 9 

< () 0.2 0.9 0.8 1.0 1.2 

Z 1.0 0.9-4] 0.792 (), 607 0.475 0.36 

Ff’ 1.0 0.94 0,782 0.604 0.474 0.357 
Continuation 

' | 

° 1.5 1.7 2 2.2 2.4 2.477 

Z 0,21 0.13 0.00 0.022 0.0005 () 

fF’ 0.2 0.125 ().046 0.02 0) 8) 


Let us calculate the integral appearing in Eq. (3.26) 


1 
b, \- F'Zd1, 0.695. 


0 


Substituting the numerical value of this integral into Eq. (3.26), we obtain an expression 
for the axial velocity of a plane-parallel jet 


Um 1.2 


, (3.39) 


analogous to the Prandtl—Tolmin formula [17]. It can be seen from Table 9 and formula (8.39) 
that the agreement of our solution with the Prandtl—Tolmin solution is satisfactory. 


2. Main Region of a Plane-Parallel Submerged Jet of carbon dioxide at a pressure of 
p=981 -104 N/m? and temperature T)=248°K. 


The parameters of the surrounding medium are: T,,=908°K, V}y,=1.82 -1073 m°/N, 


if7=137-104 J/kg. 


In the calculations, the turbulent Prandtl number is taken equal to 0.5 and the coefficient 
ais taken to be 0.1. Let us take several values of the specific volume Vy, on the jet axis. We 
will determine the function 6'/9 by means of Eq. (3.28a). 


We first of all solve the system of equations (3.26) and (3.27). Todo this, we take a 
series of values of ay and determine the values of A* , AZ, A* Ay and BF , BF , BF ; B* 
from Eqs. (3.29) and (3.30) for each value of dy selected. 


We then construct a graph of bs =f@,) for Eqs. (3.26) and (3.27) and from the intersection 
of the two curves we obtain the solution of the system of equations. 


The results of the calculations are given in Table 10. 


TABLE 10 
ee 
Section I 9.5-10 4 18.8 90 248 0.7 1.44 
Section II 12.5-10 4 14.28 120 298 0.7 1.6 
Section III €0-10 -4 2.98 180 388 0 ] .64 
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TABLE 11 
y, 0) 0,1 (),2 0.3 O.4 0.5 0.6 ().7 0,8 0,9 1.0 
Z ] 0,958 0,882 0.779 ().657 0,521 0,38] 0.242 0.123 ().037 () 
_, e : ] 0.979 0.939 0.883 0.810 0.722 0.617 0.492 ().35 0,192 0 
5 0 | l 1.05 1.105 1.37 1.9 3,79 G.42 9.05 12.0 15.3 18.8 
S| gui | i 
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Using Eq. (3.26), we calculate the velocity profile (see Fig. 18), 


and from Fq. (3.18) we 


calculate the enthalpy profile (see Fig. 19). Making use of Eq. (3.31), we find the value of the 
enthalpy for the given transverse section and from the graph of Fig. 2 we determine the value 


of the specific volume. 


We now find the improved specific-volume profile (Fig. 20) 


Incompre s- 
sible fluid 


Fig. 18. Velocity profiles in 


the main region of a flat jet. 
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Fig. 19. Enthalpy profiles in 
the main region of a flat jet. 
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Fig. 20. The specific-volume 
profiles in the main region of 
j | } Incompres- 
a flat jet. ) — sible fluid 
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Incompres-_ 
sible fluid Fig.22. Variations of velocity and 


enthalpy along the axis ofa flat jet. 


Making use of the equation of state (1.1), we find 
the function 9'/g from Eq. (3.33). Using Eqs. (3.20), 
(3.37), and (3.18), and (3.33), we determine total thrust 
and enthalpy profiles (see Fig. 21) 


Fig. 23. Variation of the relative juste (yy _ Pea 
width of a flat jet along its length. Pm Hn omimrim 


=f (%). 


The results of the calculation are given in Table 11. 
By a graphical method [32] we evaluate the integrals 


( : F'\Zdy and ((- P’)edy, 
; 


x 


then, using Eqs. (3.38) and (3.36), we determine ax/by and uy,/Ug. The results of these calcula- 
tions are given in Table 12. Figure 22 shows the fall of the relative velocity u;,/up and the rela- 
tive enthalpy Ai,,/Aiy along the jet axis as functions of the coordinate ax/b). The variation of 
relative velocity along the axis of an incompressible-fluid jet is shown in Fig. 22 for com- 
parison purposes. 
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TABLE 12 
Section I Section II Section III 
So l 1.32 6.3 
Om 
¥ 
as ().872 0.62 
Alp 
7 
\(5- F’) Zar 0.3218 ().278 0.279 
Q * 
1 
\( Ved 0.354 (301 0.317 
ae 1.78 3,22 97.9 
bg | 
Lm 
— 1.1 0.945 0.703 


Figure 23 shows the relative width of the boundary layer of a jet of carbon dioxide and the 
relative width of a jet of incompressible fluid. 


It can be Seen from Figs. 18 and 19 that as the density ratio py,/p) increases, the velocity 
and relative-enthalpy profiles become more convex. The decrease in the velocity along the 
axis of a jet of carbon dioxide is less rapid than that in thecase of a jet of incompressible fluid 
(see Fig. 22), while the width of the boundary layer increases more slowly in the case of a jet 
of carbon monoxide than in the case of an incompressible fluid (see Fig. 23). 


4. The Main Region of an Axially Symmetric Turbulent Jet 
of Real Gas 


Formulation of the Problem. Let us consider the case of the outflow of a jet 
of a real gas from a turbulent point source (See Fig. 17). As in the case of the main region of 
a plane-parallel jet, we place the coordinate origin at the source singularity [17]. We will use 
the same notation as in the case of the plane-parallel motion. We will match the initial and 
main regions of the jet in the transitional transverse section. 


We will subdivide the jet into intervals Ax along its length and we will assume that the 
density remains constant (p,,=const.) in each of these intervals. We also assume that m=uy/um 
is constant in Ax. 


We will make use of the following equations to solve the problem: 


1. Equation of motion 


— ou — ou >» 1 @ —/ Ou ) | 
tl —~ + VU = — C°X- — —— | Ve | —— : 4.1 


2. Continuity equation 


ou, Oy WY _g (4.2) 
OX T oy + y 
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3. Energy equation 


— Oi »»o 1 @O —~ ou Oi 
| 0 |: (4.3) 


U—— == — NOX 
+ pv — ay Oy 


— di 
Olt 
Ox oy yoy 


The averaging of these equations was carried out according to Reynold's method in the 
same way as in the case of plane-parallel flow. 
4. Equation of conservation of excess momentum 


b—_— — 9 
Qn \ pu (u—u,) ydy=o,u, (u,—U,) aR. (4.4) 


a 
OC ye 


5. Equation of conservation of excess enthalpy 


Yb 
Qe! ot E—iy) YY = pylly (ip — by) =Ro. (4.5) 


0 


6. Equation of state (1.1). 


Transformation of the Equations. Let ustransform Eqs. (4.1), (4.2), (4.3), 
(4.4), (4.5) and (1.1) into the coordinate system x, g=y/ax with the help of the functions (3.3). 


3). , 
Here,a=Vc* is an experimental constant. 


Let us transform the equation of conservation of excess momentum (4.4) making use of 
functions (3.3) and the new variable g~. We will obtain this equation in the following form: 


Pol (fy — Uy) \ Ro 


aab 
frat (im — tn) (“J | s@2Z@)ede=05. (4.6) 
0 
From Eq. (4.6) we can obtain the values of u,, and Aum =Um-— Up 


M MVvVi— 
Ll. = Au, =- V J (4.7) 


x 


(y44( — bby R On 
M= / ae tl) Polly (Ho Hu) Ro _ f= pu 
aV om Omim 
> ” 5(e) 209): F429 


Using the continuity equation (4.2), we can express the current density components in 
terms of the current function 2: 


where 


Oy. Loy (4.8) 


1 
ell — U 
yoy y Ox 


From the first equation of system (4.8), we find the current function to be 


ve ( Pld WAY = Pll VX? \ I (¢) d¢. 
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Replacing in this expression uy, by its value (4.7) and denoting the integral by F(g), we obtain 


Al 
——_——. a-xF (9). (4.9) 


Y= 0m 
} l—wm 


From Eq. (4.9) we then find the main component of the current density 


| F' (9) 
——. (4.10) 


Using the second equation of system (4.8) and Eq. (4.9), we find that the transverse component 
of the current density is 


70 =f, —a ~|P@-2 |. (4.11) 


This expression is valid in an interval Ax where Py, and m are independent of x. 


Let us write down the equation of motion (4.1) in the form 


— O(u—U — oO(u—u 
1 Oke ae) Oe a)? 


19 (=f _O(u—ay) 7°) 4.12 
Ox ‘ Oy y dy \>! | | . ( ° ) 


oy J 


We can transform this equation with the help of expressions (3.3) (4.10), and (4.11), as 
well as the expressions describing the equalities of the derivatives of the velocity with respect 
to x and y, 


0 (u — uy) __ ? d(u — Uy) 
a FE 


Ox x dy 


O(u—uy,) | d (u — Uy) 
oy ax dp ) 


As a result of this, the equation of motion will become 


= =(1—-m)0(Z’)?, (4.13) 


The relation between the velocity component u and the current density pu can be found 
from 


Pll = pll + ait’. 


Substituting this expression into Eqs. (3.3) and (4.10) and making use of the Prandtl] turbulence 
model, we obtain an equation connecting the functions Z and F 


f= 8 [m+ (1 —m) Z] 4 nad'Z'(1—m). (4.14) 
» 
Eliminating F from Eqs. (4.13) and (4.14), we obtain the equation of motion with respect to the 


variable o 


224. 28 (ZZ —22"Z'Z (ZY —(Z')'Z — na — 2/22 =0. (4.15) 
u ¥ 


l—m 


Replacing in the equation of motion (4.15) the current density function ¥ by the specific volume 
function 6(% = 1/9) and introducing the independent variable n(n = o/b,), we obtain an equation 
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of motion which is valid over an infinitesimally small section Ax where We can assume that 
Pr =const. and m=const. 


L (Z)=—— Z°b, + 256, +— (Z'PZ—2Z"Z'Z4+(Z')" —Z(Z'’ = — na — Z'27b, =0. (4.16) 


In solving this equation, we will make uSe of three boundary conditions analogous to condi- 
tions (3.11), 


1. On the jet axis 7 =0(g=0); 


Z (0) =1, (4.17) 
2. On the jet boundary 7 =1(gp=Dx«); 
Z(1) =90. (4.17b) 


Let us investigate the equation of motion (4.16) on the boundaries. In order for the solu- 
tion of the equation of motion to satisfy the boundary conditions on the axis, the following equal- 
ity must hold on the jet axis: 


3 


b 
; — = 27" (0) Z” (0). (4.18) 


Since Z'(0) =0 on the jet axis, while the left-hand side of equality (4.18) is not equal to 
zero, then Z"(0) + %. On the jet boundary (7 =1), the solution of the equation of motion (4.16) 
will satisfy the third boundary condition (4.17a) provided that 


Z’ (1) =0. (4.18a) 
Let us transform the equation of conservation of excess enthalpy (4.5) with the help of the 


functions (3.3) and the variable g~. We obtain 


Hh >, 
Pm&m (Lm — bn) AX 2 ad a ro(fon — 
tig ip = ln) ( R, | i, i, I (9) edo =0.5. (4.19) 
Let us introduce the abbreviations: Ai, =i,—iy, Ai=i-iy, Aip=ip—iy, and e =Ai/Aim. 
We can find Ai,, from Eq. (4.19), 


Ai, =~ Vi-mt, (4.20) 
M Xx 
where 


K _ 0. Dag lly Aly 
Pb 
a \2 
Ro Cm \ e(p) f (9) ¢de 
0 


0 


Let us write the energy equation (4.3) as follows: 


oa la) FO) gaye | 2 yp ee |. (4.21) 
Ox oy y oy oy Oy 
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~] 


vl 


Let us transform the energy equation (1.21) with the hclp of expressions (3.3), (4.10), (4.11), 
and (4.20) and it will become 


Fe ny — m)dZ's', (4.22) 
5 
where n=1/Pry7. 


Dividing the energy equation (4.22) by the cquation of motion (-1.13), we obtain the energy 
equation 


Z g’ a 
Pr, 7 —_——— . (-4 2 Oo) 


For solving Eq. (4.23) we will use boundary conditions analogous to boundary conditions (3.17): 
T. On the jet axis m=0 (n = 0); Z(0)=1; ¢(0) =1. 
2. On the jet boundary g= Py (7 = 1); Z( Oy) =0; (gp) =9. 


The solution of the energy equation (4.23) together with the boundary conditions is of the 
form 


Pr 
ea *. (4.24) 
Equation (4.24) allows us to calculate the enthalpy profile across the jet. 


Solution of the Equation of Motion. Wewill solve the cquation of motion 
(4.16) with the main boundary conditions (4.17), (4.17a), and (4.17b) and the supplementary condi- 
tions on the jet boundary. 


In Eq. (4.16) 66 is an unknown function and 7 a quantity which will be determined from a 
simultaneous solution of the equations of motion and cnergy, togethcr with the equation of state (1.1). 
We will seek the general solution of the equation of motion by Galerkin's method [380]. Lect us 
write down the solution in its general form as 


Z,= S aye, ({=0, 1, 2, 3..... hk). (4.25) 


The system of functions 0: Will be chosen to satisfy the boundary conditions (4.17), (4.17a), and 
(4.17b), as well as the supplementary boundary conditions (4.18) and (4.18a). Let us consider 
the following system of functions: 


— l, 2, — 37° — Ay), 9% —= (] a mh 5 

3 (k-1) 

ys ” 2 

Gy (1 —4)?y", 22. p= — Ya (4.26) 


If we set @=1 and @,=1, then the sum of functions (4.26) will satisfy all boundary condi- 
tions for any values of the coefficients a, a, ...,@. We will seek the solution of the cquation 
of motion (4.16) with the help of the first threc functions of system (4.26). The solution will be 
of the form 


Z»= 1+ (877 — 44°") aa (1 — 2x), (4.27) 


In this case, the solution of Eq. (4.16) reduces to the determination of the two unknowns 
a, and bx. We can write down two equations according to the Galerkin method for the determi- 
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nation of these unknowns, 


J L ~~ 
1. f L (Z.)9,d7,=0, 2. | L (Z») ¢.4%4=0. (4.28) 
( 0 


0 


Transforming Eqs. (4.28), we obtain two equations for the determination of the unknowns A and 
bx; 


1. Aj+ A3;—6°Bi — 6, B3=0, | 


x * 3 * * (4.29) 
2. AS+Ai+0°B:— 6,B;=0. 
In the first equation of system (4.29), the functions A* and B* are given by 
Aj = 1.18697 +0.10579a, +0.103C6a3 + 0.004425, 
As= Ay +@.A;+a34,+43Ay, 
* in m ~ 
=/0. 0.01275 +. 0.01884 
By (0 03591 + 0.0613. “)+( 01275 +.0.0188 ae + 14.30) 
+ (0.00219 +-0,00204 a3 -+ 0,00018a3, 
—m 
B3=an(B,+a.B, +05B,+38,), 
where 
A, =324J,—1512/u:, + 2628/, — 2016/4, + 684/, — 360J, + 396/, — 144J4,, 
A;=486Su;,— 1746, + 1320S, + 1884, —2796Ju;, + 1104S), — 12/, — 186 Jy), — 126/, + 72S», 
A, = = Jy — 588J1), ~ 383 4+ 1848, — = J,— 2161S), + 
~ 7 
4760S ¢-} 1172S, — 2 J = 322Iui,-+ 6344 + 60S, + Js — We, 
A.= — Ju, — 499 — = Joop, +238, + a Jo, —471Jg— 
— 363 Jin, +484), A Jn = 27g Sis + 8g -+ Sn 
B,=162/, — 810S13,,-+ 1512S, —_— 1248 Ju, + 492 J, — 396F +- 480/, — 192/;,-+ 18/5 — 42/s,+ 24/., 
(4.31) 


B,= ve J, —774J 4 + 633 Jy, +916, — J, — 106 Sg + 
. 694Jw), + 120/, — Se 170, + 81S: ,-+ 20/4-+ oJ ~6)., 
Bg=81S yp —210S wy, — 152+. 744J 2), —8 1g — 1000/5, + 
+-336S_ + 648Jis), — 217 Sg — 234 uy, + 245+ 64s, +9, —12Jn, 


By = Jey, = 14S) 87), 16S 99-+  Je— 1T0Sy — 


—150Su ,-+ 200+ Je, — 130/,—33 J, + 44/4-+ - Ju, —6J 5 
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In the second equation of system (4.29), the functions A* and B* are given by 
Ax = — 0.25095 +- 0.008362. — 0.01454a3 — 0.00005a2: 


Ay— A, +a, A,+ 03d, + agg: 
Bi =( 0.00772-+-0.01198 = ) + (0.00237 + 0.00303 m \x 


1— m 1. m (4.32) 
X Gob (0.00035 +0.00027 ) @3-+.0.00003a3; 
— In 
By=an(B,+a,8,+438,+a38,), 
where 
A, =108/s,— 360; + 180s, +576, —648Ju,—144/, -+360Js, —36/>, — 72, + 36J;,, 
Ag= 162Jy— 366 ir, — 372/ + 1296S, -+ 78, — 173 4d, 
4360/5 1104Ju),-— 258, — 378», + 12, + 96J::, + 18/4 — 18s: 

Ag => Soy, — 91S 9 = > Jn, $466 $A Jin, = 981g Jy, + 1084, 

+711 So,—661J_— = Ju, + 210F, 4 ne Ju, —21F ,— = J.,+ = Je, 
Ag == Sy = 84S 1 +O Sy — 416 Sy + Sg 2976S, + 1094, — 24, + Is 
B, = O4 S17), — ] 98), +- 132Siy, + 300/, — 390S sp, — 90/ + (4 .33) 


+-216Ju), + 72S, — 54Jy, — 90, + 36Js),-+ 12S, -+ 6s), — 6/9" 
By = Sy) 168m, — 148, +624), — 904 Sy, + 192J, 
+ 672Jis;,— 105, —31 2S), — 365 + 112Je;,-+ 36, — 24S), — 8, + =- Jy 
By = 27J xy, — 34) — 150S 2, + 188) go -+ 352), +430) ) — 458), + 520g +370 y, — 350), — 202J m4 
+124, -4+80Ju, — 18S; —22Je, +3Jy, 
Ba= - S13 — 26S 42 + 84S — 154) yp 175g — 126 y + 56/5 — 14 ¢-+ ~ Jy. 
The integral J; appearing in these functions is given by 


1 
Jy=\ tay (b=1, Yay 2, Yay 3...) (4.34) 
0 


The first equation of system (4.29) can be written as [33] 


b°4+-3P,b, +2Q,=0, (4.35) 
where 
B* A* 4 A" 
pat SB, gat. 
B; 2B; 


The solution of this equation yields the value of b,. Since according to its physical mean- 
ing b, Should be a real positive number, we can immediately reject irrational roots to obtain 


b, =e + Ue, (4.35a) 
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where 


mV. Q+] GLP, Us] = V- —) i+ Pi. 


The second equation of system (4.29) is solved in the same manner. This equation is 


b?4+-3P.b, +2Q.—0, (4.36) 
where 
Py == — Bs » Qye= Ane A 
3 Bs 2B, 
The real root of this equation is 
b= le -+ Ung, (4.36a) 


where 


m=V —-Q4+V G43, m=V —-@—|] +P. 


The system of equations (4.35) and (4.36) is solved graphically, i.e., we determine the point of 
intersection of the two curves. 


The Calculation of an Axially Symmetric Jet of Real Gas. Wetake 
several values of the density py, (or the specific volume V,,) on the jet axis, assuming that the 
density is constant over an infinitesimally small region of the axis. 


We will perform the calculations for each value of the density py, =const. 


In the first approximation, we assume that the variation of the specific volume across a 
section of the jet is as follows: 


Via \n 
p= (7) =o (4.37) 


In this case, the known function 9'/g@ is of the form 


, 


6 
— = In 6,=const. (4.374) 


Then, the functions A¥ , B¥, Af , and Bf of system (4.29) will become 


A3= In 0, (— 0.21797 — 0.01916a — 0.00756a3 — 0.000923); 
B3 =an In 0, (0.0833 + 0.01664a, +-0.00204a2 + 0.00012a3); 
Aq=In 6, (0.03954 — 0.00114a,-+ 0.00065a5 + 0.000123); 


Bi, =an In 0, (— 0.01663 —0.00206a, — 0.000143). 


Taking several values of @,, we now solve the system of two equations (4.35) and (4.36) and mak- 
ing use of expressions (4.30) and (4.32) for the functions Af, Bt , A¥, B} and expressions 

(4.38) for A* , BY, Af , BY, we determine b+ . From Eqs. (4.27) and (4.24) we find the veloc- 
ity and enthalpy profiles in the boundary layer of the jet. Knowing the function ¢(n), we can 


(4.38) 
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easily find the enthalpy iat any point of the cross section of the jet boundary layer 


mer (4.39) 


In Eq. (4.39), the enthalpy of the parallel flow is known, while the value of iy, can be easily 
found from the known value of py,(Vyp)) on the jet axis by means of the graph of i=f(V). Knowing 
the enthalpy across the jet, we can determine the values of the specific volume from the graph 
of i=f(V) and thus we can find 9 = V/V. 


In other words, we find the new improved specific-volume profile over the cross section 
of the jet. 


Equation (1.1) can be written as 


= (ig — by) + ig = AV 6+ B = 
v2 02 
Differentiating this equation, we find an expression for 9'/@ 
Oe im i) Hn fw —~ (4.40) 
AV 7) — v0 V2 03 


Using Eqs. (4.40) and (3.34), we can calculate the function 9'/@ and hence evaluate the 
integral (4.34) by the trapezoidal method [32]. Todo this, we subdivide the range of n(0—1) in 
the graph of 9'/6 =f(n) into r subintervals each of length 1/r. The trapezoidal rule for the 
evaluation of the integral J;, is 


I Q’ ’ Q’ 0)" 
n= [(T) a t(G) et HF) $05 (FE) (4.40 
Having obtained the values of the integrals J}, we calculate the functions Aj, A», ..., By, 
Bo, ..- from (4.31) and (4.33) and then find the values of Af, A}, ..., Bf, BY from Eqs. (4.30) 


and (4.32). 
Solving the system of equations (4.29), we find new values of @, and bx. 


Having calculated several cross sections in this manner, We can proceed to the determina- 
tion of the variations of the jet parameters along the axis. 


With the help of expression (4.7), we can write down the equation of conservation of excess 


momentum (4.4) as follows: 
— My PO 
ln 0. 0.70/ ae V — VY -, 


“uo 
ob, F'Zdy 
Ry ; 


where mp =Uyj/Up, and py are parameters at exit from the nozzle. 


(4.42) 


In this equation, F' is given by 


Fl = [m+(1—m)Z] 4—an = Z'(1—m) —> 4, (4.43) 


« 
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We can find the value of ax/R, from the equation of conservation of excess enthalpy (4.5) 


; l 
| — , a f  % 
0.707 / a \/ \ Fiza | —_ 
1 — mm, , vm 


a (4.44) 
Ry 


Atm 
— 6, \ Fredy, 
Aly e 
) 


We can now uSe Eqs. (4.48) and (4.44) to determine the variation of velocity along the jet axis 
and the coordinate of the given cross section. If the jet is emitted into a parallel stream, then 
the calculation is performed by means of a trial-and-error method involving a variation of m. 


In order to match the initial and main regions of the jet, we calculate the transitional sec- 
tion of the jet where py=fPo, and Ai =Aiy. The value of uy /uy in this transitional section is 
larger than unity or, in other words, the thermal transitional section is situated closer to the 
jet nozzle than the dynamic transitional section. From the intersection of the curve Up/Uy = 
fiax/R)) with the line uy)/Up=1, we find the coordinate of the dynamic transitional section. 


In this way We can obtain all of the parameters of the jet and their variations along the 
axis and cross section of the jet and, finally, we can match the initial and main regions of the 
jet. 


Sample Calculations 


1. Incompressible Case. We determine the velocity profile in the main region of an axial- 
ly symmetric jet of an incompressible fluid. 


Let us solve the system of equations (4.29) taking into account the fact that (V) =const. 
and uzyz=0. We find the values of the coefficient a, and the dimensionless width of the boundary 
layer b, to be 


ay9=—0,34; b, =3,32. 
Substituting the value of the coefficient a, into Eq. (4.27), we determine the velocity pro- 
file across the cross section of the jet. 
The results of the calculation are given in Table 13. 


Let us compare this approximate solution for the velocity profile of an axially symmetric 
jet of incompressible fluid with the more accurate Prandtl—Tolmin solution. The relative width 
of the boundary layer in the latter solution is b, =3.4 [17]. The velocity profile u/uy, = FY/7 
obtained from the Prandtl—Tolmin solution is given in Table 13. 


Let us determine the variation of velocity along the axis of a jet with a circular cross 
section. To do this, we first of all evaluate the integral 


] 
b* \ Z24,d1,=0.506. 
Q 


Substituting the value of the integral into Eq. (4.42), we obtain the formula for the axial 
velocity in the main region of a circular jet of incompressible fluid 
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TABLE 13 
: 7 
r (0 p. 0,2 0.3 lod 08 0.6 07 0.8 jO.9 11.0 
Z | Q.8953)0.7436.0, 5846 0.437 (0.305 |0.1978'0.1101)0.0473 0.0171 0) 
Ul | | 
—_ =- ~~ 1 (0.90-E 0.259 0-59 O.44 10.302 )0.187 [0.1 0-4-4 0.017 | () 
mt 
| | | | 
TABLE 14 
| Vin | Ore | im Dm | a 6, 
| | 
Section | 9.5-1074 | 18.8 CO 2-48 1.55 | 0.2] 
Section II 12.5-10-4 14.28 120 998 1.1 | 9 3 
Section III 20.5-1074 8.71 140 318 1.0 9.45 
Section IV 40-10 4 4.46 160 334 0.6 2.7 
2.87 


Section V 60-104 2.98 180 388 0.3 


The formula for the axial velocity obtained by Tolmin [17] is of the form 


Um __ 0.96 
Uo ax 
Ro 


A comparison of the approximate solution for the case of an incompressible fluid with 
the more accurate Prandtl—Tolmin solution shows that the agreement between the two solutions 
is satisfactory. 


2. The Main Region of An Axially-Symmetric Jet of a Real Gas. Let us calculate a jet 


of carbon dioxide at p=981 -10' N/m’ and an initial temperature T,)=248°K. The parameters of 
the surrounding medium are: Ty=908°K, Vy=1.82 -107? m°/N, ig =137-104 J/kg. We assume 
that Pry =0.5 and a=0.1. 


We take several values of the specific volume on the jet axis, V;,. Using these values of 
Vm, we find from Figs. 2 and 3 the values of i,, and Tm. As a first approximation, we take the 


Z 
"ST LTT 


\ Wii Section I a 
RY [section im [| | 
Nidxq Section V Pf | 
PET MAL TE ET 
CORR ECE 
Incompres- ONG 
sible fluid NN 
See eeeENe 
CCCs 


0.5 


~ 


0 10n 107 
Fig. 24. Velocity profiles in Fig. 25. Enthalpy profiles in 
the main region of an axially the main region of an axially 


symmetric jet. symmetric jet. 
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TABLE 15 
, 0 0.1. 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Z 1.0 | 0.944 | 0.852 | 0.737. | 0.609 | 0.472. | 0.338 |} 0.21 0.101 | 0.0333 | 0 
_ : 1.0 | 0.972 | 0.923 | 0.859 | 0.78 | 0.687 | 0.581 | 0.458 | 0.318 | 0.183 | 0 
: 8 1.0 | 1.05 116 | 1.47 | 2.26 | a.7t | 7.168 | 9.64 | 1205 15.25 | 18.8 
5) 2 1 | 0,868 | 0.665 | 0.444 | 0.313 10.0775 | 0.0248 | 0.0075 | 0.0015 | 0.0001! 0 
Om te, | | | 
— .. 
a “on 0.721 | 0.518 | 0.401 | 0,113 | 0.0428 | 0.0164 | 0.0049] 0.0005 0 
1 OmUm* Alm | 
Z 0.93 | 0.82 | 0.693 | 0.559 | 0.423 | 0.297 | 0.181 | 0.086 | 0.029 | 0 
: 0.065 | 0.906 | 0.834 | 0.748 | 0.651 | 0.546 | 0.426 | 0.252 | 0.169 | 0 
= 0 1.22 | 1.83 | 2.54 | 3.22 13.95 9 4.71 | 5.61 6.64 | 7.56 | 8.71 
c Ou 
2 : 0.814 | 0.434 | 0.222 | 0.1183 | 0.0547 | 0.0232 | 0.0077 | 0.0016 | 0.0002 | 0 
® Om: ny 
yu ME | | 
a 0.815 | 0.479 | 0.208 0.1586 , 0.0842 0.0427 0.0181 | 0.0054 | 0.0011, 0 
OmUm* Alm | | 
| 
Z 0.9117 | 0.78 | 0.636 | 0.495 | 0 362 | 0.045 | o.144 | 0.056 | 0.023 | © 
: 1 0.955 | 0.884 | 0.797 | 0.704 | 0,602 | 0.496 | 0.379 | 0.256 | o.151 | 0 
. 0 1 1.088 | 1.24 | 1.415 | 1.6 1.8 2.03 | 2.26 2.5 2.7 2.98 
6 OU-u . 
2 s ! 0.796 | 0.52 | 0.307 | 0.166 | 0.0796 | 0,033 | 0.0105 | 0.0021 | 0.0002 | 0 
$ Om Um 
1p) 
ou Mi | | | | 
i 0.834 | 0.588 | 0.381 | 0.236 0.133 0.067 | 0.0278 ' 0.0082! 0.0016! 0 
Om ms Aim | ! | | 


function 9'/@ to be given by Eq. (4.37). 
the coefficients A* and B* corresponding to these @, from Eqs. (4.30), (4.32), and (4.38). 


Taking several values of a, , we determine the values of 
From 


the point of intersection of the two curves b+ =f@), we find the solution of the system of equa- 
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Fig. 26. Specific-volume 
profiles in the main re- 
gion of an axially sym- 
metric jet. 


tions (4.35) and (4.36) (see Table 14). From Eqs. (4.27) and 
(4.24), we find the velocity and enthalpy profiles in Figs. 24 and 
25. The results of these calculations are given in Table 15. 


Using Eq. (4.39), we find the variation of enthalpy across 
the jet and from the graph of Fig. 2 we determine the value of 
the specific volume vy and the quantity 9= V/Vy (see Fig. 26). 


Using the equation of state (1.1), we find the improved val- 
ue of the function 9'/6 from Eq. (4.40). From Eqs. (4.27), (4.43), 
and (4.24) we determine the total-head profiles and the excess- 
enthalpy profiles (see Fig. 27), 


aS (), 
Pmt 
pu-di 

=] (7, 
Ta ai, 


The results of these calculations are given in Table 15. The 
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TABLE 16 
| Section I | Section | Section H/seetion III | Section 1V|_ Se Section II|/Section If Section IV| Section V 
ty) | 
— l 1.32 2.16 4,2] 6,3] 
Om 
Ato 
rr 0.874 0.788 0.704 0.62 
\F’ Zdy 0.0538 0.03423 | 0.03302 | 0.03705 0.0411 
| 
JZ’ cdr 0.0647 | 0.04161 0.041 ().04783 0.0545 
6 ! | 
an | | 
— 1.145 | 1.8 2.38 3.07 3.73 
Ry 
sm 1.205 1.055 0.945 0.900 ).818 
iy 
calculation of the integrals 
I 1 
puru ( F’Zdy, and \ F’ed7, 
oe é 
) Cae is carried out graphically. 
The dimensionless quantities ax/Ry and 
Section I 0 
aN ection aes Um/Uy, calculated from Eqs. (4.44) and (4.42), are 
an (Ree shown in Table 16. The results of the calculations 
as WN sotion aan were uSed to construct graphs showing the vari- 
Cc . . . 
TT AYN TTT ation of the relative velocity Um/ Uy and the rela- 
rT TAWA YT PTT tive enthalpy Aim/Aiy along the jet axis (see 
. Fig. 28). 
Ree | 
Pt dt PN] The same figure also shows for comparison 
0 Qs 10n the variation of the relative velocity um/u, for a 
Fie. 27. Total-head files ; jet of incompressible fluid. For each section we 
1B: . ° a ca prow es in determine the relative width of the jet boundary 
the main region of an axially layer 
Symmetric jet. 
Pa ax 
— b, 
Ro 
Um 
i0 ~ The variation of the relative width b along the 
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Fig. 28. Variation of velocity and en- 
thalpy along the axis of an axially sym- 


m 


etric jet. 


jet is shown in Fig. 29. The relative width of a 
jet of incompressible fluid is alSo shown in the 
same figure. 


The results of the calculations show that 
as the ratio p,,/py increases, the velocity (Fig. 
24) and enthalpy (Fig. 25) profiles become more 
convex. In the case of the carbon dioxide jet, 
the velocity falls more slowly along the axis 
(Fig. 28) and the boundary is less steep (Fig. 29) 
than in the case of a jet of incompressible fluid. 
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In the calculation of a jet of carbon dioxide, we 
have replaced the true variation of the function 9'/9 ap- 
pearing in the equation of motion by the constant 


{y’ 
op ine 


Let us carry out a more accurate calculation (second 
approximation) for the third cross section (Aiy,/Aiy = 
0.788) of the jet by introducing the function 9'/9@ obtained 
in the first approximation. Todo this, we use Eq. (4.34) 
Fig. 29. The variation of the to evaluate the integrals J, and Eqs. (4.31) and (4.33) to 
width of an axially symmetric jet. find the values of the coefficients A and B. 


Solving the system of equations (4.29), we find new values of the quantities a, and by, 
namely, 


Using the same equations as in the first approximation, we determine uy)/Up and ax/R, 
(Up /Up = 9.957, ax/R,y =2.24). A comparison of these values with the results obtained in the 
first approximation shows that the first approximation is satisfactory, while it is much simpler 
than the second approximation. 


Therefore, for future calculations of jets of real gases, it is recommended that the first 
approximation be used. 
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THE CALCULATION OF THE SHAPE 
OF AN ISOBARIC MIXING CHAMBER 


O. V. Yakovlevskii 


One of the main problems facing designers of ejectors, as well as various mixing cham- 
bers, is the determination of the optimal length and shape of the channel in which turbulent inter- 
mixing takes place. Some authors have uSed for this purpose the laws of propagation of free 
turbulent jets in a parallel stream. 


Two papers by A. Ya. Cherkez [1] and M. V. Polikovskii [2] have been devoted to the cal- 
culation of flow in a mixing chamber. However, as was pointed out in [3], the method of Cherkez 
based on the experimentally observed analogy between the velocity fields in a cross section of 
the gas stream in an ejector and in a free turbulent jet is only of limited validity. The limita~ 
tion arises becauSe the nominal velocity of the parallel stream introduced into the calculation 
is found to be imaginary under some conditions. 


The method proposed in [2] for the calculation of the inlet (initial) region of the ejector 
involves the determination of Successive approximations to the stream parameters for a num- 
ber of transverse sections. This method leads to laborious calculations and does not allow us 
to make an analysis of the influence of hydrodynamic and geometric parameters on the proper- 
ties of flow in the mixing chamber. 


The attempt made by I. S. Bogolyubov [4] to determine the shape of the isobaric inlet re- 
gion of an ejector has been based on the theory of intermixing of parallel streams developed by 
A. Kuethe. However, A. Kuethe's theory contradicts experimental data on free jets (for example, 
see [1]), so that the uSe of its results for the calculation of flow in a mixing chamber has no 
practical meaning. 


Moreover, the papers [1, 2, 4] listed above do not allow us to calculate the turbulent jet 
flow of fluid in a channel when the active stream, i.e., the stream possessing the higher mo- 
mentum is the peripheral (annular) stream and not the central jet. 


An analysis of the ejection properties of turbulent gas jetsled the author in 1960 to formu- 
late the hypothesis of the universality of these properties [5]. Calculations carried out on the 
basis of this hypothesis were found to be in satisfactory agreement with the corresponding 
experimental data for turbulent jets propagating under various conditions (jet of hot gas, includ- 
ing plasma, in a parallel stream; supersonic jet under various conditions of outflow; gas-liquid 
double-phase jet). 


It was shown in [5] that the method based on the universality of the ejection properties of 
gas jets can be extended to the solution of other jet problems. 
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In the present paper, this method is applied to the solution of the problem of determining 
the shape of a mixing chamber which results in the intermixing of fluid streams under isobaric 
conditions. As is known, other conditions being equal, an isobaric mixing chamber is always 
more effective than a cylindrical chamber, i.e., it secures lower total-pressure loss for the 
mixture. 


It should be noted, that for the calculation of flow in an isobaric mixing chamber and the 
determinination of its shape we could uSe the well-known theory of turbulent jets. However, as 
has already been noted in [5], the theory of jets propagating in a parallel stream [1] does not 
agree with experimental data over a certain range of the principal parameters. It was for this 
reason that a new approach to the calculation of gas jets propagating in parallel streams became 
necessary. On the other hand, the application of the laws of turbulent-jet theory [1] to isobaric 
mixing chambers does not allow us to obtain the required relations in explicit form, which con- 
siderably complicates the analysis of the dependence of the chamber shape on the principal flow 
parameters. 


Therefore, the method of calculating turbulent jets based on the hypothesis of the univer- 
sality of ejection properties of jets [5] has been used for the solution of the above problem. 


1. Description of Flow, Principal Assumptions, 
and the System of Equations 


Let us consider flow in an axially symmetric channel (see Fig. 1) in the initial region of 
which there are two streams with unequal velocities U;) (central jet) and Us) (annular jet) sepa- 
rated by the wall of the central nozzle. Because of vortices that are present on the unstable 
boundary separating the two streams, an intense turbulent intermixing of the central and annular 
jets takes place. The region in which mixing takes place widens with distance from the initial 
section and penetrates more and more deeply into each of the two streams. As a result of this, 
in some transverse Section (which we will call the transition section), one of the streams be- 
comes completely involved in intense turbulent mixing and the unperturbed core of this stream 
disappears. If the central nozzle is of small dimensions (by comparison with the transverse 
section of the whole channel), the unperturbed core of the central jet ceases to exist first and 
the inner boundary of the mixing zone intersects the axis of the channel. On the other hand, if 
the transverse dimensions of the central and annular jets are comparable in the initial region 
(see Fig. 2), the unperturbed core of the annular jet is the first to disappear and the outer boun- 
dary of the mixing Zone touches the channel wall. The flow depicted in Fig. 1 will be called for 
brevity Scheme A, that depicted in Fig. 2 will be called Scheme B. 


Downstream from the transition section, 


} ° F if further broadening of the mixing zone leads to the 
} | disappearance of the core of the second stream 
| Pe also (the annular jet in the case of Fig. 1 or the 
| U, central jet in the case of Fig. 2) and this occurs 
Ko j at a certain distance from the nozzle at the posi- 
| ; tion of the section labelled V. In section V the 
R, Ry y, entire flow is found to be completely subject to 
Primary one! Secondary | x strong turbulent intermixing. Downstream from 
Zone section V— in the main region of flow -— a further 
Initial region |_Main region decrease in the nonuniformity of the velocity pro- 
file takes place due to the action of the process of 
Fig. 1. Flow associated with the propa- turbulent diffusion. Flow in the main region has 
gation of a jet in a restricted parallel been studied in detail in the author's papers (3, 6] 


stream according to Scheme A. in which formulas for the calculation of the stream 
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parameters for an arbitrarily shaped channel 
were derived. 


In the present paper we will consider the 
calculation of flow in the initial region of the mix- 
ing chamber between the cross Sections labelled 
Oand V. In the case of an ejector, this region is 
: called the inlet region; it is in this region that the 
Primary 1 Secondary | ejected medium becomes entrained by the active 

Zone zone stream in a process which is completed by sec- 


| 
| 
| 
| 
| 
| 
| 
| 


: Initial region | Main region tion b 
Fig. 2. Flow associated with the propa- The problem can thus be formulated as 
gation of a jet in a restricted parallel follows: To find the shape of the initial region 
stream according to Scheme B. of the channel which with the given geometric and 


hydrodynamic parameters of the flow in the initial 
section would result in a constant static pressure along the stream (p=const.). 


We will assume that the velocity distribution of the two streams are uniform in the initial 
region and that the magnitudes of the velocities in the central and annular regions are U4) and 
Uso, respectively. The densities of both streams will be assumed to be constant (po =const.). 
With pressure constant, the velocity in the unperturbed core of each stream will not vary. By 
analogy with a free turbulent stream, we will assume that the profile of the dimensionless ex- 
cess velocity in the mixing zone is Self-similar. In the system of coordinates illustrated in 
Figs. 1 and 2, this assumption can be expresSed as 


_ R—-R; 


’ U—wU 
AU%'—_—+_ =o (8), == , 
U,—U, ? G) Ry—R, 


(1.1) 


Here U is the velocity of the stream at an arbitrary point on the X axis; the indices 0, 1, 
and 2 correspond to the initial section, the central jet, and the annular stream; R is a transverse 
coordinate measured from the axis of the channel; R,; and R» are the inner and outer boundaries 
of the mixing zone; Rc is the channel radius. 


It should be noted that in the primary zone (between Section Oand V), the velocities U, and 
U, are constant and equal to Uj) and Us), respectively; in the secondary zone, either the velocity 
Us, = Usg (in the case of flow according to scheme A) or the velocity U, =U, (in the case of flow 
according to scheme B) remain constant. 


Thus, to solve the problem formulated above, we have to determine the X-dependence of 
the channel radius Re, as well as of two other quantities: either R, and Ry of the primary zone, 
or U,; and Ry, or U, and R, for the secondary zone, or U,; and Uy, for the initial region. In order 
to determine these three unknowns, We have three equations at our disposal, namely, the equa- 
tions of continuity and conservation of momentum, as well as the ejection equation based on the 
hypothesis of the universality of the ejection properties of turbulent jets [5]. 


The flow rate and the momentum of the liquid in the initial region at constant specific 
gsravity are given by 
Gy = ROU 1 aa (Reo —_ Ro) TO x 
and 


Ky=aRo jot = (Reo — Ro) 120. 
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Por converfence tn subsequent (ransformattans, wo will Introduce the areca-averaged 
veto tty go in the initial region of the channel 


Ciy 
Lea) im ’ 
Ve 
aa NY 


os wellas (he flow-rate averayed velocity ky In the initial seetion of the channel 


K\, 

ky N 
ry at 

In the following we will also use the quantity i) Ky/pfo", the so-called field coefficient [3], 

whieh characterizes the kinematic nomimiformity of the flow in the initial region of the channel. 


As apreometric parameter of the problem characterizing the How conditions In the initial 
reyion, We choose fhe dimensionless quaontity a defined as the vatia 


Ry ) 
(i. 
itis not ditientt to see that all of the dimensionless flow parameters in section o will be 
completely determined by two dimenslouless quantities which we take to be the parameter @ 


and the relative velocity of the annular streany in the Initial section, wag Us9/pg. It ts ensy to 
show that the other characteristics of the Tow can be expressed In terms of Usg and @was follows: 


| Uy ; t rm 
Uy day | y, ty 1 (lo da), 
( it 
t a) t yu 
Ny "3 . (1.2) 
iv | Hoy (| t) 
N Chay a w(t 1) 
Chyy { Moy (| a) 


The equatton for the conservatton of flow rate in the Initial region of an axially syminctric 
Channel can be written as 


Ie 
Ryo RU ALLY CRIR | (RE RYU, (1.3) 


‘ 
sy 


While the equation of conservation oF momentun will be 


he , 
hyo RIEL EM RAR | (Re RDU. (1.4) 
hy 


HW in the (ntegrands ot Mqs. (1.3) and (1.t) we replace U and R by AU® and £ with the help 
of formula (1.1), then afler a series of transformations these equations become 


nei at LW ane a | (W75 AM, (i 1) 


nnd 


tory, MOaryradi ears Mas (1.6) 
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Here, we have introduced the following abbreviations 


l—riu i, ~ Que to reu5 

2 yo yt 9 

My== 5, Myse 
Uy —- Me (iy > Un) 
ial R., Re 

Pym yO PQ 
Rev Rew Rew 


aq,=1+A,—2A b, -A,,--A,, C=A_, 


b,=B, -BL, o=B,, 


+ 


wo) 


If, for the description of the velocity profile in the mixing zonc wc makcusc of Schlichting's 
formula [1] 
eQ=(l-E Ge?) 


which is used extensively in the theory of turbulent jets, then the cocfficicnts listed above will 
have the following numerical values: 


A,.= 0.450, B,.=0.316, 
a, =0.357, b, =0.193, c, =A, = 0,257, (1 .8) 
a, = 0.502, b, =0,182, cy B,=0.134. 


AS the third equation required for the determination of the thrce unknowns, wc uSc the 
ejection equation. According to the hypothesis of the universality of the cjcction propcrtics of 
jets [5] mentioned above, the rate of increase of the relative mass of the jet is indcpcndcnt of 
the relative velocity of the surrounding medium,! but is completely determined by thc ratio of 
the densities of the surrounding medium and the jet. Since we are considcring the mixing of 
incompresSssible-fluid jets at constant density, then, conSequently, the cjcction cquation will hc of 
the same form as that for a submerged jet. An analysis of the cxpcrimecntal and thcoretical data 
obtained by various authors [1, 5, 8] shows that the variation in the rclative flow rate with in- 
creasing distance from the initial section of the jet can be describcd by the following function: 


lt+px for x<X,, 


0 a 
Ge (x)= w(x+6) for XZ X,. 


(1.9) 
Here, Ge =Ge/Gj9, X=X/Ro, while p, w, and 6 are constant coefficicnts. Thc cocfficicnt 

w according to {1, 5, 8, 9] is a constant and equal to 0.155. The magnitudc of the cocfficicnt p, 

as Well as the values of 6 and xy will, in general, depend on the initial turbulencc of thc jet. 

However, for the case of turbulent channel flow most frequently cncountered in practice (turbu- 

lence intensity of 1-5%), experiments [1, 5] show that the coefficient p can be taken as 0.072. 

The remaining constants in Eq. (1.9) then are x, =12.4 and 670. 


In order to make use of the ejection equation when uy) < 1, we will writc the cquation of 
conservation of flow rate in the channel as 


Gy=G.+%(Re — R2) 1 29. 
T This assertion is valid as long as the jet is active, i.e., as long as its momcntum pcr unit 


cross sectional area is high than that of the parallel strcam. If the parallcl strcam beccomcs 
active, then all arguments refer to it. 
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Putting this equation into dimensionless form, we obtain 


2 2 
0 L = ty (rg — 75) 


Ady 


which taken together with the ejection equation (1.9) allows us to determine the variation of all 
variables of interest to us along the channel 


] -- a Go (x) 


NS 
») 


| 
. 


ow 


ay (x). (1.10) 


Un) 


In the case where the annular and not the central jet is active in the mixing process 
(Uy) >1), the ejection equation with the help of the hypothesis of the universality of ejection 
properties [5] can be written as 


Cy — G, = w/LU.y2cR,X, 


where G, is the flow rate of the fluid in the unperturbed core of the central jet and y' is a coef- 
ficient which, according to the data of [10], is equal to 0.035. Since G, =m RyyUj, this equation 
transformed to dimensionless form becomes 


9 Uy — Qu UaXr 
PE it (1.11) 
uy 


For the analysis of the primary zone of the initial region where the unperturbed cores of both 
jets still exist, we obviously have U,= Uj) and U, =Us9, so that formula (1.11) can be simplified to 


ri=2 (1 — 2u! log x| — by (x). (1 12) 


2. Determination of the Boundaries of the Flow Regimes 
In the (u»9,, a) Diagram and the Solution of the System 
of Equations 


Before we proceed to the solution of the problem, i.e., the calculation of the function re(x) 
for various values of the parameters Us, and a, let us attempt to establish the conditions for the 
existence of the flow regimes A and B as functions of the parameters Uy) and a. In this con- 
nection, aS a criterion characterizing the transition from one flow scheme to the other we will 
use the simultaneity of the disappearance (in the one transverse Section) of the unperturbed 
cores in the central and annular streams, namely, 


rj=0, r=h. (2.1) 


Turning to Eqs. (1.5)-(1.6), we see that their 
simultaneous solution together with conditions (2.1) 
yields the connection between Uy, and a which must 
be satisfied by the geometric and kinematic charac- 
04 06 08 0 12 846 16 18 Uy teristics when conditions (2.1) hold: 


Fig. 3. Regions of existence of flow 


, {one naan enw 2 _ 
according to schemes A and B as a [/ 0,508 — 0.9942ta) + 0, 552u59 — 1 + 1.23712 (2.2) 
. CSS oo eeoooooowrrrre_e_e__c cc ln ° 
function of the parameters Us, and a. 2th 
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The curve representing dependence (2.2) is shown in Fig. 3 (curve 1); it subdivides the 
whole domain of the parameters uy, and wv into two regions: A and B, the first of which corre- 
sponds to channel Now according to scheme A, the scecond to flow according to scheme B. In the 
shaded portion of the usy, w plane, the values of the paramcter uy are negative, so that this re- 
gion of parameters uy, and v must be excluded from consideration (curve 2 corresponds to uy, 7 0). 


Jtaving determined the boundary between the two low regimes, we will proeced to the so- 
lution of the system of equations describing the flow of intermixing jets in the channel. 


The flow of Muljd in the main region of an isobaric mixing chamber (sce Figs. 1, 2), as 
well as the shape of the channel can be determined by the method described by the author in [6]. 
As shown by calculations, the varsation In the form of the Jsobaric mixing chamber in this re- 
gion js small by comparsson with that in the initial region. In the present paper we give the so- 
lution for the Initial region of an isobaric chamber, while for convenicnee and for greater defi- 
nitencss of the resulting formulas, we wil) consider separately the primary and secondary 
zones of the Initial region. 


The Primary Zone of the Jnitlal Region. This region is characterized by 
the existence of cores of unperturbed flow in both the central and peripheral jets. In view of the 
condition p const., the velocity In cach core Is constant and equal to the corresponding initial 
value (Uy Uyg, Uy Uy). The unknowns to be determined are ry, ro, and re. 


The simultancous solution of Mys. (1.5) and (1.6) allows us to express the ordinates of the in- 
ner and outer houndarics of the mixing zone r; and ry, In terms of the current chamber radius 
ro and the Initial values of the kinematile and geometric parameters uy and «a 


2 1, V0 WY, 


ry hy (i l, z), (2 3) 
where 
lr, A, Kiri, Q, (Aiyo— By yore)’, | 
” 2 
A (I 1,72 PV May) L(V May)? | 9h 0), (2.4) 
Ly, ad [ss My) WM) (/ !, 2, 3, 4). | 
(1 — Wy,)* 
HNere, we have 
8, OMG; s, 0,225: s, 0.52); S, °° 1 A060, 
{, 0,048; L, ONO; tl, f, 1.0; 
bh, — 0.055, hy, - O08), 


Since, as has been shown In Section 1, the cjectfon equations are different for uy < 1 and 
Uy) > 1, we will consider these two cases separately. When ugg~ 1, we express r, in terms of 
ro and 74(%) and then substitute this expression Into formula (2.3) with 1=2 to obtain an equation 


whose solution for ré yields 
ré (x) L |/ Za (2.5) 
1) fs 


hy + 2B, 4 B%, | 
Ke AGA RAS (x) 4 My (4), | 
fe Ay} AB, } (My | Bx) py (4). 


where 


(2.6) 
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Fig. 4. The variation of the radius of the isobaric 
mixing chamber at the end of the primary zone as a 
function of the parameters uy) and a, where Up)< 1. 


In order to obtain the function r,(x) for the case when the velocity of the peripheral stream 
is higher than that of the central jet (u.)> 1), we make use of Eqs. (1.12) and (2.3) with i= 1, the 
simultaneous solution of which is 


; K / LB; 
To (x)= B2 ( I+ | 1— a ) (2.7) 
3 


K =A;B;— By (x), | 
L = Aj — 2A yj» (x) +492 (x). 


where 


(2.8) 


Expressions (2.5)~(2.8) obtained above allow us to determine the shape of the isobaric 
mixing chamber from the initial values of the kinematic and geometric parameters uy) and a 
and then, with the help of formulas (1.1) and (2.3), to determine the fluid flow during the process 
of intermixing of jets in the channel at constant pressure. It should be noted that the calculation 
of flow and channel shape by means of formulas (2.3)~(2.8) is only valid for the primary zone of 
the initial region, i.e., provided that unperturbed-flow cores exist in both the central and peri-~ 
pheral jets. As soon as one of the unperturbed cores ceases to exist, the primary zone of the 
initial region terminates. It is therefore clear that formulas (2.5) and (2.7) can only be used 
provided that re > r,,, where Lop is the radius of the mixing chamber at the point corresponding 
to the end of the primary zone. 


05 
0 Qi 0.2 03 04 rod 


Fig. 5. The variation of the radius of the isobaric 
mixing chamber at the end of the primary zone asa 
function of the parameters uy) and a, where uy) < 1. 
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The magnitude of rep is completely determined by the initial values of the kinematic and 
geometric parameters Uy) and qa. It can be found from the system of equations (1.5)-(1.6) sub- 
ject to the following conditions: 


r,=0 for flow according to scheme A (see Fig. 1); 
ry =rq for flow according to scheme B (See Fig. 2). 
The corresponding formulas for the determination of the function Lop(go, a) are 


for scheme A 


A. 
Fepl(ttay 2)—= (2.9) 
for scheme B 
) P / CA 2.10 
Top ieo» 1) = Cc l —| ] — p: ) ( ° ) 


where 


C = hy + 2By+ Bi, 
P=A,+A,B,. 


The function Lep(Ugo> a) given by relations (2.9) and (2.10) is shown in Fig. 4 (for uy) < 1) 
and in Fig. 5 (for uy7>1). The dot-dash curve in both figures gives the boundary between the 
flow regimes A and B (see Fig. 3), the curves lying to the left of this line corresponding to 
scheme A and those to the right to scheme B. 


The Secondary Zone of the Initial Region. The shape of the mixing cham- 
ber and the flow parameters for the secondary Zone of the initial region are found by a method 
analogous to that used above. The only difference is that instead of one of the boundaries of the 
mixing Zone (the inner one in the case of scheme A and the outer one in the case of scheme B), 
another unknown makes its appearance — the velocity either on the axis (scheme A) or on the 
chamber wall (scheme B). 


Let us consider the case that is most frequently encountered in practice, namely, flow ac- 
cording to scheme A with uyy7< 1. To find three unknowns rg, r,, and uy, we have a system of 
three equations (1.5), (1.6), and (1.10), which under the conditions indicated above (r;,=0, Uy =Up9) 
can be written as 


2.2 
c (2.11) 


With the help of the first two equations of this system, we can express u, and Yo’ in terms 
of the chamber radius and the initial parameters 
by — Qian + ru 


c 
My = Uy 4-— 


5 ’ 
C9 ] — Tyla 


(2.12) 


2 
Co (1 — routeo)? 


22 
lo — QU) -F F409 


The simultaneous solution of the last equations of system (2.11) and (2.12) allows us to deter- 
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mine the variation of the radius of the mixing chamber along its length in the secondary zone of 


the initial region 
ro (x). - 1_( / pt) (2.13) 
Mu. | J? 


where 
Ja tn =" 4 (x), 
QU 2 
Ha= 1 4-1 (tg 29) 21 OX); 
C7 


The radius of the isobaric mixing chamber in the final section of the initial region is found 
2 bo 
1 = Q.486 —— 
- (Loy 
ne ~ —~} 1— 


from Eq. (2.13) with 7 ,=90, 
cs os (2.14) 
ayy ly . 
nn 
Mn) 


The position of the final section of the initial region is also found from the condition 7 ,=0 


bo 


Ns a 
Awl yy w [1 — (| _ a) uy] 


or, alternatively, if we make uSe of the concept of ejection coefficient N, 


3. Discussion of the Results 


The relations obtained in Section 2 allow us to determine the shape of the initial region 
of an isobaric mixing chamber from the given values of the kinematic and geometric parameters. 
An examination of Figs. 4 and 5 shows that the radius of the isobaric mixing chamber should 
decrease downstream, and this decrease will occur more rapidly as the deviation of u., from 
unity increases and as the value of wdecreases. The reason for this is that these conditions 
correspond to an increase in the nonuniformity of the velocity profile in the initial section of 


Ud; 
of 


: “ 


Fig. 6. Cross section of an isobaric mixing cham- 
ber for ug9=0.7 and ~w=0.05 (m)y=0.105, N=1.98, 
ip) =2.71). 
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Fig. 7. Cross section of an isobaric mixing cham- 
ber for Uy) =1.5 and a =0.5 (m)=3.0, N=3.0, ig = 
1.25). 


the chamber. However, it is well known (for example, See the author's paper [3]) that as the 
kinematic nonuniformity of the stream in a channel of constant transverse section increases, 
there is a greater static-pressure increase associated with the smoothing out of the nonuniform- 
ity. Consequently, if we impose the condition that pressure in the channel is to be constant, we 
find that its cross Sectional area (or radius) should decrease downstream and that this decrease will 
be greater, the larger the initial nonuniformity of the velocity profile, i.e., the smaller the 

value of a and the greater the deviation of uy) from unity. 


In conclusion, let us consider two concrete examples of a numerical calculation of the 
shape of an isobaric mixing chamber, taking u,)=0.7, @=0.05 (Example 1) and uy)=1.5, w=0.5 
(Example 2). In the first example, flow is according to scheme A, in the second according to 
scheme B (see Fig. 3). The values of mp, i, and N corresponding to the specified values of uy) 
and q@ are given in Figs. 6 and 7 which show the longitudinal cross section through the initial 
regions of the mixing chambers, as well as outlines of the mixing zones. It can be Seen that the 
maintenance of isobaric flow requires an appreciable decreaSe in the radius of the mixing cham- 
ber (especially in the first example). 


The method for determining the shape of an isobaric mixing chamber and the flow param- 
eters presented above can be used for the selection of optimal ejector parameters (for example, 
the velocity of the ejected stream u,, and geometric parameter a for a given value of the ejec- 
tion coefficient). However, the solution of this problem was not an aim of the present paper. 


The author would like to express his gratitude to G. B. Krayushkina for her help with the 
calculations. 
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THE PROPAGATION OF A TURBULENT JET 
IN AN OPPOSING STREAM 


A. N. Sekundov 


The investigation of the propagation of jets in opposing streams is acquiring great impor- 
tance at the present time in view of many technological applications. For example, the possibil- 
ity of the use of aerodynamic stabilizers in afterburners of turbojet engines is pointed out in 
[1]. The properties of this type of stabilization of combustion in the zone of development of a 
jet in an opposing stream exhibit significant advantages by comparison with the properties of 
conventional mechanical stabilizers. The design and construction of this type of device neces- 
sitates acomplete understanding of the laws of propagation of turbulent jets in opposing streams. 


One of the first to investigate this problem was L. A. Vulis [2] who made an experimental 
determination of the "range" / (the length of the circulation zone) of a jet propagating in the 
initial region of an opposing jet with a large diameter. In these experiments, the principal 
parameter m (the ratio of the jet velocity to the velocity of the opposing stream) was varied 
within the range 2.13 to 2.63. L. A. Vulis proposed a superposition method for theoretical cal- 
culations, but the method only leads to qualitative agreement between theoretical and experi- 
mental data. 


Yu. V. Ivanov, Kh. N. Sui, and E. P. Timma [3, 4, 5] have made a detailed investigation 
of the velocity field in a jet propagating in an opposing current over a wide range of variation of 
the parameter m (1 to 20). The dependence of the range of the jet on this parameter, l(m), was 
established. In varying the diameter of the jet source, the above authors noticed that the op- 
posing stream can be considered to be infinitely wide when the condition D/lI = 2.5, where D is 
the diameter of the opposing stream, is Satisfied. 


A paper by A. S. Ginevskii [6] has appeared recently containing a theoretical calculation 
of the whole flow zone. A. S. Ginevskii notes that there is good agreement between the results 
of this calculation and the experimental data given in [7]. However, these results differ from 
those of [3] and [4], for example, by approximately 20-30% with regard to the quantity 1. More- 
over, in the theoretical analysis of [6], several assumptions have been adopted without sufficient 
substantiation: Thus, the author assumes that the thickness of the boundary layer at the nozzle 
section is equal to zero and that the pressure throughout the flow zone is constant (See Fig. 1). 


Thus, the problem of the propagation of a jet in an opposing stream has not yet been 
solved unambiguously, since there are considerable variations between the results of different 
authors. 
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Fig. 1. The propagation of a jet in a cylindrical 
chamber. 


In the following, an approximate solution of the above problem is obtained as a special 
case of a more general problem of the propagation of a jet in an opposing stream flowing in a 


cylindrical chamber. 


1. The Propagation of a Jet in an OppoSing Stream 


Let us consider the propagation of a jet in a cylindrical chamber (see Fig. 1). A jet is- 
sues with velocity uy from a nozzle of radius ry and propagates along the axis of a cylindrical 
chamber of radius Rj; the jet expands as it is retarded by the opposing stream whose velocity 
Uco is assumed to be less than the velocity uy of the jet. There exists a closed region bounded 
by the surface r(x) in which the jet flow retains its initial direction. The boundedness of this 
region is the most typical feature which distinguishes the flow scheme under consideration 
from the usual cases of jet propagation. The surface r(x) separates regions of unperturbed and 
perturbed flow in the opposing stream (the unperturbed flow has a uniform velocity field) and 
r;(x) is an analogous surface within the main jet. Experiments show that the whole region of 
flow in a cylindrical chamber can be subdivided in a longitudinal direction into two zones. In 
zone I, which adjoins the nozzle, the propagation of the jet is qualitatively similar to that of the 
usual submerged jet, namely, the flow rate and the thickness of the jet r.(x) increase with dis- 
tance away from the nozzle, the boundaries of the jet are practically linear and subtend an angle 
a =arctg 0.22 with the chamber axis, and the static pressure is approximately constant. In 
zone II, the magnitudes of r.(x) and static pressure change sharply — the final retardation of the 
jet by the opposing stream takes place. Let us adopt as the transverse section separating these 
two zones the section in which the jet velocity uy;, becomes equal to the velocity u, of the op- 
posing stream. 


On the basis of the data given in [4], we can show that in this section the thickness of the 
jet rp also reaches its maximum value. 


We will use a method described in [8] to obtain a closed system of equations. Let us 
write the equation describing the variation of the thickness of the jet in the main region as fol- 
lows: 


Ae I) py, (1.1) 


where v' 1s the transverse fluctuation velocity, ugy the average velocity of the jet, and k(x) a 
correction factor which takes into account the various features of jet propagation in an opposing 
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stream. Let us discuss the meaning of this coefficient in greater detail. As was already pointed 
out above, in section P (see Fig. 1) we have dr,/dx =0, whereas in the case of jet propagation, 
for example, in a parallel stream we have dr,/dx> 0. This behavior of the derivative dr,/dx 
can be explained as follows. In the general case, the increase in the thickness of a turbulent 

jet occurs because particles of the jet being carried by turbulent fluctuations into the surround- 
ing medium entrain the matter of this medium and impart motion to it, as a result of which mat- 
ter from the surrounding current becomes added to the jet. If the external current moves ina 
direction opposite to that of the jet, this process can only occur if there are particles in the jet 
whose velocities are higher than those of the particles in the external medium. It follows from 
this, that the broadening of a jet propagating in an opposing stream is governed by that compo- 
nent of the jet whose velocity is greater in magnitude than the velocity of the opposing stream 
uy. Let us denote the ratio of the "active" mass of the jet to the total mass by 


ra 


5 Qi i urdr 
Q) 
+ == ——__—_-— 1,2 
G r: , ( ) 
20 | urdr 
0 
where ra, is the radius at which the jet velocity is u=u,. Let us further assume that 
‘GCG n 
R \a=-{ —2t 1.3 
(x)= (2) (1.3) 


where n is an empirical constant whose value according to preliminary determinations is 0.25. 
Following the method described in [8] and assuming that |v'|~u,, + u, and ugy~uUp, we finally 
find that 


Ory Um (2) (1.4) 


ax Um G 


The constant c can be found if we proceed to the limit ux, —0 in expression (1.4) 


so that, as in this case the jet reduces to the usual submerged jet for which 


fs 9,99) (1.5) 


ax 
it is reasonable to assume that c=0.22. 


The law governing the broadening of of the turbulent boundary layer in the initial region 
of zone I can be obtained in the same manner as in the case of the main region, but the varia- 
tions of the parameters uy), G, and Gg in the initial region are such that the right-hand side of 
expression (1.4) is close to the constant c and since we have c=0.27 for the initial region of a 
submerged jet, we can assume that 


MM =) 997. (1.6) 


ax 


In order to be able to calculate the jet in an opposing stream, we must know the velocity 
profile in the zone of intermixing of the jet and the opposing stream. Several empirical expres- 
sions of varying accuracy are given in [3-6] for the description of the true velocity profiles. 
Calculations have shown that in view of the fact that the equations of motion are in integral form, 
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the final results depend only weakly on the expression used for the velocity profile. Therefore, 
to simplify calculations, we will use a simple velocity profile in what follows; this consists of 
two linear sections as follows: The velocity profile from the jet axis to ro(x) is 


Mo Pe (1.7) 


Um Ps 


while the velocity profile in the opposing stream from ry(x) to r3(x) is 


— To 
He FT . (1.8) 
il Py Fa 


* 


We will assume that in the initial region the velocity profile is also linear 


ear (1.9) 


{Ly ry — Fy 


2. The Propagation of a Turbulent Jet in a Cylindrical Chamber 
and an Infinite Opposing Stream 


Let us now proceed to the solution of the problem of the propagation of an axially sym- 
metric jet in an opposing stream inside a cylindrical chamber (see Fig. 1). The opposing stream 
in flowing around the broadening jet is accelerated to a velocity ux. According to Bernoulli's 
equation, the pressure in it decreases according to 


Po Pp, = (11° —u’). (2.1) 


Measurements of the velocity ux and pressure px have shown that within zone I these 
parameters vary only slightly, so that we will assume that velocity and pressure are constant 
Within Zone I. 


The equations of conservation of flow rate and momentum for the contour bounded by sec- 
tions O and V can be written as follows: 


fay 
| 2 I 9 . I 2 
> Moro + =u. Ro—\ urdr— ——U, (Ro — ro3)= 0, (2.2) 


re 


l ] C ° | 9 9 

3 Mara — Wi Ro — \ u-rdr We (Ro ros) +E Ri (2.3) 
Where Ap=p,. — px and rpg is to be evaluated in section 0, i.e., at x=0. Transforming these 
equations together with Eqs. (1.8) and (2.1), we obtain 


_ 3 (m + R?) 
ss 


7 — = , (2.4) 
BRA — (ra +73 + 1) 


= Qm? + Re (u? — | 
roa= —0.333-++ YV — 0.2924" 7 SoM (2.5) 
u 


* 


where 


Uo > R 
m==— and Ay=—. 
oe ry 
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For given values of m and Ros this system of equations can be 
solved for u=u, /u,, and rpg =Yo3/r, by the method of successive ap- 
proximations. Calculations have shown that uy increases with in- 
creasing m and decreasing Ry, whereas ry3 increases with increas- 
ing m and Rp. 


In addition to the obvious condition 


L<r3(m, Ro) < Ro, (2 .6) 


the admissible values of the parameters m and Ry are also subject 
to one other condition which follows from the fact that the flow 
scheme under consideration cannot exist if the total pressureinthe 
jet 


Fig. 2. The region of 
admissible values of ou? 
the parameters m and Po=P, +A 
Ro. : 


is less than the total pressure in the opposing stream 


9 
Aue ou; 


QO ’ renN 


This yields the following restriction on the possible values of the parameters m and R;: 
u,(m, Ro) <m. (2.7) 


Relations (2.6) and (2.7) allow us to determine completely the region of admissible values 
of the parameters m and R, (see Fig. 2). The solution of the system of equations (2.4) and (2.5) 
simplifies considerably if Ry ~ ~, which corresponds to the case of propagation of a jet in an 
unbounded opposing stream 


= 1, (2 .8) 


Fogz=) 6m(m+1)+1. (2.9) 


In particular, it follows from expressions (2.8) and (2.1) that in the propagation of a jet 
in an infinite opposing stream, the pressure in zone lis py =p,. An "infinite" opposing stream 
is produced in practice by a jet with a sufficiently large diameter D,). The magnitude of this 
diameter on the basis of the flow scheme under consideration should satisfy the inequality 
Dy > 2Y¥o3- This condition, together with relation (2.9) and expression (2.32) obtained below can 
be written as 


“) > 2. (2.10) 


It should be noted that condition (2.10) is in good agreement with the experimental result 
obtained in [4]. 


Let us now consider the main region of zone I. The equations of conservation of flow rate 
and momentum for a contour bounded by a section through the main region and section V can be 
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written as 
r2 ry 


. | 2 ; ] 9 a 
\ urdr - ~- . Ro -- \ urdr — a, (Ao — 3) - 0, (2.11) 


0) Is 


f2 ry 
yo, }) 39? oy 1) 2 yp 9 A 
\ aordr -- “Ou i \ u-rdr — ou, (Rg — 3) +> f 


0 2 


FR, (2.12) 


4 
I 


Let us transform Eqs. (2.11) and (2.12) with the help of relations (1.7), (1.8), and (2.1). 


into 
- Bake 
N? 3-N 4-3=—_ , (2.13) 
= 
(3-6) NPA (2 b)N-=3(a—b)= 0, (2.14) 
where 
V=- ry ; ee 4 je Q— Uy — | h Heed (2.15) 


Po u, i, uy ul 


in which U,, =U,,/uU.. is the dimensionless velocity on the jet axis. 


This system of equations allows us to determine r,(Um, Ro, Ux) and Fs (Uy, Ro, Ux), where, 
as was pointed out above, u, is found from relations (2.4) and (2.5). If in these functions we 
set Um =m, we will obtain the values of r, and ry in the transition section H between the initial 
and main regions: Yo} =Y(m, Ro, Uy) and Y3,, =F3(m, Ro, Uy). 


If we now make uSe of Eq. (1.6), we can easily find the length of the initial region of zone 
I, namely, 


] 


9 2 ie 0.97 Poff. (2.16) 


Since the parameters in the initial region vary approximately linearly, then, knowing their 
values at the beginning and end of the initial region, as well as the length of the latter, we can 
find them for any intermediate section. 


In order to calculate the main section, we must solve the system of equations (2.13) and 
(2.14) for um(Lo, Ro, Uy ), substitute into the resulting expression the value of u,(m, Ro), and 
then integrate Eq. (1.4) which, together with expressions (1.2) and (1.7), can be rewritten as 


are tm + ty t ~3( cP )+2(2))" (2.17) 
ax Um Um 


Um 


Since it is impossible to obtain an explict expression for Um as a function of r,, the inte- 
gration of this equation is difficult even numerically since the right-hand side of (2.17) tends to 
zero aS Um — U,. AS a rough approximation we can take 


fe og, (2.18) 
ay 
since the right-hand side of expression (2.17) is close to c over most of the range of variation 


of Up, from m tou ,: In addition, an indirect confirmation of the validity of this assumption can 
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be obtained as follows. It follows from Eq. (2.18) that the maximum thickness of the jet in the 
transverse section P is 


Top =CXp +00, (2.19) 


where the quantity cy =0.05x,;, can be neglected when m is large. Making use of relations (2.27) 
and (2.31) which will be derived below, we can determine the relative width of the jet, which is 
found to be independent of the parameter m, 


or 


~~ . 


-=const— O.33. (2.20) 


5 


This result is confirmed by many experiments, although the numerical value of the con- 
stant differs from experiment to experiment; thus, this constant is found to be 0.36 in [5] and 
0.43 in [8]. 


Equation (2.18) allows us to find the dependence of r, on the distance from the end of the 
nozzle 


—— 


r>—-ron =0.2 2(x—x 71). (2.21) 
The system of equations (2.13), (2.14), and (2.21) which is a system of three algebraic 
equations with three unknowns ry, r3, Uy, and two parameters m and Ry can be easily solved. 


From the condition that Um = Uy at the end of zone I, we can determine the length of the 
main region 


—— Xy) => [rap —Ly)). (2.22) 
Calculations have shown that this quantity increases with increasing values of the param- 


eters m and Ro. 


The solution of system (2.13), (2.14), and (2.21) can be obtained in explicit form if we con- 
sider the case Ry ~~ which corresponds to the propagation of the jet in an infinite medium. 


Making use of relation (2.4) and eliminating the indeterminacy in Eq. (2.13) which arises 
as Ry — ©, we obtain 


—9 2 ¢ -I. : + 1) + 
rp St + Om + 24 | mim ed etm) (2.23) 
wb bum + 3 Y (Um) 
Fy= (iim I) 72 C24) 


These relations allow us to find the lengths of the initial and main regions very simply if we 


take account of Eqs. (2.16) and (2.21). 
2 (my 
/ 2.2 
0.27 ( °) 


u (mt) 


~~ (ut) ¢ (im) 2.26 
oe LW io y/o, P20 
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where 


«(m)= 6m? +-9m +2--] Om (m-——-1)+1, 
') (m)= mm? +4m+3. 


It is not possible to neglect pressure 
variations in zone II as was done in the 
case of zone I, since the pressure here is 
comparable with the velocity head. The 
velocity profiles likewise change very 
rapidly. Therefore, it is not possible to 
perform an exact calculation of the flow in 
zone IJ. Let us show how the dimensions 
of this zone can be estimated approximate- 
ly. In view of the fact that the velocities 
of the direct and opposing streams in Sec- 
tion P are equal (uy, =u, ), we can assume 
that the development of the jet in zone II 
occurs as if without energy losses. This 


Fig. 3. Comparison of calculated and experi- 
mental values of the range / for various values 
of the parameters m and R;; 


x R,)— «, according to the data of [7]; assumption allows us to consider flow in 
d R,— ~, according to the data of [4]; zone II to be potential flow and to make 

O Ro | according to the author's use of the method described in [5] for the 
® Ry=3 experiments calculation of the circulatory zone behind 
O R,=15 a badly streamlined body. As a result of 


such a calculation we obtain 


Xx — Xp ~~ 156+ rap == 1.56 / a ° (2.27) 


We will assume that this dependence remains valid when R, ~ ~, as well as when R, > ~. 


Before we proceed to the calculation of the "range" of a jet in an opposing stream, let us 
point out an important circumstance. As can be seen, for example, from formula (2.26), when 
m— 1 the width of the main region in zone J is equal to zero, while the length of the initial re- 
gion, as determined from expressions (2.16) or (2.23), will be a small but nevertheless non- 
zero quantity. However, it is obvious that with m=1, the condition u,, =u, holds directly at the 
end of the nozzle and zone J will be completely absent, so that the "range" of the jet will be 
governed only by the dimensions of zone JJ. This contradiction is apparently due to the arbitrary 
nature of the flow model used and can be eliminated through the introduction of a correction 
coeffic ent 


X=k(mM) Xj; (2.28) 
where 
A O, when m=1, (2.29) 
(m) =| when m > oo. 


Since we do not have any auxiliary conditions for the determination of the function k(m), 
we will assume that 


k(m)= mot (2.30) 


ii 
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Thus, the "range" of the jet is given by 


Lok (m) Xp 4+-(Xp — Xj)“ (X4— Xp). (2.31) 


_ Expressions for the individual terms in this equation have been derived above. The quan- 
tity / depends on the parameters m and Rj; this dependence is shown in Fig. 3. It has been found 
that for Ry) > 100, the value of J is insensitive to changes in Ry so that only one curve has been 
shown in the figure for values of Ry> 100. Thus, when Ro — oo, the "range" of the jet is governed 
only by a single parameter m, 


7 (m) = 2.16 bo (m) — (4.55 —3.71 | \/ ee | (2.32) 


mt Ww UN 


In addition to this dependence, Fig. 3 also shows the results of "range" calculations for 
R,=30, 15, and 10. It can be seen from the figure that these curves terminate abruptly at large 
values of m, since here we fall into the region of inadmissible values of m and Ry (see Fig. 2). It 
should be noted that there is satisfactory agreement between the experimental data of Kh. N. Sui 
and the author, and the results of the calculations. 


3. The Additional Resistance Created by the Opposing Stream 


Cases of practical application of the flow scheme under consideration, for example, 
aerodynamic stabilization of combustion in afterburners of turbojet engines, have been men- 
tioned above. The flow of a jet in an afterburner is of course nonisothermal. The calculation 
carried out above neglected density differences between the material of the jet and the opposing 
stream. Let us make an approximate estimate of the influence of nonisothermal conditions on 
the range of the jet. It has been shown in [9] that the following expression holds for the velocity 
variation along the axis of a submerged nonisothermal jet: 


Ll, 
~~ 


il ~, 
xv 6 


(3.1) 


m 


where 6 = Poo/Pos Po iS the density of the opposing stream, py that of the jet. In view of the fact 
that in our scheme the law of flow in zone I which forms the major part of the total length of the 
jet is similar to the law of propagation of a submerged jet, we can assume on the basis of ex- 
pression (3.1) that 


I 


1, (3.2) 
Vo 


— 
—_— 
—_ 


where L is the length of a nonisothermal jet. 
This result is in satisfactory agreement with the experimental data of Yu. V. Ivanov. 


In conclusion, let us determine the coefficient of hydraulic resistance € in a cylindrical 
chamber containing an opposing axial stream (neglecting friction with the chamber walls), 


( ou. ( sata | 
0 0 ot — {4-7 
p_ Px Pa —pa_\ie So MA 2 ua 2 (3.3) 


2 2 


2 2 


where Doo” and DA are the total pressures in sections at infinite distances to the left and right 
of the nozzle. It is assumed that the velocity and density fields are uniform in section A. In 
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order to find the unknowns appearing in expression (3.3), 
we ean write down the equations of conservation of mo- 
mentum, flow rate, and excess enthalpy for the contour 
bounded by seetions A and B, namely, 


9 y 209 2 2 2 
ayltoro — 9, URo= pala (Ro— 10) 


-(P.~ po) ro +(P, —pa)(Ro—10), (3.4) 
intl Ua RO= Fl 4 (RO— TO), (3.5) 
OU RST. To) oat (Ro—ro)(T 4-7). (8.6) 
0 0, 005 9.010 gE Simple calculations show that with typical values of 
Fig. 4. Drag in a cylindrical the various parameters: Ro =20, m=5, p,=1 atm, Woo = 
chamber as a function of the 100 m/sec, the ratios poTo/paTa and paT a/p,,T,, differ 
paramcters g, m, and 9. from unity by less than 2%. Therefore, Eq. (3.6) can be 


rewritten as 
u.Ro(1— 0)=u, (Ro —1)(1—9,), (3.7) 


where 0 =P./Po Pa =p a/Po- We can dctermine Pa> Pa» and ua by a simultaneous solution of 
Inqs. (3.4), (3.5), and (3.7). Substituting the valucs of these quantities into (3.3), we obtain 


: R- [2a bE (1 Om] — 20p. — pa) (Re —!)  (m=—0) 3 3) 
W (RW 


The last two terms in the numerator can be neglceted when Ry= 20. If we introduce a 
parameter g characterizing the rclative flow rate and defined by G,/G..=m/oR%, the final ex- 
pression for é will be 


pm (LY Ome (3.9) 
(11 — Og): 


This dependenec is shown in Fig. 4. Since the flow velocities in combustion chambers 
quoted above are small, compressiblity effects have been neglected. Therefore, the total-pres- 
sure-loss in a eylindrical chamber can be calculated from the expression 


fa 
0 
Ps l—<a 
- = ——___—_ = --——_ -— , 3.10 
Pin 1— x(a.) ( 
2p 
from which wc have 
o=1—'[1—<(A,)]. (3.11) 
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AXIALLY SYMMETRIC SUPERSONIC TURBULENT JETS 
DISCHARGED FROM A NOZZLE WITH UNDEREXPANSION 


Chiang Che-haing 


Free turbulent jets emitted from axially symmetric supersonic nozzles under nonsimilar- 
ity conditions (pressure at the nozzle orifice not equal to the pressure of the surrounding medi- 
um) have not been extensively studied although they are of great practical interest. 


If the jet flows out of an axially symmetric or planar supersonic nozzle into a uniform 
supersonic stream and has the same pressure as that of the surrounding medium, then only 
weak oscillations occur, but density discontinuities are absent. Such jets do not possess aperi- 
odic or almost periodic structure [9, 11]. 


A jet flowing out of a supersonic axially symmetric nozzle with underexpansion into a sta- 
tionary medium or a uniform subsonic stream has a fairly complicated structure. Because of 
the distortion of the boundaries of the jet and expansion waves, the compression waves become 
superposed on one another on account of reflection from the jet boundary. As a result of this, 
a system of density discontinuities appears in such waves (see Fig. 1). When the coefficient of 
nonsimilarity n=pexit/PH is only slightly higher than unity, the jet has a periodic or almost 
periodic structure. Here, pexit iS the pressure at the nozzle orifice and py the pressure of the 
surrounding medium. 


Most of the papers devoted to the problem under consideration contain few experimental 
data on the first "barrel" of the jet (the zone situated between the nozzle and the first system of 
shock waves) obtained in the case of outflow from an axially symmetric nozzle. 


One of the first such papers is by Winckler [14] published in 1948. This is a purely ex- 
perimental paper devoted to a quantitative investigation of the first barrel of a jet. The ve- 
locity fields in the first barrel of an axially symmetric supersonic jet were calculated from 
the density fields measured by means of an interferometer; in these calculations, it was assumed 
that in front ofthe central forward discontinuity, the temperature, pressure, and Mach number 
are adiabatic functions of the density. Measurements of the total-pressure rise and relative 
measurements of the velocity profiles near the nozzle orifice by means of an interferometer 
confirm the validity of the last assumption. The measurements show good agreement with one 
another and this means that the intensity of the oblique discontinuity extending from the nozzle 
orifice is low, so that the entropy of the jet remains approximately constant in front of the for- 
ward diseontinuity. 


1 This paper is based on the results of work performed by the author at the Moscow Aviation 
Institute under the supervision of G. N. Abramovich. 


111 


112 CHIANG CHE-HAING 


Fig. 1. Shadow photograph of a supersonic jet. 


This assumption has also been made in other similar experiments. 


Wang and Peterson [12] in 1958 used the method of characteristics to calculate the field 
of Mach numbers in front of a system of compression waves. BecauSe of the interaction be- 
tween the last rarefaction characteristic leaving the nozzle orifice and the first compression 
characteristic reflected from the jet boundary, a density discontinuity is produced, so that the 
exact calculation of the flow field behind the discontinuity becomes very difficult. However, an 
approximate method for the calculation of the points of intersection of the density discontinuity 
and the mesh formed by the characteristics by means of a digital computer has been proposed 
by Love [9]. If the operating regime of the nozzle differs only slightly from asimilarity regime, 
then the boundary of the jet remains almost linear and with a small error it can be considered 
to be a similarity regime for 0.8<n< 1.2. If the operating regime of the nozzle differs mark- 
edly from a similarity regime, a system of density discontinuities appears in the jet and this 
causes the jet boundary to become curvilinear. At the present time, there is still no analytic 
method for the determination of the jet boundary in the case of outflow from an underexpanded 
nozzle. The theoretical boundary obtained numerically from the method of characteristics can 
be taken as the average line of the turbulent boundary layer. However, the method of character- 
istics requires a large amount of computation. 


A particularly simple method for the determination of the curve giving an approximate 
boundary of the jet has been proposed by Love [9] on the basis of the method of characteristics. 
In this method, it is assumed that the boundary of the initial region of the jet can be approxima- 
ted by a circle that passes through the nozzle orifice and the point where the diameter of the 
boundary is a maximum. The center of this circular arc lies on the normal to the jet boundary 
at the nozzle orifice and is situated at a distance 6,=dexjt from it. Love's paper contains 
families of curves for the determination of the radius of the arc as a function of n, Me it, and k. 
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1. Qualitative Investigations of Jets Flowing out of 
an Axially Symmetric Nozzle with Underexpansion 


1. Distribution of Pressure Along the Jet Axis. ‘The experimental data 
on the distribution of static pressures along the axis of a jet flowing out of an axially symmetric 
nozzle under nonsimilarity conditions are shown in Figs. 2-5 for Meyit equal to 1.5 and 2.0, 
respectively. It can be seen from these graphs that in a jet flowing out of an underexpanded 
nozZle with n > 1, the static pressure on the jet axis falls away from the nozzle orifice to the 
discontinuity and then, after the discontinuity, it does not differ appreciably from the pressure 
of the surrounding medium. Subsequently, the static pressure on the jet axis varies according 
to a damped sinoidal curve and at some distance from the nozzle it becomes equal to the atmo- 
spheric pressure; in other words, a jet flowing out of an underexpanded nozzle periodically ex- 
pands and contracts; because of energy losses in the density discontinuities, the excess pres- 
sure in the jet tends to zero. 


It should be noticed that with n>1 these oscillations on the jet axis rapidly die out (with 
the exception of the case n= 2) as a result of the strong forward discontinuity in the first bar- 
rel. With the exception of the first phase, deviations of static pressure on the jet axis from 
the pressure of the surrounding medium are less than 10%. The larger the value of n, the more 
rapid the attenuation of the pressure oscillations along the jet axis. When n#=2, the deviation 
of the pressure on the axis from atmospheric pressure is less than 20% during the first phase. 
If the value of n is only negligeably different from unity (either greater than or less than unity), 
the deviation of pressure on the jet axis from the atmospheric pressure is also small and less 
than 10%, although the pressure oscillations extend to greater distances from the nozzle ori- 
fice. This phenomenon can be explained as follows: With n large, an intense central forward 
discontinuity arises behind the first barrel and this leads to a large total-pressure drop. If 
nis negligibly different from unity and the mean pressure in the jet is only slightly above the 
pressure in the surrounding medium, then strong discontinuities do not arise and total-pressure 
drops are small. As a result of this, the static-pressure oscillations are damped out only slow- 
ly. The case n=2 should be considered as a transition regime in which the jet already contains 
a central discontinuity, but its diameter is small and rarefraction in the jet before the discon- 
tinuity is not very great. This circumstance leads to a considerable increase in the pressure in 
the discontinuity which nevertheless decays along the jet axis at a slow rate (see Fig. 2). 


o Experiment 
— Theoretical curve 


| 
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Fig. 2. Distribution of static pressure along the 
jet axis, Meyjt=1.5, n=2. 
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o Experiment 
—— Theoretical curve 
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Fig. 3. Distribution of static pressure along the 


jet axis, Mexjt=1.5, n=2. 


The experimental curves of the distribution of static pressure along the jet axis can be 
described by expressions for damped oscillations whose period and amplitude depend on the 
Mach number at the nozzle exit, Mexjt, and the degree of nonsimilarity n. Some of these calcu- 
lated curves are shown in Figs. 2 and 3 and these were obtained from the formula 


Py = Ae-x gin wx. (1.1) 
Pu 


The following values of A,a@ , and w were chosen for the calculations: 


Mexit 


= 


1.5 | 9.0 
2 3 4.09 5 5.45 ma | 2.56 
0.55 0.22 0,22 0.22 0.22 0.3 0.24 
0.1 0.1 Q,] QO. 1 0.1 Q.1 | 0.1 
a mt mw ne = a | a 
| 2.5 | 2.9 ) 2.5 | 2.5 2.5 4 2.0 
PL. 
Pu 
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Fig.4. Distribution of static pressure along the 
jet axis, Meyit =2; Po/PH =29; Pexit/PH = 2.56. 
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We were unable to find the theoreti- 
cal expressions for the coefficients A, a, 
and w as functions of the degree of non- 
similarity n and the Mach number at the 
nozzle exit, Mexjt- 


The experimental points shown in 
Figs. 2 and 3 are in satisfactory agree- 
ment with the curves obtained from for- 
mula (1.1). 


2. Determination of the 
Length of the First Barrel 
of the Jet. The lengths of the first 
barrel of a jet were measured on 


Position of discontinuity 


x | 


5 10 IS y 
é 2 shadow photographs obtained with the 
Fig. 5. Distribution of static pressures along IAB-451 instrument for various regimes 
the jet axis, Mexjt=2.0, n=4. of a jet from a nozzle. Figures 6-8 show 


the variations of the relative lengths w as 
functions of the degree of nonsimilarity n for three values of the Mach number at the nozzle 
orifice. Here, w is defined by 


(1.1la) 


where w is the length of the first barrel of the jet. 


These figures give the experimental points obtained by the author and the results of [9]. 
It can be seen from these figures that at a given value of Mexjt, the length of the first barrel 
of the jet increases with increasing air pressure inthe receiver. The rate at which the length 
of the first barrel increases with increasing value of n grows with increasing value of the 
Mach number Meyjt at the nozzle orifice. At the present time, there is no method for the deter- 


9 Author's experiments 
e NASA experiments 
——. Empirical curve from 
Love's formula 


—— Empirical curve from 
formula (1.2) 


0 10 20 JO 40 n 


Fig. 6. Variation of the relative length w as a func- 
tion of the degree of nonsimilarity n at the nozzle 
orifice for Meyjt=1.- 
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° Author's experiments 
e NASA experiments 


_— Empirical curve from 
Love's formula 


—— Empirical curve from 
formula (1.2) 
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Fig. 7. Variation of the relative length w as a 
function of the degree of nonsimilarity n at the 
nozZle orifice for Mexjt =2. 


mination of the theoretical relation between the quantities w, n, and Mexit. However, the experi- 
mental data reveal that for n>2, the variation of w with nis linear on a logarithmic scale. The 
gradient of this line depends on M,,jt- This fact allows us to find the following empirical for- 
mula for the determination of the length of the first barrel for any values of n and Me yjt: 


when n=2 
1 1 1 
— “5 9 “> 9 oT, 9 y 
W= 3.1M3,1 (M2, —1)? —(Maxie—1) ° i 2(Mixit—1)°, (1.2a) 
when n= 2 
a 1) 
_- y 7 “) , 9 ~~ 4 t 
o—(3.1Mea (nMexit —1) > — (Mexit 1) ° | +.2(Mexit —1) J (=) | (1 .2b) 
WwW 
28 
26 
24 
22 
20 
18 
16 
14 
12 
0 Author's experiments 
6 +—< e NASA experiments 
6 —— Empirical curve from 
4 Love's formula 
> — Empirical curve from 
formula (1.2) 

0 J 10 15 20 n 


Fig. 8. Variation of the relative length W as a func- 
tion of the degree of nonsimilarity n at the nozzle 
orifice for Mexit =3. 
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Formula (1.2b) differs from (1.2a) only by multiplier (n/2)°. The magnitude of the expo- 
nent t depends on Meyit as follows: 


- 0.523 
when M exit < 1.5, {= —_—_. , (1 .3a) 
} Mexit 
when Mexit > 1.9, f==0.491 —O.016M.yi - (1.3b) 


In addition to the experimental points, Figs. 6-8 give the calculated curves for the varia- 
tion of w with n obtained from formula (1.2) for values of Mey; of 1.0, 2.0, and 3.0. These 
curves are in good agreement with the experimental points. 


An empirical formula for the determination of w as a function of n and Mex; has also been 
obtained in [9]. A comparison of the author's empirical curve with Love's empirical curve (see 
Figs. 6-8) shows that the empirical formula (1.2) is in better agreement with experimental data 
than Love's formula with the exception of outflow from a sonic nozzle. In the latter case, the 
two curves coincide when n=20, while for n>20 and Meyjt=1, empirical curve (1.2b) shows 
worse agreement with experiment than Love's curve, although formula (1.2) can be used for 
practical calculations even in this region. 


3. Determination of the Distance between the Orifice and the 
Central Forward Discontinuity in the First Barrel of the Jet. Fig- 
ures 9-11 show the experimental data on the variation of the distance between the nozzle orifice 
and the central forward discontinuity, lo =1./Rexit, aS a function of the degree of nonsimilarity 
n for the following values of the Mach number at the nozzle orifice: Mexjt=1.0, 2.0, and 3.0. If 
the operating regime of the nozzle is appreciably different from a similarity regime and there 
is no direct forward discontinuity in the jet, then the distance J, is measured from the nozzle 
orifice to the point of intersection of the oblique discontinuities (see Fig. 9). The variation of 
thé relative quantity lc=Ie/Rexit With n and Mexjt is similar to that of w. 


It was possible to obtain a simple relation holding between Jc and W for fixed values of 
n and Meyit> 


1, =0.8 w. (1.4) 


o Author's experiments 
—..— APL experiments 
—— NASA experiments 
—— Empirical curve accord- 
ing to formula (1.4) 
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Fig. 9. The relative length of the first barrel le as a func- 
tion of the degree of nonsimilarity n for Meyjt=1. 
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o Author's experiments 
— — NASA experiments 


—— Empirical curve according 
to formula (1.4) 


0 P) 10 19 

Fig. 10. The relative length of the first 
"barrel" le as a function of the degree of 
nonsimilarity n for Mexjt =3. 


o Author's experiments 
—.— APL experiments 
— — NASA experiments 


——— Empirical curve ac- 


cording to formula 
1.4 


Y J 10 Sn 
Fig. 11. The relative length of the first 
"barrel" Jo as a function of the degree of 
nonsimilarity n for Me yjt=3. 


Fig. 12. The magnitude of the ordinate 
of the inner boundary relative to the 
mixing zone for nonisothermal gas. 


The theoretical curves obtained from this 
formula are shown in Figs. 9-11. When My yjt= 2, 
the agreement between the theoretical curves 
and the experimental data is adequate. When 
Mexit > 2, the agreement between the calculated 
curves and the experimental points is somewhat 
worse, but the theoretical curves still reflect 
the general behavior of the experimental points. 


4. Determination of the Boun- 
daries and Thickness of the Boun- 
dary Layer in the Initial Gas-Dy- 


namic Region of a Jet Issuing from 


an Underexpanded Nozzle. If we are 
dealing with a perfect gas, then the boundary of 
a jet emerging from a nozzle in a similarity 
regime must be a straight line parallel to the 
axis. In fact, a turbulent boundary layer makes 
its appearance. In the case ofthe initial region 
of a submerged jet emerging from a nozzle in 
the similarity regime, the rate of increase of 
the boundary-layer thickness can be described 
by the following approximate formula [4]: 


7 1+ fs 
AO == ().27 (1.95) 


dx ? 


where b is the thickness of the boundary layer, 
py is the density of the surrounding medium, 
and Py is the density of gas in the constant- 
velocity core. 


G. N. Abramovich [4] has obtained curves 
for the determination of the boundary layer of 
a jet emerging from a nozzle under similarity 
conditions as a function of @ (here, 9=T> /TH 
where Tj is the retardation temperature in the 
receiver and Ty the temperature of the sur- 
rounding stationary air). One of such curves 
for a submerged jet is shown in Fig. 12. 


If it is assumed that pressure variations 
are small in the turbulent boundary layer of the 
initial region of a jet emerging from an under- 
expanded nozzle, i.e., that p py, then the width 
of the jet can be calculated from formula (1.5) 
written as 


*.=0.27 fo, (1.6) 
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where b is the thickness of the boundary layer of the initial region of a jet emerging from an 
underexpanded nozzle (this thickness has been measured in a direction perpendicular to the 
curvilinear boundary of the jet calculated by the method of characteristics), s is the length of 
the theoretical boundary from the nozzle orifice to the cross section under consideration, py is 
the density of the surrounding medium, and Pp is the gas density on the theoretical boundary 
calculated by the method of characteristics. 


Since 
wm 
ro] A= Me 
then (1.6) can be rewritten as 
4° 9.135 [ 1-}+-——_% 7. (1.7) 
ds k— | 9 
] + a Ms 
where My is the Mach number on the boundary 
ka1 i 
9 Po k \? 
M,— {—— Po\" ya fl’, 
lk 1 (a) I Ce) 
Since 
FP 
Po Po k—| 2 eo) 
Po yp Pen Mexiq (1.9) 
Pu Pexit 2 
> teh > 
Me= {7 | k (145 Mexit)— 1 |} - (1.10) 


After some transformations, formula (1.6) becomes 


= 0.135 [14+ —-———_ |. (1.1) 


—1 
=: k—1 
n* (1 + =F Mbyaig 


The rate of increase of the thickness of the boun- 
dary layer in the initial region of a jet emerging from 
an underexpanded nozzle depends only on the operating 
regime of the nozzle. If the regime is kept fixed, then 
the dependence of b on s is linear. 


QO 2468 W124 16 1820 22% TN In order to determine the theoretical boundary of 
Fig. 13. The dependence of the the jet, i.e., the boundary for perfect gas emerging 
arc radius for the theoretical from an axially symmetric nozzle with underexpansion, 
jet boundary on the degree of we can make use of the approximate method proposed in 
nonsimilarity n for several val- [9]. As was mentioned above, the essence of this method 
ues of the Mach number Mexit- is that the jet boundary in the gas-dynamic region ob- 


tained by the method of characteristics can be approxi- 
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mated by a circular arc passing through the nozzle 


ed) 

Q 

S Outer ary and a point at which the diameter of the boundary 

> Theoretical obtained by the method of characteristics is a maxi- 
N boundary — mum. The center of this arc lies on the normal to 
O Circular arc 

Z, 


Inner boundary the jet boundary drawn from the nozzle orifice. 


Arc center 


Fig. 14. Schematic diagram illustrat- 
ing the determination of the jet boun- 
dary. 


Outer boun- 


oy dary 

O 

ra Arc 
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wo ary 
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ZX 


—— Theoretical 
boundary 


~~~ Experimental 
boundary 


Fig. 15. The boundary layer of a 
jet emerging from a nozzle; Meyit = 
1,n=15, g=1. 


The dependence of the arc radius on the operating 
regime of the jet taken from Love's paper [9] are 
shown in Fig. 13. The following notation has been 
used: 


r is the arc radius; 


6 exit the angle of inclination of the boundary 
at the nozzle orifice calculated by the Prandtl— 
Meyer theory [7, 9]. 


A schematic diagram illustrating the method 
used for the determination of the boundary of a jet 
emerging from an axially symmetric underexpand- 
ed nozzle is shown in Fig. 14. Knowing the operat- 
ing regime of the nozzle, we can use the graph of 
Fig. 13 to find r and then uSe Fig. 14 to determine 
dexit- Using these values of r and dexit, We Can 
construct the circular arc which will represent the 
theoretical jet-boundary. We can then use formulas 
(1.11) and Fig. 12 to determine b/s and y,/b, where 
y;, is the distance between the inner boundary of the 
boundary layer and the boundary of an ideal jet 
calculated by the method of characteristics. 


AS a result of this, for each value of the arc 
length s, where S=ra (a is the angle measured 
from the first radial line), we can find the magni- 
tudes of b and y; and, consequently, we can deter- 
mine the outer and inner boundaries of the boundary 
layer of the jet. 


A comparison of the experimental and theoretical jet boundaries is shown in Fig. 15. 
This comparison shows that the agreement between the theoretical curves and experiment is 


satisfactory. 


Il. Quantitative Investigations of Jets Emerging from 


an Axially Symmetric Nozzle with Underexpansion 


1. The Geometry of a Jet Emerging from an Axially Symmetric 
Underexpanded Nozzle. Inthe usual free turbulent jets emerging from a nozzle under 
similarity conditions, the pressure in the jet remains almost constant, i.e., the jet widens only 
on account of the development of the boundary layer, whereas in a jet emerging from an under- 
expanded nozzle, pressure is higher than in the surrounding medium. If the degree of nonsimi- 
larity is not high and lies within the limits 0.8<n<1.2,then, as was pointed out above, the influ- 
ence of the pressure difference on the development of the jet is negligibly small and the jet 
Widens in the first place on account of the increaSe in the thickness of the boundary layer. In 
this case, the boundary of the jet is almost linear as in the case of usual turbulent free jets. 
When the degree of nonsimilarity is high (n> 1.2), the pressure in the first barrel is 
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much higher than that in the surrounding medium, so that the jet widens greatly away from the 
nozzle orifice; the increase in width on account of the broadening of the boundary layer is sec- 
ondary inthis case. The inner and outer boundaries are highly curved near the nozzle orifice. 


A schematic diagram of a submerged jet issuing from an axially symmetric underexpanded 
nozzle is shown in Fig. 16 for high values of n. 


The following notation has been used in Fig. 16: u is the local longitudinal velocity in the 
jet; Um is the maximum longitudinal velocity in a given cross section; uj, is the velocity on the 
jet axis; Rp is the radius of the outer boundary of the jet; Ry, is the radius of the line of maxi- 
mum velocity, while in the transition region R,, denotes the radius of the inner boundary of the 
boundary layer; and Ry is the radius of the core in which the velocity is a constant and equal to 
Up: 


For convenience, it is advantageous to subdivide the jet into four regions: 
1) the initial region whose length is w; 
2) the first transition region in which a constant-velocity core still exists; 


3) the second transition region in which the axial velocity increases in the direction of 
the jet, but the maximum velocity in any cross section through the jet does not occur on the 
axis; 


4) the main region in which the axial velocity becomes the maximum velocity in a given 
cross section. 


In the boundary layer, the velocity increases from the outer boundary to the inner one and 
it can be subdivided into two parts: supersonic and subsonic. Beyond the forward discontinuity 
the flow velocity is subsonic; the oblique-discontinuity surface makes an appearance beyond the 
point of intersection of the discontinuities, because gas above point A passes through two oblique 
discontinuities, whereas below A it passes through one direct discontinuity, so that, consequently, 
the streams have different velocities. The gas which passes through two oblique discontinuities 
above point A has a supersonic velocity. Therefore, point A of intersection of the two density 
discontinuities marks the beginning of an annular zone in which intermixing of the subsonic and 
supersonic streams takes place. From the shadow photograph shown in Fig. 1, we see that the 
mixing zone developes rapidly. The tangential velocity discontinuity delineating the separation 
between the subsonic and supersonic flows becomes rapidly smoothed out. At high values of n, 
the diameter of the forward discontinuity is sufficiently large and the boundary layer at the end 
of the first barrel of the jet is strongly developed. Therefore, at the end of the first region, 
the zone between the outer boundary of the boundary layer and the region in which the super- 
sonic and subsonic streams intermix becomes very narrow. In this zone, the velocity decreases 
from its value on the inner boundary of the boundary layer in both directions — to its value in 
the subsonic core and to zero on the outer boundary. 


Initial lst transition, 2nd transition )Main region After the initial region, the jet in 
9 general will have a curvilinear boundary, 
Fs but its curvature is small. 
& 
N 
8 2. Basic Assumptions. The 


following assumptions were adopted for the 
derivation of the equations to be used for 
the calculation of the velocities in a jet 
emerging from an axially symmetric under- 
Fig. 16. Schematic diagram of the jet. expanded nozzle: 
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1) for n>2, it was assumed that the static pressure in the jet after the first barrel is 
equal to the pressure py in the surrounding medium; 


2) the velocity of flow behind the central forward discontinuity is constant and equal to the 
axial velocity; 


3) the constant-velocity core at the end of the initial region has the same diameter as the 
zone of the central forward discontinuity, while the velocity is equal to the velocity behind the 
forward discontinuity; 


4) after the initial region, the inner boundary coincides with the line of maximum velocity; 
5) the transverse velocity distributions in the zone situated between the outer boundary of 


the jet and the maximum-velocity line, and the zone lying between the maximum-velocity line 
and the constant-velocity core can be described by Schlichting's formulas, 


ay (1 — £15)?, (2.1) 


where 
2 Km for R,<R<R, (2.2) 
Rm — Ro 
and 
—- =(1— $15) (2.3) 
Urn 
where 


for Rm =R=R,- 


6) in the boundary layer of the jet, i.e., in the zone situated between the jet outer boundary 
and the maximum-velocity line, the variation in the temperature of the retarded stream can be 
described by 


7T,—T° 
75 — Tx 
where 
ya Kb (2.6) 
Rn — Py 


Ty* is the temperature of the retarded stream in the receiver, Ty, is the temperature of the 
retarded stream in the surrounding medium, and T* is the local temperature of the retarded 
stream in the jet. 


From Eas. (2.4) and (2.6) we have 


y==1—5. (2.7) 


With the exception of the outer boundary layer (R>Rym), the temperature in the jet is 
everywhere the same and equal to T)*. 


We can make the following remarks concerning these assumptions. As shown by experi- 
ments, the static pressure beyond the initial region of the jet deviates slightly from that in the 
Surrounding medium for all operating regimes of the nozzle except the regime n~ 2, but this 
deviation does not exceed 20%; pressure oscillations along the jet axis are larger than those 
along any other line parallel to the axis, but they are quickly damped out; for degrees of non- 
similarity lying in the range 0.8=n<=1.2, flow in the jet is almost indistinguishable from simi- 
larity flow. Therefore, the first assumption will not introduce large errors. 


The final cross section of the initial region of the jet is situated close to the position of 
the central forward discontinuity, so that prior to the section 0—0 the zone of intermixing of the 
supersonic and subsonic flows which begin at the point A has not yet developed. The thickness 
of the mixing Zone in the cross section0—0Ois very small by comparison with the diameter of the 
constant-velocity core, so that the difference between the constant-velocity core and the diam- 
eter of the central discontinuity can be neglected. 


There is no boundary layer in a jet of perfect gas. The results of calculations obtained by 
the method of characteristics show that the maximum velocity in the initial region occurs on the 
jet boundary. A boundary layer is formed in a jet of real gas; it is clear that until the second 
oblique discontinuity, the maximum velocity is to be found on the inner boundary of the boundary 
layer. The second oblique discontinuity is of variable intensity. Graphs of density distributions 
obtained by means of an interferometer clearly show that the intensity of the second oblique 
discontinuity increases as point A is approached, while the boundary tothe supersonic part of the 
jet lies inside the boundary layer. On the basis of these arguments, we can conclude that be- 
yond the initial region, the velocity maximum lies on the inner boundary of the boundary layer. 


It is known from experiment that in both the initial and main regions of the jet the longi- 
tudinal velocity profile across the jet is well approximated by Schlichting's formula. As has 
already been pointed out, the zone between the constant-velocity core and the boundary layer of 
the jet is relatively narrow, so that the approximation of the true velocity profile by Schlichting's 
formula should not give rise to a large error. 


The transition region of the jet is similar to the initial region of the usual free turbulent 
jet where the temperature of the retarded stream is constant in the boundary layer. There is 
no reason for not accepting this assumption in the case under consideration. 


3. The First Transition Region. In order to determine the maximum velocity 


in any cross section of the jet, we will make use of the equation of conservation of momentum 
(between the nozzle orifice and an arbitrary cross section of the jet) which can be written as 


Rb Ro 
TRexitl Perit Ue it + Bexitl == 2t \ pu°RaAR + Qn ( pRdR--= (Ry —Rexit) Pw (2.8) 


i} 
where pis the local density in the jet, pa,i+ is the air density at the nozzle orifice, u is the local 
velocity, and Uexit is the velocity at the nozzle orifice. 


On the basis of the first assumption, we can conclude that beyond the initial region the 
local pressure p is equal to the pressure of the surrounding medium py, 1.e., p=py- 


Equation (2.8) can then be simplified to 


Rb 


Re it exit Mexit + Pexit — Pu) = 2 { ow’RAR. (2.9) 
6 
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Let R denote a dimensionless radius, so that 


R=, Rae, R,=Ae Ra 


’ ) mo , 0 . 
Rexit Rexit Rexit Rexit 


Equation (2.9) can be rewritten as 


Rb 
2 — 
exittexitl Pexit— Pye \ pu-RdaR 


or 


k—} 


Pexitl exit — 9) 7 Pall RdR 
+ Pe xit— P= ’ 
Esp ( b \rexi 


0 


where J is the mechanical equivalent of heat (J =427 kg/kcal), and g is the acceleration due to 
gravity. 


After some transformation, we obtain 


Rb 
R 9 n—l ue _— — 
aa ee es —— ad 
’ eae i | \ 2gJe,( T?—-—~ =) . (2.10) 
6 co 2eICp Rk 
where 
4 e 
nae Lckit, 
Pu 


Since the first initial region contains a central core with constant velocity uy) and radius 
Ro, Eq. (2.10) will take the form 


— 


Ro, Rp 
— | aMexit-+ A — Rome] = \ 7 RAR + | oo Rdk, 
QeIcp | T* — — 2QeJcp | T* — — . 
R, . co sales) | k Ran 6 co ( sein) | R ) (2.11) 


where M, is the Mach number in the constant-velocity core and on the basis of the second as- 
sumption it is equal to the Mach number beyond the forward discontinuity, namely, Mg. Sub- 
stituting (2.2), (2.4), (2.5), (2.6) into (2.11), we obtain 


] 


= | aMexit +. — ] —RiMo|= —___+*—_~ (R,—* (R» — R,)] (R, —R,) ait 
2eJcpT, [1 — 
oe ( ogiesTs) 
0 
1 
us _ _ __ _ _ 
Fe Rn +8 Ry — Ru] Ro — Re) 
2EI5cp T 4 | + (22 — i) oy | 
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or 


1 


amt —1 ++} [Ry 3 Rn Ro) (Rn — Ry) 42+ 
l— 


] y) 
[eMexact 


Q 


4 — | +}——_____1______][R, +:(R, —R,)|(R, —R,) a, (2.12) 


u- 
2gJepTo 


5 (1+ 4(0—1)] 


where 6=Tj /Ty; n=1-£. Using the abbreviations 


A=—!] + n ) (2.13) 
}.-—— 
2gIcpT 
B=—1i+ _ (2.14) 
u- 
I-— 
IeJcpT. 
——~ (l+5@—)} 
and substituting (2.13) and (2.14) into (2.12), we obtain 
_ oe 1 
—S | PMexie tS — ROMS] = Rp (Rn— Ro) \ Ade — 
0 
_ _ 1 1 
= (Rn Ro) \ ABE + Rm (Ry — Rn) (ane —R,,)° \ Bide. (2.15) 


0 0 


The expressions for A and B are very complicated and, therefore, the evaluation of the 
integrals in (2.15) is difficult. In order to carry out the evaluations, we expand the functions A 
and B in Taylor's series. Formula (2.13) can be written as 


i 


_y uw (2.16) 
AS pet oh 


Since w?/2gJepTo is always less than unity, this series will be convergent. 


Let us now put expression (2.14) into the form 
- B12! 


pay 
(2gJepTG [1 + (1 —§) (0—1)]}! 


(2.17) 
i=1 


where 


l 
1+ (1— &)(8—1) 


= $* [(1 2) (@— 1. (2.18) 
jel 
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Fig. 17. Auxiliary function. 


Since 0 < ~é <1, the series (2.18) will only converge when @ < 2. The Series (2.18) is a 
component of series (2.17) and therefore series (2.17) will only converge when 6§<2. 


The value of function A is found by the substitution of expression (2.1) into Eq. (2.16). The 


integrals appearing in Eq. (2.15) can then be written as 


1] 
pe ROI Lf R-1 w2\y yy f RAT 8 \3 
| Adi= 5 — (Ads, 1 (Gop) Adee (4S i) (Aa i (2.19) 
) 
ye k—1 .9 k- 1 .9\?_, —f k--1 298 
Aid: = a 10(A,):, + a 15 | (A,):, i+ (4 - iS | (A): 1- (2.20) 


The values of the coefficients A,, Aj, A3, Ay, As, and A, are given in Figs. 17-22 as func- 


tions of M=um/Uy and &. 
Here é is the upper limit of the integrals fADé and fAédé, while A, is the velocity coef- 


ficient in the constant-velocity core 


Lg 
= —— 


‘0 , 


ter) PET) piers 


1 
| Bede can be evaluated by the substitution of Eqs. (2.3) and (2.18) 


1 
The integrals {Bae and 
into (2.17) 0 


0 


> » . 
) Bs +( a 3) Ba, (2.21) 


(2.22) 


eo k—1 .9 kR—1 .2\2 k— 1 3 
Bidi= 2B ( ) B. ( i) Be. 
\ 6 0 st k+) 0 s+ b+] 0 ¢G 


The values of the coefficients B,, By), B;, By, Bs, and Bg are shown in Figs. 23-28 as func- 
tions of m and 9g. Substituting (2.19), (2.20), (2.21), and (2.22) into (2.15), we obtain the equation 


of conversation of momentum in the first transition region of the jet as 


k I 9 n— | F242 = > = k—1l \2 k—1 .2\2 
Sy [Mei t AG RMB] = Ry Rn Ro) [FE ANE. (GSB (Ads ch 


—] .9\3 — — . _ 
ait | (Addt] = Rn Rol? | Fo (Adit + 
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Fig. 18. Auxiliary function. 


+ (ASB Yad, (FEE BY odes] + RalRy Red 


k+ 1 
,{ R—1 42 k--1 .9\2 kR—1 .9 
x| k+l 0B, +( k+l ‘0 | Bs + R41] is ) Ba | + 
DH HD k—1 \.2 ; k—] .9\- k—1] .2\8 
+ (A, —RnP| bal OB, +( rarer ‘| B,+( EI is | B,,. (2.23) 


In Eq. (2.23), the values of 6 and A, depend on the operating regime of the nozzle. The val- 
ue of @ is known for a given operating regime. The value of A, can be determined from the for- 
mula Ajg=Ac, Where Ag, the velocity coefficient behind the forward density discontinuity, is de- 
termined from the gas-dynamic calculation of the region of the first barrel. 


The value of M, is given by 


Formula (2.23), therefore, contains only four unknowns: m, Rm, Ry, and Rp. AS was 
mentioned above, the pressure in the jet beyond the initial region is almost equal to the pres- 
sure of the surrounding medium py. This means that, starting with the first transition region, 
the widening of the jet occurs only through the increase in the thickness of the boundary layer. 
Let us assume that the widening of the boundary layer in the first transition region follows the 
same law as that holding for the initial region of a usual free turbulent jet emerging from a 
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Fig. 19. Auxiliary function. 
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Fig. 20. Auxiliary function. 


OA, 


nozzle in the similarity regime. We can now write down the expression governing the increase 
in the thickness of the boundary layer of a submerged jet as 


aa 
Rb = Rm) 0.22 ta, (2.24) 


ax 2 


where py, is the density on the maximum-velocity line, and py is the density on the outer boun- 
dary of the jet. 


It is known from turbulence theory that the instantaneous velocities at any point of the jet 
can be decomposed into time-averaged and fluctuation components 


u=utu’, v<=THV’, 


where u is the longitudinal, i.e., the main velocity component, and v is the transverse velocity 
component. 


The rate at which the thickness of the boundary layer increases is proportional to the 
fluctuation component of transverse velocity v'. On the basis of experiments reported in [4], 
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Fig. 21. Auxiliary function. 
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Fig. 22. Auxiliary function. 


we know that the distribution of the dimensionless rms fluctuation component of transverse 
velocity, Vv'2/Um, and of the average velocity gradient 


oT (2.25) 
of 7 


across the jet are similar to each other. Therefore, the mixing length will also be approximate- 
ly constant across the jet: //b=const. In the zone situated between the outer jet boundary Rp 
and the maximum-velocity line Ry, and in the zone situated between the maximum-velocity line 
Rm and the boundary of the constant-velocity core Rg, the dimensionless velocity profiles are 
similar. Since the zone situated between the radii R,, and Ry lies immediately next to the zone 
situated between the radii Rp and Rm, we have grounds for assuming that the mixing lengths in 
thetwo zones will be equal. On the basis of the similarity of the velocity profiles in these zones 
and the equality of the mixing lengths, we can conclude that the fluctuation components of the 
transverse velocity will be the same and, consequently, the increase in the width of both zones 
will be subject to the same laws, which can be written as 


_d(Rm— Ro) (Rb = Rin) 
ax (LX 


(2.26) 
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Fig. 23. Auxiliary function. 
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Fig. 24. Auxiliary function. 


Moreover, it is known frorn experiment that the outer boundary of the jet remains linear 
and then we can write down the following equation for it 


aR —C __ bm 
ax b y , (2.27) 


where Cp is an experimentally determined coefficient Which for an air jet issuing from an 
axially symmetric underexpanded nozzle is equal to 0.117. 


Taking into account equalities (2.24), (2.26), and (2.27), we find after some transformation 
the following equations for determining the various boundaries: 


Pr 
14 
aR _ -o2a| Pm ) (2.28) 
ax 2 , 
ip py Se (2 .29) 
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Fig. 25. Auxiliary function. 
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where 
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Fig. 26. Auxiliary function. 
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Substituting Eq. (2.31) into (2.28), (2.29), and (2.30) and integrating, we obtain 


R= — 0.3234 + Roo: 
Rz= —0,103)-LR 
R, =0.1179+R 
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Fig. 27. Auxiliary function. 


13] 


In these equations PHP is a function of the parameters A), 9, and m and can be written 


(2.31) 


(2.32) 
(2.33) 
(2.34) 


(2.35) 
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Fig. 28. Auxiliary function. 


Rmo, Roo, and Rpp are the values of the radii in the sectionO- Oof the first transition region of 
the jet; and x, is the distance from the nozzle orifice to the sectionO—O, x)p=w 


After substituting expressions (2.32), (2.33), and (2.34) into (2.23) and performing the nec- 
essary transformations, we have 


k—1 yee 
{ ‘ 


where a= —(0,0227 (a, +6,) — 0.0484 (a, — b,) + 0.0261 ; 


Jaay + bY +, 
(2.36) 


2 2 ° 
b= (O.1I7 Rpg F-9.103R oq) A, — 0.44 (Ring — Rog) do + (0.323 Any —0.103R 6, + 


(k — 1) Roo 
k+] k—| , 
2 2 


—_— 


C=Ring (Aino 7 Roo) a; — (Ano 7 Rog)? Qs + Ring (Ayo — Ro) b, + 


k— | Reo 


+.0.44(Rpy — Roy) 6. —0.1615 


+ (Ryo Rmo)” b, + 


; k—| 2 k—1 .9\ 
FONG ~ 13 (Ade t+ (S=-% ) (Anke 


a= ——-13 (A): aaa 15 ) (A; det(F 1a) (Ad 


k+] 


—l1.,. 3 
i) Ba 


k— | 2 k— ] 9 - p— } 2 3 
= 2B. 2) B+ Ms B.. 
b. 0 +S | a | (<a 0) 


For the first transition region we have £,=1, while the value of 7 corresponds to a defi- 
nite value of m. From expression (2.36) we can obtain two solutions for y for each value of m. 
However, the value of zy in the initial cross section 0-0 of the first transition region is equal to 
zero. 


R+] 


It can be seen from (2.34) that y increases continuously along the jet. Therefore, we will 
choose the smallest of the two possible values of y. 


Substituting the value of 7 obtained from Eq. (2.36) into Eqs. (2.32), (2.33), and (2.34), we 
can calculate the values of Ry, Rm, and Rp corresponding to the given value of m. For m=m,, 
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we have 


Ro=0 (2.37) 
(m, is the value of m in the end section I—I of the first transition section). 


Next, we must determine the value of x corresponding to the given value of m, where m is 
an unknown function of x. It is very difficult to obtain an exact evaluation of integral (2.35), so 
that an approximate method must be used. 


With this aim, we subdivide the first transition region into several intervals for each of 
Which Eq. (2.35) becomes 


i+-6(1— — 16 m*)| dx + 
\ + 


i 
l k— | 2.9 
+ fise(ie 2m?) | dx = 
+ fies ( sh el 
*i-\ 
*; 
4 ] kR—1 ,2 4 
= 014 >[t+8(1- — nim) | dx (2.38) 


Here i is an index which depends on the number of intervals chosen. 


Assuming that the length of each interval is small, we find that the increase in the param- 
eter m in each interval must also be small. This allows us to replace m in the integrand of 
(2.38) by its average value (mj + mj_,;)/2. As a result of this, Fq. (2.38) assumes the simpler 
form 


k—| i-1_\"] 
nado > EEO TFS [ie MSPELY | ar 


or 


x= te) tg, (2.39) 


k— | m;-+-mj;_, \- 
| Oo) ] — dy — EL 
T k+1 (. 2 y | 


Thus, we can calculate the parameters of the jet in the first transition region from formu- 
las (2.28)-(2.35), (2.39), (2.1), (2.2), and (2.8). 


4. The Second Transition Region of the Jet. Inthe second transition re- 
gion, the velocity profile becomes more flattened. The velocity on the jet axis increases, while 
at the end of the second transition region the maximum velocity becomes equal to the axial 
velocity. The results obtained for the first region can be used for calculating the parameters of 
the second region. We can use an approximate method for the calculation of the variable veloc- 
ity along the jet axis. We will assume that the flow in the zone lying between the axis and the 
maximum-velocity line in the second region is produced by a fictitious plane-parallel stream 
which has an infinitely wide core with a constant velocity uy (See Fig. 16). According to this 
assumption, the boundary of constant velocity uy can be extended as shown in Fig. 16 and the fol- 
lowing relation holds between the rates of increase of the thicknesses of the zones (Rp — Ry) 
and (Rm — Ro): 


4 (Rp —Rm) _ 4 (Rm — Ro) 


AX aX 
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Consequently, as in the first transition region, the following relations hold: 


Ry =O0.11T+Ryy, 


(2.40) 
R= —0.1039+R,1, (2.41) 
Ro=—0.323 y, (2.42) 


where the subscript 1 is used to label the parameters in the initial section of the second tran- 
sition region of the jet (i.e., the section I-I). Let us write down the equation of conservation of 
momentum for this section, 


Rm Rb 
Pit Le it + Paying Py=2 \ or RdkR+2 \ ou-Rdk, (2.43) 
0 Rm 


which, after Eqs. (2.2) and (2.4) have been substituted into it, becomes 


Eo l 
We xit + Pexit— Py=2 \ pie [Rin — & (Rm — o)] (Rin— Ro) ad; +2 Olt” [Raw “t (Ry ” --k,,)| (R, —k,,) Or 
0 


©) 


Pain Me 


(2.44) 


where £, is the dimensionless ordinate of the point in the fictitious jet situated on the jet axis. 
In Eq. (2.44), the first integral has an variable upper limit &). 


It is known from formula (2.2) that 


(2.45) 


The magnitude of £) varies from zero to unity. In the initial cross section of the second 
transition region (in section I—I), we have R,=0 and ,=1. In the end section of this region (in 
section N-N), we have Ry»,=0 and, consequently, €,=0. 


As the result of some transformations, Eq. (2.44) becomes 


k—] 


” — | ~ > _ k— --] »2 2 k—1 .9\8 
(mMexie + SS) = Re, Rn Red Sa HOAs (maby) + (G28 | As (mie) + (G's | Asm’) - 
sor 


I 


~(R,—R,)? ie IGA, (mo) + (45 a) (7m,2,)-+ 
+(47) As (sie)|-+Rn (Rb — Ral 5B, (Bm) + (4? 3) Bs (0m) +-(<— 13) By (Om) |-+ 
+ (Ry —Rn}* | $-5Bs (Oy) + (F-26) Bs Orr) +(F=—-5 J’ Be Oy) I. (2.46) 


Formulas (2.46) and (2.23) have the same form, with the exception that in (2.46) &) varies 
between zero and unity, whereas £,=1 in Eq. (2.23). 


The values of the coefficients A,, Aj, Aj, Ay, As, and A, are given in Figs. 17-22 as func- 
tions of the parameters m and £5. 


The coefficients B,;, By, Bs, By, Bs have the same values as in the first region of the jet 
and are given in Figs. 23-28. 
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Substituting Eqs. (2.41) and (2.42) into (2.45), we obtain 


— Rm (i _— §) 
0.292, + 0.103 — (2.47) 


-( 


From the system of equations (2.40), (2.41), (2.42), and (2.47), we find that 


(»Mexiet — )= [ — 0.0227 (a, + 6,) — 0.0484 (a, — 6,)] 7 


+ [0117 R12, — 0.4R gy Qy + (0.323R ,1 —0.103Rp,) 6; + 


7 7 7 oe - 7 7 (2.48) 
+ 0.44 (Rot + R ni) 42] b+ Ran (Qy —_ A»)- LR ant (Roy 7 Rimi) b, “+ (Ro, _ Rimi) Ox», 
where the coefficients a@,, a, b,, b) are given by 
k— |} k— 9 \3 
a= Pal mo A, (so) E Pa 8) A, (1455) +H(45 hi | Ax(m5,), 
a= FM Asm) (Gaz) Asm) + [ Lay (mE) 
_ ._ (2.49) 
a 0B, (71,6) + ( bal =) B,, (17,9) +(4S a B,;(8,m), 
by = A — 168, (nn, + (GS Mi ) B,(m; \+(4 = is ) B, (m9). 


Having found from Eq. (2.47) the value of ~ for each value of £,) and substituting it into 
formula (2.48), we can find the value of m corresponding to the given £). With this value of m, 
we can then determine the corresponding value of x from the formula 


2 (UY; — Yin) 
k— |} m+ m1 \° 
1+6] 1— ——]{2,, ——_——_— 
+ k+1 (% 2 | 


where xj is the distance from the nozzle orifice to the cross section I-I. 


xX;= 


+ xj-1, 


Substituting the known values of m and é, into Eq. (2.1), we can determine the velocities 
Um and up, in each cross section of the second transition region of the jet. 


5. The Main Region of the Jet. In the main region of the jet, the maximum 
velocity lies on the jet axis. In this region, the assumption that the pressure in the jet is a 
constant is close to reality. The flow pattern in the transition region affects the flow in the 
main region. In the jet region under consideration, the boundary is curvilinear and at a suf- 
ficiently great distance from the nozzle orifice the boundary becomes linear with an inclination 
of 0.22. In the initial part of the main region, the inclination of the jet boundary varies from a 
value, equal to the inclination of the boundary in section N-N, to 0.22. According to the theory 
of a free turbulent jet, the jet boundary is given by the formula 


Pr 
1+ — 
akRb -.| Pm | 


ax 


which after some transformation becomes 


aRp _ ( _ ka! 12 ‘) 
7” co] 1 +6 l ba howe ||. (2.50) 
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For the initial part of the main region of the jet, we can find the value of the coefficient cy 
as the arithmetic mean of the experimental values of the coefficients for the transition and main 
regions, 

= O.M7 + 0.22 __ 9 1685, 
2 
The equation of conservation of momentum in the main region can be written as 


Rb 
= Py=2\  pu?RaR. (2.51) 
8) 


f ie L 
Yexit “exit | Pexit 


After a number of transformations, this equation becomes 
k 2 n— | 
f= [mMexi + “—) (2.52) 


Differentiating Eq. (2.51), we obtain 


_ _ Oy db, 2.53) 
dR, = ——_—? ( 
Rp 26. dm din 


Substituting expression (2.53) into (2.50), we can obtain the following formula for the de- 
termination of x: 


n fo dhe 


Xj; = Va dm +X, 
0.3376.{1 4.01 —-2=* 22m: (2.54) 
mn. OL | + — kal Loin ] 
where Xn is the distance between the nozzle orifice and the section N-N, 
Woe Bi Be (AL GP a (Aa GY oe 
dm k+ ] Y dm hel adm a } dm” 
B, ~ . ~ 
a — 20 (0.1665 — 0.17825 -+-0.111492—0,02356") m: 
aii 
1B. » ~ ~ » ~~- . . 
os 40? (0.2053 — 0.54596 + 0.92316? — 0.72426 -+ 0.35530! — 0.10179 + 0.01299") m?: 
dB, 


ne 69°(0.7248 — 3.024 -+ 7.13166? — 10.143893 +- 
Qi 


+ 9.95889! — 6.81969 +. 3.339265 — 1.079467 + 0.216% — 0.01869") m*, 


For a jet emerging from a nozzle with 9=1, we have 


<= 0.1524m; “2. 0.4992m", “Bs 1,698m' 


dim din dim 


It is very difficult to evaluate the integral on the right-hand side of (2.54) analytically, so 
that a graphical method is recommended for the solution of 


dbs 
adm 


a ne (2.55) 
Pal + of kal ine) | 


¢(m)= 


where F is the area underneath the curve g(m) between the ordinates m=m, and m=m 
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On the basis of (2.54), we can write 


n° 
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TURBULENT SUBMERGED JETS OF REAL GASES? 


G. N. Abramovich, V. I. Bakulev, 
I. S. Makarov, and B. G. Khudenko 


The experimental investigations of axially symmetric jets of liquid nitrogen at super- 
critical pressures were carried out by means of special apparatus, a schematic drawing of 
Which is shown in Fig. 1. The principal part of the apparatus is the high-pressure chamber 1. 
Its operative section was constructed in two ways. The first variant incorporated a viewing port 
and was used for visual observations and optical measurements by means of an IAB-451 instru- 
ment, while the second variant was used for the measurement of velocity-head and temperature 
profiles across the jet. The working part of the chamber had an internal diameter of 116 mm. 


The present investigation concerned jets of liquid nitrogen (critical pressure p, =34.6 
atm, critical temperature T, =126°K) issuing into the pressure chamber from a shaped nozzle 
2 into gaseous nitrogen. Liquid nitrogen was at a temperature of 80-90°K, heated nitrogen at 
a temperature of 250-420°K. The pressure in the pressure chamber p,; was kept approximately 
constant at 40 atm. [Liquid nitrogen was ejected by compressed air from a reservoir with a 
volume of 0.4 m®. Twosystems were used for the nitrogen feed. Cold nitrogen was fed directly 
to the interchangeable nozzle of the apparatus through the butterfly valve 3. Nozzles with diam- 
eters of 5.0, 3.0, and 1.12 mm were used during the experiments. In the other system, the cold 
nitrogen passed through the butterfly valve 4 into the heat exchanger 5 which was heated by 
combustion chamber 6 and then the heated nitrogen was supplied to the assembly. An adjustable 

nozZle 7 was used to maintain constant pressure 
in the assembly. 


a The following quantities were measured 
during the experiments: T,, the temperature of 
the cold nitrogen by means of thermocouple 8 
situated at the entry tothe nozzle, T., thetempera- 
ture of the heated nitrogen by means of thermo- 
couple 9 at the entry into the working part of the 
pressure vessel, P,, the pressure in this section 
of the pressure vessel by means of manometer 
10, Aps, the pressure drop in the section between 
the nozzle and the working part of the pressure 


T The Russian text of this article was published 
in Izvestiya Akad. Nauk SSSR, Mekhanika 
Fig. 1 Zhidkosti i gaza, No. 1, 1966, pp. 154-158. 
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vessel by means of the differential manometer 11, hyq, the velocity head at one point of the hot- 
nitrogen line by means of the differential manometer 12. The velocity head and temperature 
across the jet were measured by means of a Pitot tube combined with a thermocouple 13, at 
several fixed points along the jet. The output of this thermocouple was measured by a PP-1 
potentiometer. The output of the velocity-head tube was fed to two differential manometers: a 
mercury manometer 14 used for measuring small pressure drops and a low-sensitivity manom- 
eter 15 used for the measurement of high pressure drops. 


The width of the jet mixing zone was measured by means of the JAB-451 instrument with 
the first variant of the working part of the pressure chamber. Figure 2 shows shadow photo- 
graphs of a cold nitrogen jet and a gas jet with a very small temperature difference (an almost 
isothermal jet). The outer boundary of the jet is clearly visible in these photographs. The in- 


Fig. 2. Shadow photographs of (a) a cold nitrogen 
jet (Ty =80°K, T, =363°K) and (b) of an isothermal 
jet. 
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Fig. 3. The velocity-head and temperature 
profiles in the initial region of the jet. The 
experimental points (open symbols for 7,, 
filled in symbols for Q,) were obtained for 
the following values of the parameters 
(x/Yq, To°K, 6), Yo mm): 1) 10, 417, 20.8, 
2.5, 2) 16.7, 360, 18.8, 1.5, 3) 16.7, 378, 16, 
1.5, 4) 16.7, 410, 6.9, 1.5, 5) - —1.25; the 
theoretical curves are 1)Q,, a=0.0905, 

2) 71, @=0.098, 3) Q, a@=0.1. 


ner boundary of the cold nitrogen jet is not 
visible because of the large density change 
in the jet boundary layer. 


The second variant of the pressure 
chamber was used for the measurement of 
the velocity head and temperature in several 
transverse sections of the jet when the emis- 
sion velocity uy of the jet was varied in the 
interval 20-50 m/sec. The variations of 
Reynolds number across the nozzle exit were 
included in the range R=(1.7-5.8)-10°. The 
velocity of the parallel flow did not exceed 
0.3 m/sec, so that the jet could be considered 
to be submerged in all cases. The experi- 
ments were carried out for two temperatures 
of the heated nitrogen: Ty» =(250-300)°K and 
Ts = (870-42 0)°K. 


In the case of the initial region, the 
experimental data were processed with the 
help of the relations 


_ <puy uy (2) (AT) — <T)—Ty = ;(¥) 
— — —_ ; — a rr 


Pouo” x 


and in the case of the main region with the help of the relations 


{pu> <u> 
Pgs? 


Q.= 


<AT>) — <T)—T? = ;(4). 


= (4) and t= 


AT T3;—Ts 


Here x, y are the coordinates of the measurement point, yo” is the velocity head at the 
nozZle exit, p33” the velocity head on the jet axis, <pu> <u> the averaged velocity head at the 
point of measurement, <T> the average temperature at the point of measurement, p, the nitro- 
gen density at the nozzle exit, and p. the density of the surrounding medium. 


The variation of the velocity head along the jet axis has been described as follows: 


2 
P3U3 
3=-— 
Palo 


where ry, is the radius of the nozzle. 


2) 


The velocity-head and temperature profiles in the initial region are shown in Fig. 3. The 
same figure also shows for comparison purposes the velocity-head profile for an isothermal air 


jet obtained with other apparatus. 


It can be seen from a comparison of the experimental data for the initial region of the 
jet that the boundaries of a low-temperature nitrogen jet in this region are linear as in the case 
of the isothermal jet. The width of the boundary layer of the cold nitrogen jet is less than that 
of the isothermal jet and, consequently, the cold-nitrogen jet has the higher range. The open- 
ing angle of the jet is independent of the pressure drop Aps at the nozzle (emission velocity up). 
The profiles in the initial region of the low-temperature jet were found to be affine and fuller 


than those of the isothermal jet. 


142 G. N. ABRAMOVICH, V. I. BAKULEV, I. S. MAKAROV, AND B. G. KHUDENKO 


Fig. 4. The velocity-head and temperature profiles in 
the main region of the jet. The experimental points 
(open symbols for 7,, filled in symbols for Q;) were ob- 
tained for the following values of the parameters (x/rp, 
To°K, 69, r9 mm): 1) 49.4, 250, 8.09, 1.5, 2) 49.4, 230, 
7.37, 1.5, 3) 44.6, 282, 7.84, 0.56, 4) 34.7, 287, 10.4, 
1.5, 5) 34.7, 404, 15.1, 1.5, 6) — — 1, 1.5; the theoreti- 
cal curves are: 1)Q,, a=0.07,2)) T2, a=0.0665, 3) Q., 
a=0.08. 


On the basis of experimental data (Figs. 2, 3), we can say that at supercritical pressures 
the liquid nitrogen jet qualitatively behaves in the same manner as the usual gas jet [1]. In 
other words, in the case of a liquid nitrogen jet (temperature below the critical temperature) 
propagating in a medium consisting of gaseous nitrogen (temperature above the critical temper- 
ature), mixing in the nonisothermal boundary layer at supercritical pressures occurs without 
the formation of droplets. 


} & The figures showing the experimental 
b ‘a 42 e velocity-head and temperature profiles for 
rn aa a 42 the initial region of the jet also show the cor- 
4 Ne} ,2 | = 09 responding theoretical profiles calculated 
os ee Ser from the equations given in [2]. 
\ 
. The theoretical profiles were fitted to 
a\e —_—_ | the experimental points with the help of the 
, oy NN | experimental constant a. The turbulent 
ae . l o SS | Prandtl] number was chosen such that the 
a | | ep ANN values of the coefficient a for the dynamic 
0 004 008 012 y and thermal profiles were as close as pos- 
r sible to each other. It was assumed for the 
Fig. 5. The velocity-head and temperature isothermal jet that a=0.1. In the case of the 
profiles in the main region of the jet. The cold nitrogen jet with Prandtl number P =0.87, 
experimental points (open symbols for 75, the value of the coefficient a for the velocity- 
filled in symbols for Q.) were obtained for head profile (a=0.0905) was only a little dif- 
the following values of the parameters ferent from its value for the temperature 
(x/Y9, To°K, 6, Yr mm): 1) 132, 389, 1.97, profile (a=0.093). With these values of the 
0.56, 2) 92.8, 408, 3.2, 0.56, 3) 66.7, 200, coefficients, the theoretical profiles were in 
2.71, 1.5; the theoretical curves are: 1) satisfactory agreement with the experimental 


Qo, @=0.067, 2) To, a=0.064. points. 
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The values of the experimental coefficient a 
obtained above indicate that its magnitude does not 
remain constant, but varies somewhat with a suf- 
ficiently large change in the relative density 9, = 
Po/P2, Which in our case changed from 9, =1 for the 
isothermal jet to @) =20 for the cold nitrogen jet. 


The velocity-head and temperature profiles 
obtained for the main region of the jet are shown 
in Figs. 4 and 5. The velocity-head profile for an 
isothermal jet is given in Fig. 4 for comparison 
purposes. 


As can be seen from these data, the velocity- 
head profiles for the main region of the jet are 


Fig. 6. Variation of the velocity head 


along the jet. The experimental points not affine. 

are: 1) p./py=1, T, =280°K, 2) 0.061, According to the experimental data, the main 
280°K, 3) 0.039, 400°K; the theoretical region of the cold nitrogen jet can be roughly sub- 
curves are: 1) p)./p,=0.061, a=0.072, divided into two subregions: the first region has a 
2) 0.039, 0.072, 3) 1, 0.083. density ratio §,=15-8 (at the beginning of the main 


region; See Fig. 4) and the second has a density 
ratio 6,=4-2 (at a greater distance from the nozzle; see Fig. 5). The profiles for the first 
regions are fuller by comparison with those for the second region. This can be seen most 
clearly from the temperature profiles. 


The relative width of the main region of the cold nitrogen jet is less than that of the iso- 
thermal jet. 


On the basis of the experiments performed, we can say that the profiles change along the 
jet and at a sufficiently large distance from the nozzle they approach the profiles for an iso- 
thermal jet. The relative width increases with distance away from the nozzle and also ap- 
proaches its value for an isothermal jet. 


Thus, the experimental data confirm the assumption normally made for theoretical calcu- 
lations that the profiles vary with position along a jet of real gas [2, 3]. 


The velocity-head and temperature profiles for the main region of an axially symmetric 
cold nitrogen jet were calculated on the basis of [3]. The comparison of the theoretical pro- 
files with the experimental ones has been made with the help of the experimental coefficient a 
in the same manner as was done for the initial region. For a jet with Prandtl number P =0.87, 
the values of the coefficient a for dynamic profiles (a=0.07-0.067) and for thermal profiles 
(a=0.0665-0.064) are found to be approximately the same. The theoretical profiles obtained 
with the given value of a are in satisfactory agreement with experimental data. The value of a 
chosen for the isothermal jet was a=0.08. As in the initial region, the coefficient ais not a 
constant but varies within relatively narrow limits when there is a significant change in the 
initial density ratio. This variation in the value of the coefficient @ as a function of the initial 
temperature ratio T,/T, can be explained if it is assumed that the approximate theory given in 
(2, 3] does not take account of all factors influencing the deformation of profiles and the position 
of the jet boundary. 


The variation of the relative velocity-head along the main region of a cold nitrogen jet is 
shown in Fig. 6 for several temperatures of the surrounding medium, T,, together with the 
relative velocity head in isothermal jets. As can be seen from this figure, the range of the 
cold-nitrogen jet is greater than that of the isothermal jet. The smaller the temperature ratio 
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T/T», the smaller the opening angle of the jet and, therefore, the greater the range of the jet. 
The opening angle is independent of the emission velocity up. 


The comparison of the experimental data with the results of theoretical calculations car- 
ried out on the basis of [3] shows that the values of the experimental coefficient a for the veloc- 
ity head along the jet axis are little different from the values of the corresponding coefficient 
obtained from measurements across the jet. The theoretical variation of the velocity head along 
the jet axis is in satisfactory agreement with experimental data. 
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